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Simulating  the  Effects  of  Laser  Damage  to  the 

Retina 

A  Detailed  Review 


1.  Introduction 

The  Air  Force’s  interest  in  eye  damage  caused  by  exposure  to  intense  light  began  with  the 
dawn  of  the  atomic  bomb.  Atomic  bomb  flashes  emitted  a  bright  light  Aat  caused  visual 
damage  to  observers.  Presently,  eye  damage  due  to  intense  light  exposure  is  more  likely 
to  be  caused  by  lasers  than  by  atomic  bombs.  This  is  because  lasers  have  become 
prevalent  in  the  military;  they  are  used  for  both  aiming  devices  and  weapons.  Both 
accidents  and  direct  enemy  action  result  in  pilots  being  particularly  vulnerable  to  eye 
damage  from  accidents  and  direct  enemy  action. 

Laser  damage  to  the  eye  can  result  in  a  number  of  injuries.  Hemorrhages  can  occur  if 
arteries  are  hit  by  lasers  and  even  blindness  can  occur  if  there  is  strong  laser  blow  to  the 
to  the  optic  disk  which  destroys  the  optic  nerve.  Since  the  topic  of  laser  damage  is  so 
vast,  this  Phase  II  effort  focused  on  retinal  receptor  drop  out  due  to  heat  or  occlusion  by 

blood.  It  was  assumed  that  no  underlying  nerves  were  affected  and  that  retinal  loss  was 
“all  or  nothing.” 

The  present  Phase  n  combined  Air  Force,  university  and  CMI  to  research  the  effects  of 
laser  eye  damage  on  pilots’  performance.  The  research  produced  exciting  experimental 
results,  innovative  models  (biologically  plausible,  visualization  and  signal  processing) 
and  two  commercial  products.  The  experimental  results  came  from  two  sets  of 
psychophysical  experiments  performed  at  Conceptual  MindWorks,  Inc.  (CMI)  and  the 
University  of  Texas  at  San  Antonio  (UTSA),  VHP  experiments  performed  at  the 
Umversity  of  Texas  Health  Science  Center  at  San  Antonio  (UTHSC)  and  biological 
experiments  performed  at  UTSA.  The  biologically  plausible  model  resulted  from  a 
collaboration  between  Dr.  Wilson  Geisler  at  the  University  of  Texas  at  Austin  (UTA)  and 
CMI.  The  visualization  model  was  a  result  of  a  collaboration  between  Brooks  Air  Force 
Base  and  Conceptual  MindWorks,  Inc.,  and  the  signal  processing  model  resulted  from 
work  performed  at  UTSA.  The  combination  of  data  collection,  modeling  and  commercial 
products  has  resulted  in  the  following  Phase  II  report. 


2.  Evolution  of  Phase  II 

Phase  I  resulted  in  a  model  that  demonstrated  the  decrease  in  acuity  due  to  a  laser  lesion 
in  the  foveal  region.  The  model  was  called  the  “Geisler”  VEP  model  and  was  based  on 
Visual  Evoked  Potential  (VEP)  data  from  rhesus  monkeys  with  laser  lesioned  eyes. 
Phase  n  proposed  extensions  of  the  “Geisler”  model  and  a  new  signal  processing  model 
from  Phase  H.  The  proposed  extensions  included  incorporating  VEP  data  from  humans 
with  normal  vision,  VEP  data  from  humans  with  eye  diseases  whose  effects  were  similar 
to  those  of  laser  damage  and  VEP  data  from  humans  whose  vision  was  altered  by 
anificial  scotomas.  Another  proposed  extension  was  to  develop  a  visualization  tool 
programmed  in  C  that  allowed  pilots  and  vision  researchers  to  view  what  a  scene  would 
look  like  if  he/she  had  incurred  laser  eye  damage.  This  tool  was  also  used  to  produce  a 
commercial  visualization  product  that  allowed  patients  and  their  families  to  view  how  the 
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diseased  eye  viewed  a  scene.  The  proposed  signal  processing  model  was  to  investigate 
the  modeling  of  retinal  neurons  and  their  signal  processing  properties.  This  was  to  result 
in  a  better  understanding  of  how  damage  to  specific  photoreceptors  in  the  retina  has 
secondary  effects  on  the  surrounding  sensory  outputs. 

Over  the  course  of  the  Phase  II  project,  a  number  of  changes  were  made  to  the  original 
proposal.  These  changes  were  due  to  advancements  in  knowledge  and  technology  as  well 
as  the  changing  needs  of  the  Air  Force.  All  of  the  changes  were  made  with  approval  and 
feedback  from  the  Air  Force  contract  monitor.  First,  the  “Geisler”  VHP  model  was 
replaced  by  the  Inhomogenous  Retino-Cortical  (IRC)  model. .  Psychophysical 
experiments  were  collected  for  the  Inhomogenous  Retino-Cortical  models  (IRC).  VHP 
data  of  diseased  and  normal  subjects  were  collected  in  addition  to  the  psychophysical  data 
collection.  See  Appendix  IV  for  the  results  and  conclusions.  Second,  the  visualization 
model  (the  AVS/Concentric  ring  model)  was  implemented  in  AVS  so  that  the  tool  could 
be  incorporated  into  a  larger  set  of  Air  Force  tools  called  ILPEM.  Also  three  of  the  four 
research  projects  at  UTSA  were  oriented  to  aid  in  the  new  model  development.  Third, 
the  Inhomogenous  Retino-Cortical  (IRC)  model  was  used  to  estimate  the  detectability  of 
objects  such  as  airplanes  in  the  sky. 

Other  deviations  from  the  Phase  11  proposal  were  in  the  form  of  additions.  A  second 
commercial  product,  called  the  VisionProbe,  was  added  to  the  Phase  II  work. 


3.  Psychophysical  Data  Collection  for  Model 

The  psychophysical  experiments  were  designed  to  generate  a  complete  set  of  contrast 
sensitivity  functions  (CSF)  and  to  test  performance  effects  of  laser  eye  damage.  Two 
sets  of  stimuli  were  used,  the  first  was  a  set  of  Gabor  patches  and  the  second  was  an 
airplane  picture  that  was  perturbed  by  different  artificial  scotoma.  Both  sets  o  stimuli 
were  presented  at  various  eccentricities,  the  goal  of  which  was  to  determine  the  visual 
capabilities  at  various  eccentricities.  The  Gabor  patches  were  also  used  to  measure  the 
spatial  summation  of  the  third  harmonic  of  gratings  with  a  without  the  presence  of  a 
fundamental  frequency.  A  more  detailed  description  of  the  psychophysical  experiments 
can  be  found  in  Appendix  II.  The  results  of  the  psychophysical  experiments  were  used  to 
calibrate  both  the  AVS/Concentric  Ring  Model  and  the  IRC.  Results  from  the  IRC  and 
the  psychophysical  experiments  were  then  compared  in  order  to  validate  the  IRC  as  a 
model  of  visual  performance. 


4.  Model  Development 

4.1  VEP  Model 

The  visual  damage  modeling  effort  began  in  Phase  I  with  the  “Geisler”  VEP  model.  The 
model  was  based  on  VEP  data  collected  from  rhesus  monkey’s  who  had  their  eyes 
damaged  by  controlled  laser  blasts.  VEP  measurements  were  taken  before  and  after  the 
blasts.  The  model  basically  divided  the  retina  into  a  series  of  imaginary  concentric  rings. 
Each  ring  is  assigned  a  mathematically  derived  weight  value  based  on  ganglion  cell 
densities.  Initially,  all  the  weights  were  set  to  one.  When  a  pattern,  such  as  a  sine  wave 
pattern,  is  passed  through  the  model,  the  first  operation  is  to  calculate  the  pattern’s 
contrast  sensitivity  for  each  ring.  Then,  the  responses  for  all  the  rings  are  summed.  Since 
there  is  considerable  variation  between  one  subject’s  and  another’s  VEP  data,  the 
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summed  ring  responses  are  multiplied  by  a  subject  weight.  The  model  is  trained  for  each 
subject  using  a  minimization  of  least  squares  algorithrn.  Please  refer  to  Appendix  II  for 
the  results  and  conclusions  of  the  Phase  I  effort. 

4.2  A  VS  Concentric  Ring  Model 

In  Phase  II,  a  visualization  tool  was  developed  to  display  the  maximal  amount  of 
information  available  at  each  point  in  a  receptive  field  for  a  person  with  a  laser  eye  injury 
or  a  disease  that  causes  photoreceptor  drop  out.  To  display  the  image,  the 
A VS/Con centric  Ring  Model  was  used.  For  displaying  a  scene  as  an  undamaged  eye 
would  see  it,  the  user  must  first  choose  a  foveation  point  in  the  image.  Then,  like  the 
Geisler  VEP  model,  a  series  of  rings  representing  different  eccentricities  from  the 
foveation  point  are  superimposed  on  the  scene.  A  contrast  sensitivity  function  is  defined 
for  each  ring  based  on  either  the  literature  or  the  data  collected  at  CMI  and  described  in 
the  Psychophysical  Experiments  section.  A  low  pass  filter  that  corresponds  to  the 
contrast  sensitivity  function  is  applied  to  the  pixels  in  the  corresponding  region. 
Therefore,  the  blur  increases  ring  to  ring  as  the  image  moves  radically  outward  from  the 
foveation  point.  To  simulate  the  effect  of  laser  damage,  total  photoreceptor  dropout  is 
assumed.  To  simulate  this,  the  damaged  region  is  replaced  by  the  average  illumination  in 
the  total  image.  This  basic  technique  may  also  be  applied  to  multiple  foveation  points. 
For  a  more  detailed  description  of  the  model,  please  refer  to  Appendix  I.  A  user’s  guide 
is  presented  in  Appendix  V. 

4.3  IRC  Model 

The  Inhomogenous  Retino-Cortical  (IRC)  model  is  based  on  Ideal  Observer  theory.  Ideal 
Observer  theory  takes  well-established  scientific  facts  about  vision  and  puts  them  into 
one  model.  The  facts  include  physical  properties  of  light,  the  optics  of  the  eye  and 
absoiption  of  spectra  of  pigments.  The  resulting  Ideal  Observer  depicts  what  the  physical 
properties  of  the  eye  are  capable  of  seeing  or  detecting  and  illustrates  where  information 
is  lost  in  the  visual  system.  A  d’  measure  is  used  to  quantify  the  ability  to  see  or  detect  an 
object.  The  model  does  not,  however,  do  a  good  job  of  modeling  what  is  actually  seen  by 
a  person.  This  is  thought  to  be  due  to  the  fact  that  a  human’s  vision  is  based  on  not  only 
the  physics  of  the  eye  but  also  cortical  functions.  In  layman’s  terms,  your  brains  are 
constantly  changing  and  adapting  the  information  our  eyes  send.  Therefore,  what  our 
eyes  see  and  what  we  perceive  are  two  very  different  things. 

The  IRC  takes  the  ideal  observer  theory  and  adds  floating  parameters  so  that  actual  visual 
performance  of  a  subject  can  be  modeled.  The  parameters  are  fit  using  a  set  of  contrast 
sensitivity  curves.  To  this  end,  we  collected  a  series  of  contrast  sensitivity  curves  using 
psychophysical  experimentation,  see  section  2.0  Psychophysical  Experiments.  Once  the 
parameters  have  been  fit,  the  IRC  model  can  be  tested  to  discover  how  well  it  modeled  a 
subject’s  ability  to  detect  an  object.  For  a  more  detailed  description  of  the  model,  please 
refer  to  Appendix  II.  A  user’s  guide  is  presented  in  Appendix  VI. 

5.  Commercial  Products 

5.1  I'lsionProhe 
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5.1.1  Description 

The  VisionProbe  is  a  high-speed,  high-resolution,  computer-based  graphics  subsystem 
especially  configured  for  applications  in  vision  research  that  can  reside  on  any  AT 
compatible  computer.  It  is  capable  of  producing  color  and  monochrome  patterned 
stimuli.  A  separate  graphics  co-processor  is  used  to  produce  the  graphics.  It  resides  on  a 
separate  board  that  plugs  into  the  AT  compatible  bus.  This  is  an  advantage  because  once 
the  graphics  processor  has  been  instructed  on  the  type  of  stimulus  patterns  to  display,  the 
host  computer  can  devote  all  of  its  resources  to  other  tasks  such  as  data  acquisition  and 
analysis  while  the  graphics  board  produces  the  graphics.  Because  of  this  independence 
VisionProbe^^  integrates  smoothly  with  the  DataWave  Systems’  Experimenter’s 
WorkBench  so  that  each  product  can  perform  all  of  its  functions  at  full  capacity.  This 
enables  an  experimenter  to  collect  and  analyze  data  on  the  same  machine,  a  feature  that 
no  other  visual  stimulator  on  the  market  shares. 

5.1.2  System  Configuration 

VisionProbe  contains  the  following  components  when  assembled  in  standard 
configuration: 

•  Graphics  processor  board  (OmniComp  TEXAN  ET) 

•  I/O  adapter  board  (Conceptual  MindWorks  ViPSIA/1 ) 

•  High-resolution,  2 1 "  color  CRT 

•  STIM  software  which  generates  the  stimulus  patterns  and  controls  the 
graphics  processor 

•  ViP  software  which  is  the  command  level  interface  by  which  the  user 
communicates  with  the  graphics  board. 

•  Integrated  photometer  (UDT  Model  265  or  equivalent) 

The  graphics  engine  itself  is  the  Texas  Instmments  34020.  The  graphics  processor  board 
includes  2  Mb  of  video  RAM  for  holding  the  video  images  and  4  Mb  of  dynamic  RAM 
for  holding  the  STIM  software.  Video  output  to  the  CRT  is  accomplished  through  three 
8-bit  video  digital  to  analog  converters  (DACs).  This  permits  each  of  the  three  color  guns 
(red,  green,  and  blue)  in  the  display  to  be  set  to  one  out  of  a  maximum  of  256  luminance 
levels.  The  maximum  and  minimum  luminance  levels  are  set  by  the  operator.  The 
ViPSIA/1  adapter  board  provides  for  input  to  and  output  from  the  stimulator  system. 
These  signals  allow  for  synchronization  of  the  stimulator  with  external  data  acquisition 
and  display  devices. 

An  optional,  although  strongly  recommended,  component  of  the  VisionProbe  is  the 
mtegrated  photometer.  When  present,  the  photometer  permits  VisionProbe  to  perform 
automatic  gamma  correction  of  the  CRT  display.  Gamma  correction  is  necessary  to  make 
the  luminance  gain  function  of  the  display  linear  and  to  balance  the  output  of  the  three 
color  guns  so  that  colors  are  produced  correctly.  The  user  is  also  able  to  use  the 
photometer  to  verify  brightness  of  the  display.  If  an  external  photometer  is  available,  the 
display  brightness  can  be  measured  with  that  instrument,  and  gamma  correction  can  be 
accomplished  manually. 
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Figure  ] .  VisiouProhc  Standard  Configuration  Components 


5.1.3  VisionProbe  Development  Tasks 

5. 1.3.1  Background 

The  VisionProbe  was  originally  developed  in  the  late  1980s  by  Krug  Life  Sciences,  Inc. 
to  provide  a  source  of  visual  stimuli  for  vision  researchers.  The  product  was  sold  to 
several  vision  researchers  and  is  still  being  used  by  some  of  them.  This  product  was 
designed  based  upon  then-state-of-the-art  graphics  processors.  Since  that  time  the 
graphics  card  (National  Design,  Inc.  1288,  NDI  1288)  has  become  obsolete.  National 
Design,  Inc.  has  gone  out  of  business  and  the  CPU  that  the  software  was  designed  for  has 
been  enhanced  (The  TMS34010  has  been  replaced  by  the  TMS34020).  The  I^g  version 
of  the  VisionProbe  was  integrated  with  a  data  acquisition  system,  DataWave 
Technologies  (formerly  BrainWave,  Technologies)  Experimenter’s  WorkBench.  The 
unavailability  of  the  VisionProbe  has  left  the  Experimenter’s  WorkBench  without  an 
integrated  visual  stimuli  system.  Our  research  into  the  world  of  vision  research  has 
determined  that  there  is  still  a  need  for  a  product  similar  to  the  vision  probe,  especially  in 
situations  where  the  Experimenter’s  WorkBech  system  is  utilized  for  data  collection  in 
vision  experiments.  The  following  sections  discuss  the  tasks  that  were  performed  to 
return  the  VisionProbe  to  a  production  state. 

5. 1.3.2  Evaluation  Of  Design 

The  first  step  to  determine  if  the  VisionProbe  Design  could  be  produced  with 
commercially  available  parts  was  to  perform  an  evaluation  and  get  familiar  with  the 
existing  VisionProbe  Design.  The  graphics  cards  were  contacted.  OmniComp,  Inc 
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located  in  Houston,  TX  had  a  second  generation  product,  the  Texan  ET  which  would 
support  the  vision  probe  function. 


5.1. 3. 3  Existing  VisionProbe  Evaluation,  Assessment  And  Familiarization 

The  VisionProbe  source  code  and  users  manuals  were  reviewed  and  the  responsible 
engineers  became  familiar  wit  the  design  of  the  system.  Design  notes  from  the  latest 
versions  of  the  VisionProbe  were  compared  to  the  source  code  to  determine  the  last 
production  baseline  of  the  VisionProbe  product.  This  baseline  would  be  the  basis  for  the 
new  VisionProbe  product  to  be  developed.  Enhancements  and  added  capabilities  would 
be  considered  as  time  and  resources  permitted.  The  source  code  was  compiled  into  an 
executable  and  operated  on  the  hardware  to  ensure  that  all  software  components  were 
accounted  for. 

5.1.3. 4  VisionProbe  Update  Planning 

The  VisionProbe  can  be  divided  into  three  components  for  planning  purposes.  The 
software  which  was  originally  written  for  a  TMS34010  based  processor  required 
modifications  due  to  processor  differences.  A  large  portion  of  the  software  could  be 
transferred  to  the  TMS34020  system  with  only  small  revisions  for  Input/Output  and 
memory  usage.  The  graphics  board  was  purchased  from  a  vendor  and  configured 
properly  for  the  VisionProbe  application.  The  daughter  card  (ViPSIA)  was  redesigned 
based  upon  the  new  graphics  card  interfaces,  new  physical  envelope,  and  the  latest 
components  commercially  available  (some  of  the  ViPSIA  parts  have  become  obsolete 
since  the  original  design  was  completed). 

5.1.3. 5  Software  Development 

Software  development  for  the  VisionProbe  focused  on  the  conversion  of  the  software 
written  fro  the  NDI  1288  board,  TMS34010  processor  to  the  OmniComp  Texan  ET, 
TMS34020  processor.  The  software  required  revisions  for  the  following  functional 
areas: 


•  3401 0  versus  34020  differences. 

•  Graphics  Card  differences  such  as  VRAM  wrap,  bulk  clearing  of  VRAM,  etc. 

•  Palette  chip  differences 

•  ViPSIA  design  differences  including  addresses  and  provisions  for  expanded 
signals 

•  Host  Communication  Differences  (see  section  5. 1 .3. 13) 

5.1.3. 6  Hardware  Development 

Hardware  development  for  the  VisionProbe  consisted  of  producing  the  VisionProbe 
Stimulus  Interface  Adapter  (ViPSIA)  circuit  card.  The  block  diagram  for  the  ViPSIA  is 
shown  below. 
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Figure  2:  ]  iPSlA  Block  Diagram 


5.1.3. 7  Define  Interface  Between  The  Texan  ET  And  The  ViPSIA 

The  NDI  1288  based  VisionProbe  interfaced  with  the  ViPSIA  through  a  memory 
connector  provided  on  the  NDI  1288  board.  This  connector  was  not  available  on  the 
OmniComp  Texan  ET  board.  Several  alternatives  were  considered:  connect  through  the 
PCs  motherboard,  modify  the  Texan  ET  to  provide  data  and  address  connectors,  or  utilize 
the  floating  point  unit  (FPU)  socket  that  is  provided  on  each  Texan  ET  for  a  co-processor 
option.  The  FPU  socket  option  was  selected  due  to  its  suitability  to  the  design.  The  PC 
bus  option  was  ruled  out  due  to  the  overall  requirement  to  allow  the  VisionProbe  to  be 
co-located  with  the  data  acquisition  software.  The  FPU  socket  provided  all  the  necessary 
signals  required  fro  the  ViPSIA  with  the  exception  of  three.  These  signals  will  be  routed 
from  elsewhere  on  the  Texan  ET  to  FPU  socket  unused  pins,  power  and  ground  for  the 
ViPSIA  will  be  provided  from  the  FPU  socket  as  well. 

5.1.3.8  Define  Schematic  For  ViPSIA 

The  schematic  for  the  updated  ViPSIA  was  derived  from  the  Krug  ViPSIA.  The  basic 
desipi  for  the  analog  Input/Output  portion  of  the  card  remained  the  same.  The  digital 
portion  of  the  design  was  repackaged  into  a  Field  Programmable  Gate  Array  (FPGA)  for 
space,  power  consumption,  component  availability  and  flexibility  reasons.  The  modified 

schematic  was  developed  to  suppon  procurement  of  the  circuit  card  as  well  efforts  to 
prototype  the  card. 

5. 1.3. 9  Breadboard  ViPSIA 

Validation  of  the  design  concept  was  demonstrated  by  creating  a  wire  wrap  version  of  the 
ViPSIA  digital  section.  Discrete  components  were  utilized  with  wire  wrap  sockets  to 
create  a  logical  representation  of  the  digital  portion  of  the  circuit  card.  Once  the  design 
was  operating  as  required  the  refined  schematic  was  delivered  to  the  ViPSIA 
manufacturer  for  card  layout  and  construction. 

5.1.3.10  Develop,  Design  and  Manufacture  The  ViPSIA 

Estimates  from  three  sources  to  develop  the  layout,  design  the  circuit  card  and 
manufacture  initial  products  were  received.  A  local  company,  Galactic  Technologies, 
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Inc.  was  selected  based  upon  their  cost,  hardware  design  experience  and  desiun 
recommendations,  to  build  the  ViPSIA.  The  layout  for  the  board  was  completed,  desi^ 
for  the  FPGA  was  generated  based  upon  the  logical  design  contained  in  the  design  input, 
the  printed  circuit  board  was  procured  and  the  board  was  manufactured  and  tested. 

5.1.3.11  Test  The  ViPSIA 

Once  the  ViPSIA  was  manufactured  the  circuit  card  was  integrated  with  the  Texan  ET 
and  tested  to  ensure  that  all  functions  operated  correctly.  The  configuration  was  tested  in 
TBD  different  configurations  to  demonstrate  compatibility  of  the  VisionProbe  with 
various  manufacturer’s  PCs. 

5.1.3.12  Cable  Design  And  Manufacture 

Cables  for  the  VisionProbe  products  were  designed  and  manufactured.  A  Cable  is 
required  between  the  ViPSIA  PC  connector  and  the  data  acquisition  system  and  an 
internal  cable  is  required  between  the  ViPSIA  circuit  card  and  the  PC  connector  plate. 

5.1.3.13  Data  Acquisition  Interface  Development 

Discussions  between  DataWave  Technologies  and  CMI  engineers  revealed  that  the 
Experimenters  WorkBench  software  included  some  VisionProbe  routines  which 
supported  direct  communications  between  the  two  products.  Examination  of  this 
software  determined  that  it  communicated  with  the  NDI  1288  board  directly  through  it’s 
TMS3401 0  host  interface.  This  type  of  host  interface  is  not  supported  in  the  TMS34020. 
Recall  that  one  of  the  primary  design  requirements  for  the  VisionProbe  development  of 
this  product  was  to  not  require  DataWave  to  have  to  make  changes  to  their  product. 
Initial  descriptions  of  the  VisionProbe  stated  that  the  Experimenter’s 
WorkBenchWisionProbe  communication  was  performed  via  an  ASCII  or  binary  file.  We 
began  thinking  of  ways  to  simulate  the  old  host  communications,  but  this  approach  would 
cause  problems.  The  problem  is  that  the  host  communication  with  the  34010  is  tied  in 
with  the  hardware,  (the  GSP  could  detect  a  read  or  a  write  to  it’s  HSTDATA  register)  this 
means  that  the  PC  could  write  or  read  as  fast  as  it  wanted  to.  A  software  solution  would 
introduce  errors  because  we  would  not  be  able  to  tell  when  the  PC  has  read  or  written  to 
the  GSP  memory,  thus  resulting  in  data  loss.  The  solution  to  this  problem  is  to  use  the 
new  host  interface  and  make  the  appropriate  changes  to  the  DataWave  code.  These 
changes  will  not  be  major  because  they  all  occur  in  functions  that  are  easily  replaceable. 
The  new  functions  were  developed  by  CMI  so  that  it  will  recognize  TMS34010  or 
TMS34020  based  VisionProbe  and  set  the  communications  protocol  accordingly. 
DataWave  must  only  recompile  their  Experimenter’s  WorkBench  product  to  contain  the 
new  source  code.  The  interface  has  been  tested  at  DataWave’s  facihty  and  is  operating 
correctly.  The  VisionProbe  configuration  that  utilizes  the  DataWave  data  acquisition 
system  is  showm  in  below. 
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Fi^rc  3:  Swept  StimuliLS  Operation  Example 
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5.1.4  Commercialization 

The  sale  of  the  VisionProbe  to  vision  researchers  will  be  conducted  with  two  primary 
emphasis.  The  first  emphasis  will  be  to  market  the  VisionProbe  product  with 
DataWave’s  Experimenter’s  WorkBench.  The  VisionProbe  is  fully  integrated  with  the 
Experimenter’s  WorkBench  and  already  has  developed  a  list  of  potential  customers  based 
upon  their  use  of  the  previous  version  of  the  VisionProbe.  Plans  are  currently  being  made 
to  show  the  VisionProbe  with  the  Experimenter’s  WorkBench  at  the  Neuroscience 
convention  this  fall.  The  VisionProbe  is  also  available  as  a  standalone  vision  stimulator 
to  be  utilized  in  various  research  application  requiring  a  variety  of  programmable, 
controlled  visual  stimuli. 


5.2  Educational  Tool 

5.2.1  Description 

The  Educational  Tool  product  will  provide  the  clinician  or  technician  a  method  to  display 
the  effects  of  certain  diseases  on  the  human  vision  system  to  patients  and  their  families. 
The  purpose  of  the  educational  tool  is  to  provide  the  ability  to  view  a  scene  through  the 
eyes  of  an  individual  with  visual  degradation  due  to  disease  or  other  factors.  The  tool 
will  be  utilized  as  a  teaching  tool  to  show  people  how  the  vision  system  of  someone  with 
eye  disease  expresses  itself  It  is  based  on  the  AVS/Concentric  Ring  model  described 
earlier  in  the  report. 

The  educational  tool  shall  have  the  ability  to  display  pictures  of  representing  general 
cases  of  diseases  specified.  Their  shall  also  be  options  to  enter  specific  disease 
parameters,  (i.e.  damaged  areas  of  the  retina).  The  basis  of  the  educational  tool  shall  be 
the  Vision  Model  developed  for  the  A  VS  environment  for  the  SBIR  Phase  II  project. 

The  diseases  listed  below  will  be  modeled  within  the  educational  tool.  Specific  details  of 
each  of  the  models  were  developed  based  upon  discussions  with  subject  matter  experts. 

•  Glaucoma 

•  Photo  Toxicity 

5.2.2  System  Configuration 
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5.2.2  System  Configuration 

The  Educational  tool  is  an  PC  application  software  package  under  Windows  3.1  or  higher 
environment.  No  special  hardware  or  system  configuration  is  required  to  install  or 
operate  the  tool. 

5.2.3  Educational  Tool  Development  Tasks 

5.2.3. 1  Porting  the  AVS  Model  to  the  PC 

The  Vision  model  developed  for  laser  damage  in  the  SUN/AVS  environment  was  re¬ 
hosted  to  the  PC  environment  to  allow  a  broad  group  of  users  access  to  the  program.  All 
of  the  display  and  data  entry  routines  were  either  developed  or  converted  from  the 
SUN/AVS  software  developed  for  this  project. 

5.2.3. 2  Define  Disease  Parameters 

In  order  to  adequately  model  the  various  diseases  of  the  eye  in  the  educational  tool, 
interviews  and  research  was  conducted.  Vision  experts  at  the  University  of  Texas  Health 
Science  Center  were  contacted  and  provided  tutorials  to  the  design  engineers  regarding 
the  diseases  of  interest,  how  they  affect  an  individuals  vision  and  how  the  disease 
progresses.  Descriptions  of  each  of  the  diseases  to  be  modeled  were  provided  to  the 
software  engineers  for  use  in  the  implementation  of  the  models.  The  models  used  for  the 
diseases  largely  resembled  the  AVS/Concentric  Ring  model  presented  in  this  report. 

5.2.3.3  Design  and  Implement  the  Disease  Models 

Each  of  the  disease  models  were  translated  into  a  software  design.  The  design  was  coded 
in  C++  and  integrated  with  the  AVS  converted  software  (see  1. 2.3.1  above). 

5.2.4  Remaining  Development  Tasks 

The  following  tasks  remain  to  be  completed  in  order  to  make  the  Educational  Tool  a  fully 
marketable  commercial  product. 

•  Full  testing  must  be  conducted  on  the  software. 

•  Detailed  Disease  descriptions.  Help  features,  and  Users  Guides  must  be 
completed. 

•  Additional  diseases  may  be  included  into  the  tool. 

•  User  evaluation  and  recommendations  must  be  sought  in  order  to  improve  the 
product. 

6.  Concurrent  Research  Topics 

As  a  part  of  the  Phase  11  project,  basic  multi-disciplinary  research  at  UTS  A  was 
performed.  Dr.  Desai  of  the  electrical  engineering  department  developed  a  signal 
processing  model  for  the  effects  of  laser  damage  on  acuity.  Dr.  Tsin  of  the  biology 
department  conducted  experiments  on  the  effects  of  a  simulated  retinal  hemorrhage  on 
the  level  of  photopigment  bleaching  in  light  adapted  bovine  eyes.  Dr.  Dykes  of  the 
psychology  department  study  the  effect  of  simulated  retinal  hemorrhage  in  the  form  of 
electronic  noise  on  psychophysical  detection  experiment  performance.  Dr.  Gokhman  of 
the  mathematical  department  investigated  the  generation  of  Wavelets  as  possible  stimuli 
(other  than  Gabor  patches)  and  looked  into  the  modeling  of  the  effect  of  blood  as  a 
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scattering  media  in  the  vitreous  humor.  All  the  reports  have  been  included  in  Appendix 


7.  Summary 

The  present  Phase  D  brought  together  researchers  in  vision  and  signal  processing  from 
UTSA,  UTA,  UTHSC,  CMI  and  Brooks  Air  Force  Base.  This  research  resulted  in  four 
sets  of  experimental  results,  three  models  and  a  visualization  tool  which  may  be  used  at 
Brooks  Air  Force  Base  for  modeling  laser  eye  damage.  In  addition,  two  commercial  tools 
were  developed  and  a  vendor  has  been  established  for  one  tool  while  a  vendor  is  being 
sought  for  the  other. 
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The  Visualization  Tool  Description 


GOAL: 

This  model  was  conceived  as  a  visuali2ation  tool  to  illustrate  the  effect  on  a  scene  as 
viewed  by  a  person  with  normal  vision,  a  laser  lesion,  or  disease.  For  the  Air  Force  this 
model  has  been  implemented  with  A  VS  (so  that  it  will  interact  with  the  existing  ILPIM 
model  at  the  laser  lab)  and  for  the  commercial  package  it  has  been  implemented  on  a  PC 
with  windows. 

The  model  will  be  explained  in  a  serips  of  three  stages.  The  first  will  explain  how  we 
perturb  a  natural  scene  with  a  series  of  low  pass  filters  to  resemble  what  a  person  sees 
who  has  normal  vision  and  only  focuses  at  one  point  on  the  scene.  The  second  stage  will 
explain  how  we  further  degrade  the  image  to  reflect  the  damaged  retina  (laser  lesion  or 
disease).  The  third  stage  will  assume  that  the  user  has  selected  several  foviation  points  and 
explain  how  a  resultant  image  is  arrived  at  that  represents  the  total  information  viewed  by 
the  person. 

STAGE  1: 

A  digital  image  is  obtained.  The  user  selects  any  pixel  in  the  image  and  assumes  this  is  the 
center  of  the  region  in  the  field  of  view  that  is  focused  on  the  fovea  (herein  called  the 
foviation  point).  The  user  then  selects  up  to  15  eccentricities  that  may  extend  up  to  30 
degrees  from  the  fovea.  The  distance  between  the  retina  and  the  ring  is  called  eccentricity. 
A  source  of  contrast  sensitivity  functions  is  selected,  either  from  the  literature  or  the  CMI 
experiments.  A  genetic  algorithm  is  used  to  interporlate  the  contrast  sensitivity  data, 
generating  a  3-D  surface  with  contrast  sensitivity  as  the  amplitude  and  eccentricity  and 
frequency  as  the  independent  axis.  For  each  of  the  eccentricities  selected  by  the  user,  a 
contrast  sensitivity  function  is  arrived  at  fi-om  the  surface  that  was  just  fit  to  the 
experimental  data.  Next  a  low  pass  filter  is  developed  that  has  the  frequency  response  of 
the  corresponding  CSF  curve.  For  each  of  the  eccentricities  chosen,  a  copy  of  the  original 
image  is  filtered  with  the  low  pass  filter  that  corresponds  to  the  CSF  function  for  that 
eccentricity  (this  eliminates  any  boundary  effects  that  would  occur  if  one  to  filter  only  part 
of  the  image  with  the  corresponding  low  pass  filter).  Now  a  composite  image  is  formed. 
Each  pixel  in  the  original  image  lies  is  in  one  and  only  one  eccentricity  and  there  is  one  and 
only  one  filtered  image  that  corresponds  to  that  eccentricity.  Therefore  the  ij*  pixel  in  the 
original  image  can  be  replaced  by  the  ij*'’  pixel  in  the  filtered  image  that  corresponds  to  the 
correct  eccentricity.  This  results  in  an  image  that  is  increasingly  blurred  outward  fi-om  a 
user  selected  foviation  point  and  represents  the  loss  in  resolution  of  the  retina  due  to 
decreased  sampling  density  with  distance  fiom  the  fovea.  This  scene  is  as  would  be  seen 
by  a  person  with  normal  vision. 
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Figure  2.  Filtered  Scene  Due  to  Natural  Degredatton.  (X  indicates  fixation  point.) 


STAGE  2: 

When  the  image  of  a  retina  is  displayed  and  the  person  selects  the  eccentricities,  they  are 
also  ask  to  select  damaged  regions  on  the  retina  with  the  mouse.  These  may  be  user 
selected  or  the  retinal  image  may  be  of  a  damaged  retina  with  the  user  only  outlining  the 
actual  damaged  region.  For  laser  lesions  and  certain  diseases  (macular  degeneration, 
phototoxicity,  glaucoma,  etc.)  the  damaged  region  is  assumed  to  have  total  photoreceptor 
drop-out  and  therefore  the  pixels  in  the  damaged  region  are  replaced  with  the  average 
illumination  in  the  total  image.  For  Glaucoma  we  use  the  Humphry  field  map  to  arrive  at 
decreased  amplitudes  and  cutoff  points  for  the  low  pass  filters  corresponding  to  the 
effected  regions.  This  technique  is  easily  extended  to  any  disease  for  which  one  can 
determine  the  change  in  the  CSF.  The  result  is  a  scene  as  viewed  by  a  person  with 
damaged  vision  who  has  only  focused  his/her  eyes  at  one  point. 
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Figure  3.  Retinal  Image  with  Marked  Damage  Region 
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Figure  4.  Output  Scene  due  to  Retoial  Damage 


STAGE  3: 

This  next  stage  is  meant  to  represent  the  total  amount  of  information  available  to  a  person 
with  normal  or  damaged  vision  if  they  are  allowed  more  time  to  view  the  same  scene  and 
therefore  have  several  foviation  points.  We  assume  the  information  gained  at  each 
neighborhood  would  be  the  maximum  amount  from  the  different  foviations,  i.e.  the 
neighborhood  that  has  been  blurred  the  least.  Therefore  in  this  stage,  stages  1  and  2  are 
repeated  for  each  foviation  point  and  the  resultant  image  is  formed  by  selecting  pixel  ij 
from  the  scene  in  which  the  neighborhood  of  ij  was  filtered/blurred  the  least.  This  results 
in  a  scene  which  assumes  that  the  person  retained  the  maximum  possible  information  at 
each  foviation  point. 

Inputs:  Natural  Scene 

Number  of  Eccentricities 
Regions  of  disease  or  lesion 
One  or  more  foviation  points 

Output:  Scene  as  viewed  by  the  person  normal  or  afflicted 
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Appendix  II 


Detailed  Description  of  the  Inhomogenous  Retinal 

Cortical  Model 


Appendix  II 
Table  of  Contents 


I.  Introduction .  27 

1 . 1  Purpose  of  Study .  27 

1 .2  Measures  of  Visual  Performance .  27 

II.  Literature  Review .  28 

(Section  II  of  Appendix  II  is  available  for  review  at  AFRL/HEDO,  BAFB,  TX) 

2.1  Introduction 

2.2  Methods  of  Measuring  Visual  Performance 

2.3  Relevant  Details  of  Psychophysics 

2.4  VEP 

2.5  Basic  Vision  Facts 

2.6  Artificial  Scotoma  Literature 

2.7  Models  of  the  Vision  System 

III.  Research  Performed .  91 

3.1  VEP  Model . 

3.1.1  Introduction .  91 

3.1.2  Methods .  -  91 

3.1.3  Results .  95 

3.1.4  Conclusions .  97 

3.2  Best  Remaining  Model .  97 

3.3  Interpolation  of  Contrast  Sensitivity .  93 

IV.  Models,  Experiments,  and  Comparisons .  100 

4.1  IRC  Model . .  .  .  .  .  .  .  .  .  .  .  .  .  .  .  . 

4.1.1  Relationship  with  Ideal  Observer .  Iqq 

4.1.2  Summation  ofd’ .  lOi 

4.1.3  Computer  Simulation  of  IRC  Model  for  Gabor  Patch  CSF  Curves . 101 

4.2  Psychophysical  Experiments .  IO2 

4.2.1  Design  Details  and  Experimental  Protocol . 102 

4.2.2  Method .  IO2 

4.2.3  Airplane  Studies . 106 

4.2.4  Results . ' .  IO9 

4.3  Alternate  Model .  1 1 1 

V.  Conclusions . 

VI.  References .  216 


Simulating  the  Effects  of  Laser  Damage  to  the  Retina 

Final  Report 

List  of  Figures 

1:  Sample  psychometric  function  for  2AFC  experiment .  39 

FIGURE  2:  Sample  data  for  method  of  constant  stimuli  and  fitted  Psychometric  Function . 33 


P^^URE  3;  Decision  Space  for  Normal  Distributions  of  Equal  Variance . 37 

FIGURE  4:  Geometric  Explanation  for  equation  [5] . 3g 

FIGURE  5:  Sample  Experimental  Points  with  Fitted  Weibull  Function . 40 

FIGURE  6:  Reproduction  roulette  wheel . 42 

FIGURE  7:  VEP  Averager  as  digital  filter . 44 

FIGURE  8:  Magnitude  response  spectra  for  sample  VEP  averagers . 46 

FIGURE  9:  Steady  state  VEP  before  and  after  OS  filtering . 47 

FIGURE  10:  FFT  of  filtered  output  of  signal  shown  in  Figure  9 . . . 48 

FIGURE  II:  Photo  receptor  and  ganglion  cell  densities . 51 

FIGURE  12:  Fourier  Transform  of  Rectangular  Window  Function . 53 

FIGURE  13:  Gabor  Envelope,  Real  and  Imaginary  Terms  of  a  Gabor  Function . 55 

FIGURE  14:  Hilz  &  Cavonius  CSF  Family . .....57 

FIGURE  15:  Banks  CSF  Curve  Family . 58 

FIGURE  16:  Banks  Curve  Family  Linear  -Linear  Plot . 59 

FIGURE  17:  Banks  Family  of  Curves.  Linear  Log  Plot . 60 

FIGURE  18:  Stimulus  Response  Curves  for  Different  Levels  of  Background  adaptation . 62 

FIGURE  19:  Illustration  of  Vernier  Acuity . 65 

FIGURE  20:  Illustration  of  output  of  center  surround  stimulated  by  a  passing  edge . 67 

FIGURE  2 1 :  Summation  Charts .  m 


25 


FIGURE  22:  CS  for  a  stabilized  and  unstabilized  grating . 

figure  23:  Functional  block  diagram  of  the  vision  system  with  ideal  observer 

figure  24:  Foveal  center  in  Monkey,  Macaque  fascicularis . 

FIGURE  25:  Spatial  response  of  Difference  of  Gaussians . 

FIGURE  26:  Fourier  Transform  of  Difference  of  Gaussians . 

FIGURE  27:  Sample  cell  responses  for  simple  cortical  cells . 

FIGURE  28:  a.Texture  with  Gaussian  Biobs.  b.  Texture  with  Gabor  Blobs . 

FIGURE  29;  Sutter  Graham  Nonlinear  Model . 

FIGURE  30:  CSF  Curve  Built  up  From  SF  Channels . 

FIGURE  3 1 :  Ideal  Observer  Calculation  of  D’ . 

FIGURE  32:  Data  flow  through  VEP  model . 

FIGURE  33:  Steps  in  VEP  Model . 

FIGURE  34:  Results  from  the  VEP  model . 

FIGURE  35:  Family  of  Curves  from  Rovamo.  2D  and  3D  plots . 

FIGURE  36;  Genetic  Algorithm  fit  of  Banks  Data . 

FIGURE  37:  Timing  of  the  two-interval,  forced-choice  experiment . 

FIGURE  38:  Airplane  Stimuli . 

FIGURE  39:  Third  Harmonic  curv^es . 

FIGURE  40:  CSF  Curves . 

FIGURE  41:  Airplane  curves . 

FIGURE  42:  Spatial  Frequency  Channel  Model . 

FIGURE  43.  Over\dew  of  Calculations-Spatial  Frequency  Channel  Model . 

FIGURE  44:  Post-Channel  Calculations . 


. 72 

. 74 

. 77 

. 79 

. 80 

. 81 

. 82 

. 83 

. 84 

. 86 

. 93 

. 94 

. 96 

. 99 

. 99 

. 103 

. 106 

. 108 

. 109 

. no 

. Ill 

. 112 

. 113 


26 


1.1  Purpose  of  Study 


The  purpose  of  this  project  is  to  calculate  loss  of  visual  acuity  due  to  retinal  receptor  dropout.  The 
idea  of  modeling  visual  acuity  has  grown  from  military'  and  clinical  necessity.  Since  the  dawn  of 
the  Atomic  Age,  the  Air  Force  has  been  interested  in  eye  damage  caused  by  intense  light  which 
originally  came  from  atomic  bomb  flashes.  Later,  however,  a  much  more  common  threat 
appeared,  laser  light.  Lasers  have  become  widely  used  in  the  military  both  as  aiming  devices  and 
weapons.  There  is  constant  eye  danger  due  to  deliberate  enemy  action  as  well  as  accidents.  Even 
relatively  weak  laser  light  is  concentrated  on  a  small  part  of  the  retina  by  the  optics  of  the  eye  and 
can  cause  bums. 

Laser  damage  can  take  many  forms.  For  example,  hemorrhages  can  result  if  arteries  are  hit.  A 
strong,  direct  hit  to  the  optic  disk  could  presumably  destroy  the  entire  optic  nerve  blinding  the 
eye.  Since  the  topic  of  laser  eye  damage  is  so  vast,  the  scope  of  this  study  will  be  limited.  This 
study  assumes  retinal  loss  due  to  heat  or  occlusion  by  blood  but  without  underlying  nerves  being 
affected.  Furthermore,  retinal  loss  has  hitherto  been  treated  as  either  “all  or  nothing,”  though  some 
of  the  models  being  developed  would  allow  partial  losses  to  be  incorporated  without  modification. 

The  possibility  of  modeling  partial  damage  leads  to  another  use  for  the  models  being  developed: 
the  calculation  of  acuity  loss  due  to  various  eye  diseases  which  damage  the  retina  e.g.  glaucoma, 
retinitis  pigmentosa,  etc.).  Presently  no  disease  modeling  has  been  done;  however  the  possibility 
is  open  for  future  work. 

1.2  Measures  of  Visual  Performance 

Modeling  the  effects  of  retinal  damage  requires  metrics  for  the  quality  of  vision.  We  have  already 
chosen  visual  acuity  as  our  means  of  measurement,  but  unfortunately  there  is  no  single  definition 
for  this  term.  Two  important  types  of  acuity  are  Snellen,  or  the  ability  to  read  an  eye  chart,  and 
grating,  the  ability  to  detect  sinusoidal  intensity  variations  on  a  background  of  otherwise  constant 
intensity.  These  and  other  kinds  of  visual  acuity  are  described  below  in  section  [2.5.4].  The 
important  point  here  is  to  distinguish  the  concept  of  acuity  as  a  metric  for  visual  performance  with 
the  actual  measurement  of  acuity.  This  measurement  sounds  like  a  simple  process,  as  easy  as  hav¬ 
ing  a  patient  read  an  eye  chart.  Unfortunately,  this  is  not  nearly  accurate  enough  to  be  of  use  in 
this  study.  Two  general  methodologies  are  available  to  measure  acuity  (among  many  other 


things);  the  visual  evoked  potentials  (VEP)  and  psychophysics.  Both  arc  complex  and  have 
strengths  and  weaknesses. 

VEP  are  a  class  of  weak  electrical  signals  that  are  emitted  by  the  brain  in  response  to  visual 
stimuli.  It  has  been  known  for  many  years  that  the  brain  emits  tiny  electrical  signals 
summed  from  a  great  numbers  of  neurons;  some  of  the  signals  are  related  to  vision.  The 
VEP  are  taken  via  external  electrodes  placed  on  the  back  of  a  subject’s  skull  near  the  visual 
cortex  while  the  subject  is  watching  a  timed  image  of  varying  stimuli.  These  signals, 
though  a  very  noisy  sum,  can  be  related  to  the  image  the  subject  sees,  and  used  to  give  a 
measure  of  visual  acuity  by  testing  a  series  of  patterns,  and  detennining  which  give  a 
response  and  which  do  not. 

The  other  methodology  to  measure  visual  performance  is  psychophysics.  Psychophysics  is 
a  very  bioad  term  meaning  the  study  of  quantitative  relationships  between  stimulus  and 
perception.  It  can  pertain  to  measurement  of  any  of  the  senses.  For  example,  seminal 
psychophysical  experiments  were  used  to  study  the  ability  of  an  observer  to  distinguish 
between  two  different  weights.  Generally,  psychophysics  requires  the  use  of  cooperative 
subjects  able  to  understand  language,  implying  humans  at  least  a  few  years  old.  However, 
animals  with  special  training  have  also  been  used. 

II.  Literature  Review 

This  section  (sixty-six  pages)  has  been  withdrawn  from  this  report.  It  is  available  for 
review  at  the  AFRL/HEDO  offices  at  Brooks  Air  Force  Base,  Texas. 


III.  Research  performed 


3.1  VEP  Model 

3.1.1  Introduction 

In  a  previous  study,  a  VEP  based  model  of  retinal  damage  measurements  was  developed.  The 
model  used  data  taken  in  previous  research  performed  by  the  Air  Force.  In  a  series  of  experi¬ 
ments,  rhesus  monkeys  had  their  eyes  damaged  by  controlled  laser  blasts.  The  monkeys’  retinas 
were  studied  and  photographed  before  and  at  various  times  after  the  exposures,  in  increasing  peri¬ 
ods  from  days  to  months.  A  secondary  set  of  experiments  endeavored  to  measure  the  monkey’s 
visual  acuity  before  and  after  the  laser  exposures.  As  explained  in  section  [2.2],  the  only  way  to 
measure  an  animal’s  vision  without  special  training  is  through  the  taking  of  VEP  signals.  VEP 
studies  were  made  at  approximately  the  same  time  as  the  photographs,  using  the  counterphasing 
grating  stimuli  described  in  equation  [11].  Several  studies  were  made  on  each  monkey  before  the 
laser  exposure,  to  establish  a  baseline  for  normal  vision.  Each  study  consisted  of  a  series  of  VEP 
readings  taken  at  differing  spatial  frequencies.  Acuity  was  estimated  by  measuring  the  highest 
spatial  frequency  which  gave  measurable  VEP  signals. 

A  simple  model  was  developed  to  simulate  the  acuity  loss  and  thus  perhaps  alleviate  the  need  for 
monkey  experiments.  The  original  goal  of  the  VEP  model  was  to  calculate  loss  in  visual  acuity 
from  the  retinal  damage  as  judged  from  the  photographs.  Since  the  only  available  way  to  estimate 
acuity  was  based  on  VEP  amplitudes,  the  goal  was  later  simplified.  It  became  the  calculation  of 
VEP  amplitudes  in  a  damaged  retina  and  comparing  them  with  actual  recorded  values. 

3.1.2  Methods 

The  easiest  way  to  explain  the  VEP  model  is  to  first  show  the  methodology  for  a  normal  eye,  then 
see  the  effect  of  retinal  damage.  The  theoretical  basis  for  the  model  is  the  fact  that  VEP  are  a 
noisy,  approximate,  sum  of  signals  in  the  visual  cortex.  Mathematically,  the  model,  shown  in 
block  diagram  form  in  figure  [32],  is  similar  to  a  simple  neural  network.  Functionally,  it  divides 
the  retina  into  a  series  of  imaginary  concentric  rings.  Each  ring  is  given  a  weight,  the  determina¬ 
tion  of  which  will  be  explained  below.  The  terms  called  “ring  fraction  good”  are  all  set  to  one.  A 
pattern  such  as  a  sine  w'ave  grating  being  passed  though  the  model  would  have  its  contrast  sensi¬ 
tivity  calculated,  (generally  by  interpolation  as  explained  in  section  [3.3]),  for  each  ring.  The 
responses  for  all  the  rings  are  then  summed.  VEP  is  known  to  have  considerable  variation 
between  different  subjects,  so  the  summed  ring  responses  are  multiplied  by  a  subject  weight.  The 
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system  is  trained,  independently  for  each  subject,  with  a  simple  minimization  of  least  squares 
algorithm.  Input  spatial  frequency,  measured  VEP  amplimde,  and  a  subject  ID  are  entered.  The 
error  term,  the  difference  between  the  output  and  measured  VEP  amplitudes  is  squared  and 
summed,  separately  for  each  subject.  Jogging  the  subject  weight  and  recalculating  gives  a  slightly 
changed  error,  allowing  the  calculation  of  a  derivative  of  error  with  respect  to  weight.  This  deriva¬ 
tive  is  used  to  recalculate  the  subject  weight  in  a  direction  that  reduces  the  error.  Baseline  data  are 
input  to  the  model  and  the  model’s  outputs  compared  with  measured  ones.  The  subject  weight  is 
adjusted  to  give  the  least  square  error. 


^  i  Z  ^^’nSri^SF.  [66] 

I  =  0 

Equation  [66]  describes  the  calculations.  S  is  the  strength  of  the  input  signal,  Wj  is  the  subject 
weight,  Wj.i  is  the  weight  for  the  ith  ring,  g^j  is  the  good  fraction  for  the  ith  ring  (set  to  1  for  a 
healthy  eye),  and  CSFj  is  the  contrast  sensitivity  function  for  the  ith  ring. 

The  model  has  certain  underlying  assumptions  and  simplifications.  Firstly,  contrast  sensitivity 
(CS)  is  used  as  an  estimator  for  VEP  response.  Theoretical  justification  for  this  is  tenuous.  One 
may  argue  that  high  CS  is  an  indication  of  high  sensitivity  of  the  visual  system  and  likely  implies 
a  high  VEP  response  and  visa  versa.  CS  was  used  because  of  the  large  variety  of  data  in  the  vision 
literature.  Secondly,  CS  is  assumed  to  be  equal  at  any  point  on  any  concentric  ring.  This  is  of 
course  not  true,  CS  is  knovvm  to  vary  at  different  positions  around  the  eye.  The  VEP  data  supplied 
by  the  Air  Force  were  so  noisy  that  the  concentric  assumption  was  perfectly  good,  and  as  will  be 
seen  below,  the  use  of  CS  curves  worked  as  well  as  could  be  hoped. 

Figure  [33]  shows  an  example  of  how  retinal  damage  is  applied  to  the  model.  The  upper  left  of 
[33]  shows  a  photograph  of  a  damaged  retina.  The  upper  right  is  the  same  picture  with  damaged 
and  blood  occluded  areas  outlined  by  hand.  The  optic  disk  is  also  blanked  out.  Its  effect  on  the 
VEP  summation  is  negligible  because  of  its  eccentricity.  Its  real  use  in  this  implementation  is  to 
calibrate  the  photographs,  which  have  no  other  means  of  relating  pixels  to  degrees.  A  computer 
program  next  separates  out  the  damaged  areas,  as  shown  in  the  lower  left  picture.  The  damage  is 
scaled  and  placed  on  the  concentric  rings  shown  in  the  lower  right.  The  inner  most  rings  are 
totally  occluded,  in  this  case  and  the  ones  further  out,  only  partially  so. 
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The  cortical  magnification  factor  or  (CMF)  (Wilson,  Levi,  Maffei,  Rovamo,  De  Valois,  1990) 


FIGURE  33.  Steps  in  VEP  Model. 


describes  the  mapping  between  the  retina  and  visual  cortex.  This  is  the  ratio  of  the  length  of  a  tiny 
portion  of  the  cortex  divided  by  that  of  the  length  of  the  retina  that  maps  to  it.  In  some  literature, 
CMF  is  defined  as  a  ratio  of  areas  instead  of  lengths.  The  CMF  varies  with  eccentricity,  or  dis¬ 
tance  from  the  fovea,  or  center  of  the  retina,  in  degrees. 


0 


MciQ)  =M^/(1+^) 


[67] 


94 


VEP  Model 


Spatial  frequency 


V 


Low 

Eccentricity 

(Fovea) 

ium 

Eccentricity  V 


amp 


arrp 


rO 


High 

^  Eccentricity 


Amplitude  Response 
vs 

Spatial  Frequency 
for 

Different  Eccentricities 
Ring  Weights 


Amplitude  Response 


Ring  Fraction  good 


FIGURE  32.  Data  flow  through  \TP  model 


After  the  damaged  areas  have  been  superimposed  on  the  rings,  equation  [66]  is  again  used  to  give 
the  VEP  response  predicted  by  the  model.  The  only  difference  from  the  training  case  with  undam¬ 
aged  eyes,  is  that  the  gj-j  term  for  each  ring  is  set  equal  to  the  fraction  of  the  ring  remaining. 

Finally,  the  calculation  of  the  ring  weights  must  be  explained.  The  original  idea  was  to  calculate 
them  and  the  subject  weights  with  a  single  algorithm,  similar  to  a  backpropagation  neural  network 
(Widrow  &  Steams,  1985).  This  method  proved  impractical  because  the  training  data  were  too 
noisy,  as  is  seen  in  the  top  chart  of  figure  [34],  Another  way  to  calculate  the  ring  weights  was 
developed,  using  cortical  magnification  factors. 
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Equation  [67]  is  a  formula  commonly  found  in  the  literature  to  describe  the  change  in  CMF  with 
eccentricity.  Mq(©)  is  the  cortical  magnification  factor  as  a  function  of  eccentricity.  0  is  the 
eccentricity,  Mp  is  the  magnification  found  in  the  fovea.  The  parameter  called  E2  is  the  eccentric¬ 
ity  where  the  cortical  magnification  drops  to  half  the  value  at  the  fovea.  Me  is  a  dimensionless 
ratio  so  E2  must  have  the  same  units  as  0,  generally  degrees.  Unfortunately,  there  is  no  agreement 
on  a  value  for  E2,  or  even  on  a  narrow  range  of  values.  It  varies  drastically  with  the  nature  of  what 
is  being  measured.  According  to  Wilson  et  al.  (1990)  there  are  two  ranges  of  values  for  Eo, 
depending  on  what  is  being  measured.  One  range  was  from  .3  to  .9  degrees,  the  other  1.5  to  4. 

Once  a  value  for  E2  has  been  selected,  the  ring  weights  are  calculated  by  multiplying  the  area  of 
each  ring  by  the  square  of  the  cortical  magnification  factor.  This  squaring  converts  a  linear  mea¬ 
sure  to  an  area  measure.  The  effect  of  the  E2  parameter  is  to  set  the  relative  weight  of  the  fovea 
versus  the  rest  of  the  retina.  The  lower  the  value  for  E2,  the  faster  the  dropoff  with  eccentricity  and 
the  higher  the  relative  weight  of  the  fovea.  Conferring  with  Dr.  Randolf  Glickman,  (1993,  per¬ 
sonal  communication)  one  of  the  scientists  who  had  taken  the  original  VEP  measurements,  it  was 
learned  that  the  fovea  is  responsible  for  “almost  half  the  VEP  signal  strength.”  To  translate  this 
statement  into  a  value  for  E2,  Mathematica  was  used  to  calculate  the  first  thirty  ring  values  and  the 
fi'action  that  fell  inside  the  fovea,  assumed  here  to  be  a  circle  of  radius  2.5  degrees. 

3.1.3  Results 

Figure  [34]  shows  some  results  of  the  model.  The  top  graph  shows  the  fit  to  an  undamaged  eye. 
Note  that  the  peak  of  the  solid  curve  matches  the  peak  value  of  the  experimental  points.  Because 
of  the  extreme  noisiness  of  the  data  points,  as  evident  in  that  graph,  a  different  training  algorithm 
was  developed.  Instead  of  minimizing  the  squared  errors,  the  subject  weights  were  set  so  that  the 
peak  of  the  curve  for  undamaged  data  matched  the  peak  of  the  raw  data  points.  The  center  graph 
shows  a  predicted  curve  versus,  experimental  points  taken  just  after  exposure.  The  third  graph 
shows  predicted  and  actual  data  for  two  weeks  after  the  exposure,  when  blood  had  time  to  be 
reabsorbed  into  the  body,  slightly  improving  sight. 
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3.1.4  Conclusions 


The  results  show  promise,  as  can  be  seen  in  figure  [34],  but  the  work  was  hampered  by  the 
extreme  noisiness  of  the  original  data  supplied  by  the  Air  Force.  Work  on  this  model  was  dropped 
in  favor  of  the  psychophysical  modeling  described  in  the  rest  of  this  paper.  If  new  \TEP  studies 
were  to  be  undertaken,  care  should  be  taken  to  make  the  measurements  as  repeatable  as  possible. 
There  is  anecdotal  evidence  that  the  researchers  in  the  original  studies  were  not  interested  in 
amplitudes,  only  in  the  presence  or  absence  of  a  signal.  The  ring  weights  would  also  be  worthy  of 
further  study.  It  should  be  possible  to  vary  the  E2  parameter  to  improve  the  fit  and  to  include  some 
data  from  damaged  eyes  in  the  training.  One  could  also  return  to  the  original  scheme  of  calculat¬ 
ing  the  ring  weights  by  training. 

3.2  Best  Remaining  Model 

Another  way  to  look  at  the  vision  damage  problem  could  be  called  the  “Best  Remaining  Model.” 
This  model  is  based  on  the  fact  that  most  of  the  visual  acuity  is  based  on  the  center  of  the  fovea 
where  the  cones  are  most  dense.  As  the  eyes  move  around,  the  fovea  is  fixated  on  points  of  inter¬ 
est.  If  part  of  the  retina  is  destroyed  by  a  lesion,  the  loss  in  acuity  will  depend  on  the  position  of 
the  lesion.  According  to  this  model,  if  it  misses  the  fovea,  no  loss  in  acuity  will  result  If  the  lesion 
destroys  the  fovea,  the  loss  will  depend  on  the  observer’s  reaction.  If  he/she  looks  directly  at  the 
scene  of  interest,  the  best  remaining  portions  of  the  retina  will  not  be  focused  on  the  target  but  at 
some  distance  from  it.  Later,  the  observer  may  teach  himself/herself  to  look  away  from  w'hat  he/ 
she  is  trying  to  see,  thus  placing  the  best  remaining  part  of  the  retina  on  the  target  of  interest. 

This  model,  though  totally  different  from  the  VHP  model,  also  suffers  from  severe  weaknesses.  It 
implies  that  a  large  lesion  will  do  no  damage  to  vision  if  it  misses  the  fovea,  even  by  a  small 
amount.  Even  if  the  lesion  traversed  the  center  of  the  fovea,  the  best  remaining  model  would 
imply  that  vision  would  be  unimpaired.  Intuitively,  this  seems  hard  to  believe.  Furthermore,  it 
does  not  take  into  account  eye  motion  or  any  other  temporal  effects.  If  the  observer  focuses  the 
best  remaining  part  of  his  retina  on  a  target,  he  will  not  be  able  to  keep  his  eye  motionless.  It 
seems  likely  that  the  eye  will  at  certain  times  move  into  a  position  such  that  the  target  is  focused 
on  the  damaged  part  of  the  retina.  Clearly,  the  best  remaining  model  is  no  solution  to  this  prob¬ 
lem,  it  is  at  best,  a  guideline  for  attacking  it. 
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3.3  Interpolation  of  Contrast  Sensitivity 

Both  the  VHP  and  psychophysical  vision  damage  models  require  interpolation  of  contrast  sensi¬ 
tivity.  CS  measurements  are  expensive  and  slow,  so  only  relatively  few  are  done  in  any  study. 
Both  models  require  CS  points  to  be  available  for  any  combination  of  eccentricity  and  SF. 

There  are  two  basic  ways  to  approach  the  CSF  interpolation  problem.  One  is  to  try  to  find  an  ana¬ 
lytic  function  that  describes  the  shape  of  the  cur\'es,  possibly  with  theoretical  justification.  Such  a 
formula  would  have  parameters  to  be  fitted  to  the  data  by  minimizing  a  sum  of  squares,  giving  an 
equation  for  CSF  given  the  eccentricity  and  SF.  The  second  way  is  simple  interpolation  such  as  a 
spline  fit. 

One  such  function  attempted  is  shov.Ti  in  equation  [68]. 

CSF  =  [68] 

This  equation  fitted  CSF  to  SF  for  a  fixed  eccentricity.  Four  parameters,  a,b,c,d  were  calculated 
using  a  least  squares  routine.  This  worked  fairly  well  for  cuA'es  at  a  single  eccentricity.  One 
approach  to  extend  it  to  two  dimensions  was  to  calculate  the  four  parameters  for  each  different 
eccentricity,  then  interpolate  each  parameter,  for  eccentricities  not  given.  This  worked  for  one  set 
of  data,  the  Rovamo  showm  here  in  figure  [35],  but  failed  with  other  sets  of  data.  The  problem  was 
that  some  of  the  parameters  would  change  sign  during  interpolation  causing  instabilities.  This 
approach  was  finally  abandoned  for  a  double  spline  fit.  Each  SF  curv'e  for  a  fixed  eccentricity  was 
interpolated  by  spatial  frequency  using  cubic  splines,  giving  a  curve  with  many  evenly  spaced  SF 
points.  A  second  series  of  spline  fits  was  done  across  the  eccentricities  for  each  SF.  The  resulting 
data  were  put  into  a  table  and  read  into  a  program  to  give  fast  access.  Recently,  a  new  method  has 
been  tried.  The  frequency  response  of  a  Difference  of  Gaussian  filter  was  studied  and  it  was  noted 
that  it  had  a  considerable  resemblance  to  a  CSF  curve.  This  was  clearly  no  surprise,  because  the 
retina  is  believed  to  be  full  of  DoG  filters  converting  input  signals  from  the  retina  to  the  ganglion 
cells.  The  main  difference  was  that  an  unmodified  DoG  filter  is  band  pass  and  transmits  no  DC 
signal  what  so  ever.  A  CSF  cur\'e,  on  the  other  hand,  may  show  drop  off  at  low  frequencies  but  not 
cutoff.  Some  CSF  curves  show  virtually  no  dropoff  at  low  frequencies.  It  was  therefore  decided  to 
try  to  simulate  a  CSF  curve  with  a  DoG  by  shifting  the  latter  in  the  frequency  domain. 

Equation  [69]  shows  the  Fourier  transform  of  a  DoG  multiplied  by  A,  an  amplitude  parameter  to 
be  determined.  The  term  fg  represents  the  frequency  shift.  A  zero  value  for  fs  indicates  that  both 
terms  go  to  one  when  SF  goes  to  zero,  a  non-zero  value  will  cause  CSF  to  go  positive  or  negative 
(a  meaningless  case)  when  SF  is  zero. 
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[69] 


CSF  =  A(e^ 

The  parameters  fs,  si,  s2,  and  A  were  fit  with  a  genetic  algorithm,  described  in  section  [2.3.3. 1] 
using  minimizing  least  square  error  between  experimental  and  calculated  data  points.  Results 
from  fits  on  the  Banks  data  are  shown  in  figure  [36].  To  date,  attempts  to  extend  this  method  to 
three  dimensions  (CSF  versus  eccentricity  and  SF)  have  failed,  though  other  techniques  can  be 
used  to  interpolate  between  eccentricities. 


FIGURE  35.  FamUy  of  Curves  from  Rovamo.  2  and  3d  plots. 
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IV.  Models,  Experiments,  and  Comparisons. 


4.1  IRC  model 

4.1.1  Relationship  with  idea!  observer  theory 

As  explained  above,  ideal  observer  theory  does  not  provide  an  accurate  simulation  of  human 
vision,  nor  is  it  intended  to.  Yet  one  might  ponder  the  question:  “since  ideal  observer  theory  is 
based  on  the  best  known  information  about  the  vision  system  as  well  as  knowm  principles  of  phys¬ 
ics  and  statistics,  why  could  it  not  be  applied  to  vision  modeling?”  The  Inhomogeneous  Retino- 
Cortical  (IRC)  model  has  been  developed  by  professor  Wilson  Geisler  to  modify  ideal  observer 
theory  to  this  end.  The  concept  is  to  deviate  from  the  ideal  as  little  as  possible  while  allowing  cer¬ 
tain  parameters  to  float  as  necessary  to  match  actual  performance  for  a  particular  subject.  Once 
the  parameters  have  been  fit,  the  model  can  be  tested. 

The  IRC  model  in  its  current  state  simulates  the  optics  of  the  eye,  the  sampling  density  of  the  ret¬ 
ina  and  ganglion  cells,  parvo  cell  responses,  cortical  non-linearities  and  an  ideal  detector.  Each  of 
these  has  been  dealt  with  in  sections  [2.5]  and  [2.7].  The  optics  and  sampling  density  are  as 
described  by  Campbell  and  Gubish  (1966)  and  Curcio  et  al.  (1990).  The  ganglion  cell  responses 
are  described  with  spatial  DoGs  following  the  work  of  Enroth-Cugell  and  Croner  and  Kaplan 
(1995).  IRC  deviates  from  ideal  theory  here,  however.  As  mentioned  above,  there  is  a  time  delay 
between  the  center  and  surround  that  changes  the  shape  of  CSF  curves  according  to  the  temporal 
shape  of  the  stimulus.  Rather  than  complicate  the  model  by  simulating  these  temporal  effects 
directly,  the  relative  weights  between  center  and  surround  are  arrived  at  through  parameter  fitting. 
Continuous  or  slowly  changing  stimuli  show  low  frequency  dropoff  in  the  CSF  and  have  center 
and  surround  weights  close  to  .5.  Faster  changing  stimuli  show  little  low  frequency  dropoff  and 
have  center  weights  much  higher  than  surround. 

The  cortex  is  not  modeled  except  for  m'o  non-linearities.  The  first  is  rectification.  Cortical  cells 
unlike  ganglion  cells,  have  a  low  spontaneous  firing  rate  and  thus  cannot  transmit  a  negative  sig¬ 
nal.  The  other  is  the  contrast  exponent  described  above  in  section  [2.7.6].  The  exponent  is  one  of 
the  parameters  to  be  fit. 

A  d  is  calculated  with  the  same  method  as  figure  [31].  Ganglion  cell  responses  are  calculated  and 
the  cortical  factors  applied.  The  process  is  repeated  for  two  images  and  the  results  passed  to  a  d’ 
calculator.  This  d’  calculator  uses  a  noise  distribution  that  is  assumed  to  have  a  variance  propor- 
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tional  to  the  mean.  A  Poisson  distribution  is  characterized  by  variance  equal  to  the  mean.  Cur¬ 
rently,  the  variance  is  set  to  1.2  times  the  mean.  One  might  call  this  a  “Poisson  like”  distribution. 

4.1.2  Summation  of  d’ 

The  models  described  below  need  a  means  of  summing  many  values  for  d’.  d’  will  be  calculated 
by  comparing  cortical  cell  responses  for  two  images.  A  value  of  d’  is  calculated  for  each  individ¬ 
ual  cell.  These  d’  values  are  of  no  use  by  themselves  and  must  be  summed  to  provide  an  overall  d’ 
between  the  two  images.  To  sum  d’,  one  uses  its  definition  as  distance  divided  by  standard  devia¬ 
tion.  If  all  cell  outputs  are  assumed  to  be  independent,  distance  is  calculated  in  the  standard 
Euclidian  way,  the  square  root  of  the  sum  of  the  squares.  Likewise,  the  variance  of  a  random  vari¬ 
able  which  is  the  sum  of  other  random  variables  is  the  some  of  the  variances  of  those  variable. 
Equation  [70]  gives  the  final  result,  where  aj  is  the  output  of  the  ith  cell  for  one  image,  bj  is  the 
output  of  the  same  cell  for  the  other  image. 


d'  = 


+var(i>,.) 


[70] 


4.1.3  Computer  Simulation  of  IRC  model  for  Gabor  patch  CSF  curves 

A  computer  simulation  of  the  IRC  will  be  used  to  calculate  the  d’  values  for  various  patterns 
placed  on  a  simulated  retina.  This  program,  developed  in  FORTRAN  77  by  Professor  W.  Geisler, 
uses  the  IRC  model  to  calculate  ganglion  cell  outputs.  It  has  been  slightly  modified  for  the  pur¬ 
pose  of  the  present  work,  to  handle  large,  complex,  images  It  starts  by  assuming  Poisson  distrib¬ 
uted  noise  for  each  receptor  in  the  retina  as  described  in  the  Ideal  Observer  section,  [2.7.8].  A 
simulated  retina  is  used  for  the  spatial  distribution  of  the  inputs.  The  layout  is  based  on  the  data 
from  Curcio  (1990),  and  spatial  irregularities  in  the  hexagonal  sampling  matrix  are  also  modeled. 
Optical  blurring  in  the  eye  is  simulated  with  a  point  spread  function  and  ganglion  cell  outputs  are 
simulated  by  means  of  difference  of  Gaussian  formulas,  who’s  sizes  also  change  with  eccentricity. 
A  value  for  d’  for  a  stimulus  versus  a  blank  image  is  calculated,  according  to  formula  [70]. 

This  simulation  will  be  run  on  various  stimuli  being  tested  with  psychophysical  experiments  on 
human  observers,  especially  on  the  third  harmonic  cycle  summation  series  of  experiments. 
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4.2  Psychophysical  experiments 

4.2.1  Design  details  and  experimental  protocol 

Since  the  main  thrust  of  the  project  is  estimating  visual  acuity  under  performance  conditions,  the 
VEP  has  been  determined  to  be  a  poor  measure  for  the  reasons  stated.  A  new  psychophysical 
study  is  being  proposed  to  provide  performance  data.  This  study  generated  a  complete  set  of  CSF 
curves.  By  generating  a  new  set  of  CSF  curves,  our  stimulus  conditions  can  be  compared  with 
those  of  CSF  curves  in  the  literature  to  provide  some  idea  of  how  stimuli  may  effect  other  parts  of 
the  experiment. 

Several  sets  of  stimuli  were  used  to  develop  models  of  the  loss  of  vision  from  laser  damage.  This 
model  of  vision  loss  was  used  to  develop  visualization  tools  of  scenes  after  laser  damage  has 
occurred.  The  proposed  stimuli  will  be  used  to: 

•  measure  contrast  sensitivity  at  various  eccentricities  under  the  same  stimuli  conditions  as  a 
VEP  experiment,  so  that  a  comparison  of  the  two  methods  can  be  made 

•  measure  contrast  sensitivity  for  Gabor  patches  at  various  eccentricities  in  order  to  determine 
the  visual  capabilities  of  the  retina  at  various  eccentricities 

•  measure  the  ability  to  detect  a  natural  object  in  the  scene  at  various  eccentricities 

•  measure  spatial  summation  of  the  third  harmonic  of  gratings  with  and  without  the  presence 
of  a  fundamental  frequency 

4.2.2  Method 

All  experiments  will  be  conducted  in  a  quiet  room,  and  the  experiments  will  be  conducted  using  a 
VMI  monitor  driven  by  a  Micro  VAX  II  computer  with  Parallax  graphics  cards. 

An  additional  four  bits  of  resolution  have  been  added  to  the  brightness  values  of  the  stimuli, 
bringing  the  total  resolution  of  the  stimuli  to  rv.'elve  bits.  This  is  accomplished  by  building  a  cir¬ 
cuit  which  adds  the  signals  from  the  red  gun  to  that  of  the  green  gun,  which  is  the  main  output. 

Viewing  distances  have  been  determined  to  be: 

•  1 14.6  cm  for  spatial  frequencies  of  0.5,  1,  and  2  cycles  per  degree  (CPD) 

•  229.2  for  spatial  frequencies  of  4  and  8  CPD 

•  458.4  cm  for  spatial  frequencies  of  16  CPD 

•  573  cm  for  spatial  frequencies  of  32  CPD 

Subjects:  Three  volunteers,  one  female  aged  3 1 ,  and  tu'o  males,  aged  24  and  3 1  have  participated 
in  this  study.  Each  participant  was  given  an  ophthalmological  examination  to  ensure  normal 
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vision.  Each  participants  has  visual  acuity  corrected  to  at  least  20/20,  normal  visual  fields  and 
color  vision,  and  is  not  hypersensitive  to  light. 

This  section  describes  all  experiments  except  when  noted  otherwise.  Subject  has  a  chin  rest  to 
control  the  distance  to  the  screen.  The  computer  communicates  with  the  subject  via  a  voice  syn¬ 
thesizer.  At  the  start  of  each  experiment,  the  S3mthesizer  says,  “start.”  Next  there  are  two  stimulus 
intervals  of  200  milliseconds  (ms)  each  with  a  500  ms  pause  between  them.  Finally,  there  is  a 
variable  timed  response  interval.  This  experiment  is  referred  to  as  a  two-interval,  forced-choice 
experiment  (showm  in  Figure  [37]).  The  method  used  to  determine  the  contrast  threshold  involves 
a  three  down/one  up  staircase  with  reversals. [25] 

Two-Inter\'al,  Forced-Choice  Experiment  (4  seconds  per  trial  with  1  second  rest) 


time  (s) 


FIGURE  37.  Timing  of  the  two-interval,  forced-choice  experiment. 


For  both  the  first  and  second  sets  of  stimuli,  during  the  first  interval  a  pattern  of  stripes  (either  hor¬ 
izontal  or  vertical)  appears  on  the  screen.  Then,  during  the  second  interval  a  second  pattern, 
orthogonal  to  the  first,  appears.  The  third  set  of  stimuli  will  use  a  pseudo-natural  scene  (an  air¬ 
craft)  in  place  of  the  stripes  used  previously.  Just  as  before,  the  aircraft  will  appear  at  two  orthog¬ 
onal  orientations 

During  the  response  interval  for  the  first  three  sets  of  stimuli,  the  subject  presses  a  key  to  indicate 
which  pattern  appeared  in  the  first  interval.  The  voice  synthesizer  responds  with  “correct”  or 
“wrong”  depending  on  the  response  given.  The  response  interval  varies  fi-om  a  fraction  of  a  sec¬ 
ond  if  the  response  is  given  while  the  patterns  are  still  up,  to  about  1 .5  seconds  if  no  response  is 
made.  At  that  point  the  response  is  automatically  declared  wrong. 

The  third  set  of  stimuli  differs  conceptually  from  the  others.  The  fundamental  pattern  will  extend 
across  the  screen  and  will  not  change  during  the  experiment.  The  number  of  fundamental  cycles 
will  be  fixed,  but  the  number  of  harmonic  cycles  will  vary.  In  this  case,  all  patterns  will  be  hori- 
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zontal.  The  subject  will  be  only  trying  to  detect  the  presence  of  the  additional  harmonic.  During 
the  response  interval,  the  subject  presses  a  key  to  indicate  if  the  pattern  was  detected. 

There  are  three  sets  of  experimental  stimuli  which  are  detailed  below. 

The  first  set  of  stimuli  Table  1  consists  of  a  series  of  Gabor  patches  at  varying  eccentricities  used 
to  produce  CSFs  at  varying  eccentricities.  These  will  be  compared  to  the  literature  to  check  for  the 
effects  of  different  stimuli  conditions.  Gabor  patches  will  have  5  cycles  contained  in  2  standard 
deviations  of  the  Gaussian  envelope. 


Table  1  Psychophysical  Experiments  —  Gabor  Patches 


Stimulus  Shape 

Eccentricity 

Spacial  frequency 

Number  of  points 

Entire  Screen 

n/a 

1,2,4,8,16,32 

6 

Gabor  Patches  on 
Temporal  Meridian 

0,2,5,10,30 

cutoff/powers  of  2 

-22 

Gabor  patches  on 
Nasal  Meridian 

5,10,30 

cutoff/powers  of  2 

15 

The  second  set  of  stimuli  (Table  2)  will  consist  of  aircraft  images  with  and  without  missing  areas 
to  measure  the  effects  of  loss  of  parts  of  the  retina  on  the  ability  to  detect  objects  in  a  scene. 


Table  2  Psychophysical  Experiments  --  Pseudo-Natural  Scenes 


Stimulus  Shape 

Eccentricity 

Scotoma  size 

Number  of  points 

Aircraft;  1 .25°  long 

0 

0 

1 

Aircraft:  1.5°  long 

0 

0 

1 

Aircraft:  2.0°  long 

0 

0 

1 

Aircraft:  2.5°  long 

0 

0 

1 

Aircraft:  5.0°  long 

0 

0 

1 

Aircraft;  6.0°  long 

0 

0 

1 

Aircraft:  1.25°  long 

0 

1 

1 

Aircraft:  1.5°  long 

0 

1 

1 

Aircraft:  2.0°  long 

0 

1 

1 

Aircraft;  2.5°  long 

0 

1 

1 

Aircraft;  3.0°  long 

0 

1 

1 

Aircraft:  5.0°  long 

0 

1 

1 
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Table  2  Psychophysical  Experiments  -  Pseudo-Natural  Scenes 


Stimulus  Shape 

Eccentricity 

Scotoma  size 

Number  of  points 

Aircraft:  2.0°  long 

0 

2 

1 

Aircraft:  3.0°  long 

0 

2 

1 

Aircraft:  5.0°  long 

0 

2 

1 

Aircraft:  6.0°  long 

0 

2 

1 

The  third  set  of  stimuli  (Table  3)  will  consist  of  a  set  of  experiments  to  measure  grating  summa¬ 
tion  in  the  presence  of  a  fundamental.  This  will  allow  fine  tuning  of  the  current  theories  of  spatial 
frequency  channels,  and,  thus,  be  useful  in  perfecting  the  psychophysical  model  of  vision  loss. 
Cycle  summation  data  has  been  measured  for  Gabor  patches  [26],  but  to  our  knowledge,  it  has 
never  been  measured  for  a  harmonic  in  the  presence  of  its  fundamental. 

The  reasoning  behind  this  fourth  set  of  stimuli  is  based  on  the  decomposition  of  an  image  into 
spatial  frequency  channels.  Recall  that  an  edge  transforms  into  a  series  of  odd  harmonics  in  the 
frequency  domain.  Furthermore,  there  is  a  well  known  effect  of  cycle  summation.  A  grating  of 
two  cycles  is  more  visible  than  one  of  one  cycle  and  a  grating  of  5  cycles  is  more  visible  yet. 
Above  about  five  to  seven  cycles,  however,  the  grating  no  longer  becomes  more  visible.  What  is 
not  known  yet,  however,  is  if  this  cycle  summation  effect  remains  unchanged  for  an  odd  harmonic 
in  the  presence  of  its  fundamental. 

Most  natural  scenes  are  full  of  edges,  therefore,  any  model  for  detection  of  objects  that  depends 
on  spatial  frequency  channels  requires  information  about  how  odd  harmonics  sum  in  the  presence 
of  a  fundamental.  The  experiments  described  in  Table  4  should  provide  at  least  part  of  this  infor¬ 
mation. 


Table  3  Psychophysical  Experiments  —  Grating  Summation 


Stimulus  Shape 

Eccentricity 

Spatial 

Frequency 

Number  of 
Cycles 

Number  of 
Points 

Grating  + 

Third  Harmonic 

0 

1,2,4 

1,2,3,5,10 

45 

Third  Harmonic 

Only 

0 

3,6,12 

1,2,3,5,10 

45 
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FIGURE  38.  Airplane  Stimuli. 


Figure  [38]  shows  the  generation  of  the  airplane  stimuli.  The  upper  left  shows  the  original  picture 
of  the  plane.  The  background  has  been  made  black  with  manual  retouching.  In  the  upper  right,  the 
figure  has  been  centered  and  placed  in  a  background  of  intensity  equal  to  the  average  brightness 
of  the  plane.  In  the  two  lower  figures,  artificial  scotomas  of  different  sizes  are  shown. 


The  section  on  artificial  scotomas  [2.6.1]  describes  a  series  of  methods  to  stabilize  portions  of 
images  on  the  retina,  creating  the  effect  of  a  scotoma  or  dead  zone.  The  experiment  described  here 
achieves  a  similar  effect  without  any  special  equipment.  As  described  above,  it  uses  a  stimulus 
presentation  which  is  short  enough  (200  ms)  that  the  eye  will  not  move  during  that  interval,  thus 
no  stabilization  is  needed. 

The  airplane  experiments  consist  of  contrast  sensitivity  measurement  of  the  airplane  versus  a 
blank  background  of  the  same  intensity.  It  would  have  been  possible  to  do  90  degree  rotations  of 
the  airplane,  as  was  done  on  the  gratings.  This  was  not  done  for  two  reasons.  The  airplane  experi¬ 
ments  are  meant  to  be  simulations  of  natural  scenes,  whereas  the  Gabor  gratings  are  not.  Rotating 
the  airplane  90  degrees  is  artificial  and  not  compatible  with  the  concept  of  simulation  of  a  natural 
scene.  The  other  reason  is  that  rotating  a  long,  thin,  object  such  as  an  airplane  would  expose 
widely  different  parts  of  the  retina  to  the  image.  Nothing  on  the  vertical  meridian  has  been  tested 
in  these  experiments,  so  there  is  no  information  available  to  calibrate  a  model,  other  than  on  the 
horizontal  meridian. 
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LOG  CONTRAST  SENSITIVITY  CONTRAST  SENSITIVITY 


Cycle  Summation  at  0°  Eccentricity 


GABOR  STANDARD  DEVIATION  (cycles) 


GABOR  STANDARD  DEVIATION  (cycles) 


FIGURE  39.  Third  Harmonic  curves.  Dark  Lines  represent  IRC  simulations.  Triangles  represent 

presence  of  fundamental,  Circles,  its  absence. 
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FIGURE  41,  Airplane  curves.  Dark  Lines  represent  IRC  simulations.  Points  represent  experiment. 


4.3  Alternate  model 


This  alternate  model,  based  on  spatial  frequency  channels,  is  included  to  show  possible  alterna¬ 
tives  to  the  II^C,  It  was  never  implemented.  The  alternate  model  was  an  attempt  to  improve  on  the 
VEP  and  best  remaining  models,  but  to  be  simple  enough  to  be  easily  implemented  on  a  com¬ 
puter.  Like  the  other  models,  it  ignored  temporal  effects.  Figure  [42]  illustrates  an  overview  of  this 
spatial  frequency  chamiel  model.  In  a  similar  manner  to  the  IRC  model  (figure  [31]),  two  stimuli 
are  fed  to  a  routine  which  calculates  a  value  for  d’  representing  the  perceptual  difference  between 
the  two.  Figure  [42]  shows  a  block  diagram  of  the  processing  of  each  stimulus.  The  input  to  the 
retina  passed  through  non-linearities  related  to  control  of  sensitivity  (section  [2.5.2]  in  the  vision 
system.  Next  is  a  large  series  of  tuned  and  oriented  spatial  frequency  channels  (section  [2.7.7]) 
which  cover  the  entire  retina,  three  of  which  are  illustrated.  Each  channel  is  simulated  as  a  digital 
filter  which  is  convolved  with  a  portion  of  the  input  image.  These  outputs  are  processed  further  as 
shown  in  figure  [44]  before  going  to  the  d’  calculator.  Note  that  a  provision  is  made  for  interac¬ 
tions  between  other  channels,  though  early  implementations  of  the  model  will  probably  ignore 
these  interactions.  The  output  of  the  filter  convolution  is  passed  to  cortical  contrast  acceleration 
non-linearities  as  discussed  in  section  [2.7.6]. 


FIGURE  42.  Spatial  Frequency  Channel  Model 
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At  the  same  time,  (CS),  covered  in  section  [34],  is  calculated  as  a  function  of  eccentricity  and  spa¬ 
tial  frequency,  in  a  similar  manner  to  the  one  described  in  the  VEP  model.  CS  is  used  as  a  measure 
of  the  gain  of  each  channel  and  is  multiplied  by  the  rectified  channel  output. 
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FIGURE  43.  Overview  of  Calculations-Spatlal  Frequency  Channel  Model 


CS  is  used  here  is  a  simple  means  of  calculating  the  sensitivity  of  each  channel.  To  obtain  the  CS 
at  any  desired  values  of  SF  and  eccentricity,  interpolation  must  be  used.  This  is  described  below  in 
section  [3.3].  Therefore,  the  model  shows  a  CS  calculation  routine.  One  last  routine  calculates 
interactions  between  the  given  channel  and  others.  This  is  shown  as  being  multiplied  by  the  out¬ 
puts,  so  a  value  of  1  indicates  no  interaction  effect.  This  routine  would  have  been  added  if  neces¬ 
sary,  after  completing  the  rest  of  the  model. 
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FIGURE  44.  Post-Channel  Calculations. 


The  outputs  of  the  channels  will  be  sent  back  to  the  d’  calculator  in  figure  [42],  One  way  to  con¬ 
vert  the  outputs  to  d’  is  shown  in  equation  [71],  This  is  based  on  the  Quick  Pooling  formula  of 
equation  [46].  Equation  [71]  has  an  exponent,  q,  as  well  as  a  function  f.  f  would  have  been 
assumed  to  be  multiplication  by  an  adjustable  weight,  that  weight  and  q  would  have  been  set  by 
fitting  against  experimental  data. 
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To  utilize  this  model  in  a  practical  manner  to  simulate  vision  loss  due  to  retinal  damage,  the  func¬ 
tions  and  algorithms  taking  place  in  the  various  boxes  must  be  specified.  This  can  be  done  by  cal¬ 
culating  d’  values  for  a  series  of  psychophysical  experiments  being  run  and  perhaps  also  from 
experiments  in  the  literature.  Certain  effects,  such  as  the  early  non-linearities  shown  in  figure  [42] 
can  probably  be  safely  ignored,  assuming  linearity  for  small  ranges  of  illumination  changes.  Sim¬ 
ple  grating  experiments  should  excite  only  one  channel  at  a  time,  and  be  easy  to  simulate  in  the 
model.  More  complicated  patterns  could  be  run  exciting  multiple  channels,  allowing  development 
of  more  complex  aspects. 
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V.  Conclusions 


While  the  research  is  still  on  going,  the  IRC  model  shows  great  promise.  It  has  shown  good  results 
fitting  CSF  and  cycle  summation  curves.  The  latter  is  an  important  result  because  it  shows  that 
cycle  summation  is  not  due  to  some  mysterious  summation  rule  or  response  of  spatial  channel:  it 
is  due  to  the  same  d’  summation  rule  for  any  other  pattern.  The  premise  is  tantalizing  -  detection 
can  be  modeled  by  a  summation  of  information. 

The  spatial  frequency  channel  model  as  described  would  have  been  easy  to  implement  -  it  is  con¬ 
volutions,  CS  interpolations,  and  some  other  calculations,  all  of  them  fairly  simple.  The  model 
would  have  suffered  from  at  least  one  serious  flaw,  using  CS  functions  to  calibrate  the  channels 
might  work  -  as  long  as  the  same  kind  of  stimulus  was  used  to  generate  the  CS  curves  as  to  gener¬ 
ate  the  stimulus  under  test.  It  would  lack  the  flexibility  of  the  IRC  model  which  can  have  its 
parameters  retrained  for  different  kinds  of  stimuli.In  retrospect,  it  has  been  realized  that  this 
implementation  of  a  spatial  frequency  channel  model  does  not  simulate  the  chaimels  very  well.  It 
would  not  handle  the  properties  described  above  in  section[2.7.7],  because  it  does  not  implement 
charmel  masking  or  adaption.  Therefore,  it  does  not  seem  to  be  a  useful  exercise  to  implement  this 
model  without  major  design  changes. 

The  airplane  results  are  somewhat  contradictory.  One  subject  shows  excellent  agrement,  the  other 
is  about  .3  log  unit  high. 

How  far  will  the  IRC  model  take  us?  Color  effects  could  by  added  by  simulating  the  positions  of 
the  three  kinds  of  cones  and  their  connections  to  ganglion  cells.  Temporal  effects  could  be  added 
by  simulating  the  delay  between  surround  and  center.  This  should  alleviate  the  need  for  using 
parameter  fitting  to  set  the  ratio  of  weights  between  the  center  and  the  surround.  Tuned  and  ori¬ 
ented  cortical  cells  can  be  added  in.  The  computational  requirements  of  such  a  model  would  be 
considerable. 

Of  course,  the  purpose  of  this  project  is  the  simulation  of  retinal  damage.  The  results  are  inconclu¬ 
sive,  as  yet,  but  show  excellent  promise. 
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Appendix  III 


The  Effects  of  Disease  on  Signal  Processing 
Characteristics  of  the  Human  Visual  System 


TO  LASER  DAMAGE" 

M.  Dcsai 


Wc  have  developed  a  .signal  processing  model  for  ihe  effect  of  la.ser  damage  on  visual  acuity. 
Our  goal  is  to  explore  Iiow  laser  damage  over  a  small  region  of  the  retina  (e.g.’at  or  near  the  fovea) 
allects  vision  at  other  points  on  the  retina.  We  believe  that  the  destruction  of  part  of  the  retina  may 
alleet  acuity  cl.sewhere.  The  damaged  area  is  modeled  to  have  complete  lo.ss  of  vision,  (i.e,  zero 
sen.sory  reception  or  dead  /.one).  We  then  model  proce.ssing  of  a  .stimulus  with  the  la.ser-damaged 
letina  using  a  linear  I  liter  model  of  the  neural  receptive  fields.  The  receptive  field  models  are 
Incali/.ed,  hence  the  effects  ol'  processing  are  al.so  localized.  However,  it  can  be  argued  that  the 
ellects  ol  the  la.ser  damage  ellect  the  respon.ses  of  neighboring  receptive  fields,  hence  propagating 
the  ellecLs  ol  the  damage  Irom  the  damaged  area. 

Vi,„„l  Kuity  i,.  l.y  ll,c  O  ll,c  la,.cr.duma|;c  model  l„  dilTerciil  inpul.,,  A 

.,iinple  mea.,„rc  of  rc.„>l,nion  i.,  made  n.,i„g  simple  point  liiriiel.,  (impulses)  of. small  .separation  as  the 
mpu.  to  the  .systent.  The  point  la, pels  a,e  then  passed  thronph  appiopriate  linear  niters,  inelutlinp 
Gaus.s,a„  a,,d  Gahor  n„e,-s.  A  second  nteas.ne  of  acuity  is  ha.sed  on  texture  diseriminahility.  In  this 
has  •d*'"^ '  ll  to  be  a  .sum  of  two  simi.soidal  grating.s.  The  measure  of  acuity  is 

region."”  -separation)  of  the  gratings  in  the  vicinity  of  the  damaged 


Using  these  models,  we  are  able  to  show  that  the  la.ser  damage  reduces  visual  acuity  near  the 
-cad  xonc.  the  lollow.np  sense:  the  lihered  stintul,  (neuronal  responses,  a,-e  alTected  n  s 

•1  .  data  ,s  u.sed  hy  suh,set|nenl  pr„ee.s,sing  elements  may  elleel  proees.ses  such  as 

encops,,  an  mot, on  lield  pereeption  hod.  „„l,i„  Ihc  damaged  region  (where  i,  is  eradicated,,  and  il 
.  V  cm.ty  of  the  da.naped  ,-epion  as  well.  The  extent  of  the  effect  is  found  to  depend  di,ec,ly  „„  ,|,c 
Glee  ,ve  .span  (spat.al  exten.)  of  the  channel  or  teceplive  field.  This  implies  tha,  large  peripheral 
epttve  ,c  ds  tnvolvmp  a  signiltcan,  (spatial,  degree  of  .spatial  summation  over  the  relina  may  he 


Wc  are  al.so  aide  to  show  that  the  effect  is  wottsened  as  the  amhient  light  intensity  increase, 

eftos  r"  ™^"y'GS  the 

'■  ol  varying  s„c|,  parameters  as  the  .spatial  extent  of  the  teception  Held  (or  the  filter  handwidth) 


and  lliccxlcnl  ol  (he  dama 


A  more  extensive  analysis  is  then  developed  (hat  attempts  to  analy/.e  the  cffecLs  of  damage  on 
a  mulli-ehanncl  (wavelet-like)  visual  processing  system.  This  analysis  encapsulates  the  effects  of 
vaiiahle  resolution  (sampling  density)  over  the  retina,  as  well  as  an  analysis  of  how  the  errors  arc 
distnbuted  across  channels  (hence  scales)  in  a  multi-channel  scenario.  Finally,  observations  arc  made 
on  how  extensively  certain  basic  low-level  visual  tasks  may  be  effected  by  (he  la.scr  damage.  Tliis  is 

^K-compli.shed  by  examining  the  effect  of  the  Hawed,  filtered  stimulus  on  computer  models  of.  c.g., 
stcieo,  textuie,  and  motion  processing 


/.  1-1)  Sif^nal  Proccssinj^  iMixlcl 


()i.r  goal  m  this  one-dimensional  model  is  to  explore  how  la.ser  damage  over  some  region  of 
‘|•'"la  (e.g.  lovea)  allects  vision  at  other  points  on  (he  retina.  In  other  word.s.  how  does  the 
c  e.struction  ol  part  of  the  retina  may  alTeet  acuity  el.sewhere.  To  under.sland  this,  llrsl  wc  consider  the 
.simple  1-D  model,  that  is  equivalent  to  considering  the  damage  along  a  .single  line  as  illu.strated 
o.  this  analysis  assume  continuous  mathematics  and  uniform  sampling  The 
imiltire.solution.analysis  will  be  atidre.s.sed  in  the  later  .sections. 


point  targets 


Fovea  with  la.ser  damaged  illustrated  r 


as  point  targets. 
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Onc-dimcn.sional  1  liter  inode!  with  [t,  to  I2]  being  the  dead  7, one  or  the  laser  damaged  area 


This  illustrates  the  model.  First,  eliminate  a  chunk  of  a  signal,  then  pa.ss  the  .signal  through  a 
filter.  The  result  (output)  will  he  changed  (differ  from  the  ideal)  at  places  other  than  in  the  “dead” 
7.nne  [tj.  ij].  In  fact,  it  may  be  changed  everywhere  (although  likely  negligibly  with  .some  di, stance 
liom  the  damaged  point).  The  “lilier.s”  u.sed  may  be  thought  of  as  models  for  neural  receptive  fields, 
which  are  l()cali7.ed,  so  the  effect  will  be  locali/x:d,  yet  WILL  propagate  from  the  damage  area. 

Within  a  measurable  vicinity  of  the  damaged  region,  there  will  be  pronounced  effect  on  the  Filler 
responses. 
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Now  consider  an  ideal  image  signal  (visual  slimuh.s)  x(0.  Denote  the  signal  resulting  from 

retinal  damage  to  he  x'(t).  where  the  damage  results  in  loss  of  signal  in  the  interval  [l,.  hi.  We  can 
then  write: 


xu;  =  d(t-i3)  +  5(l-g  +  C 


wtoe  6(1)  i.  ,1,.  Dirac  impulse  luncii,,,,.  We  win  he  using  i,',s  in.egral  prnpenies  (in  cunvolu.ions)  ,o 
mrrdel  ,hc  respense  u,  clesely-placed  pein,  iargecs.  This  signal  is  depicted: 


The  constant  level  C  which  ntodels  a  non-veto  hackground  (constant  amhienl  light  Hux)  is  both 
ccessaD  and  real, she.  Any  realistic  .scene  has  a  background  illuminalion.  Moneover.  without  it,  (he 
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dead  zone  will  nol  show  up!  (this  has  an  important  implication). 
The  ideal  response  is: 


y(0  =  g(t  - 13)  +  g(t  - 14)  +  Cj  g(t)  dl. 

Make  the  simple  a.ssumplion  that  g  is  Gaussian  (could  be  Gabor  etc).  Then  we  may  model  the 
resolution  loss  at  the  gaussian  channel  output  (again,  no  damage  yet).  The  following  depicts  the 
cla.ssic  acuity  problem  observed  at  the  filter  output: 


Now  we  consider  the  case  ol  a  la.scr-damaged  retina.  The  si^-nal  with  damage  is 
x'(0  =  (5(t  -  h)  +  5(1  -  I4)  +  CJ  [u(t,-l)  +  u(l-  tj)] 

=  5(t  - 1^)  +  5(1  - 14)  +  C-[u(t,-t)  +  u(i-  tj)] 

(of  course  t,,  I4  arc  not  in  the  dead  zone).  This  is  depicted  below; 


The  output  is  now; 


y  (0  -  ti)  +  g(t  -  (4)+ C{  r  g(t)dt+  r  "g(t)dl} 

Thu,,  we  can  model  U,c  rc.,„ln,ion  lo.s,  a,  ,l.c  sau,,,,ia„  channel  m„pu,  for  .he  ca,,e  where  .here  has 
no  damage''"  ““  «l'cm  there  is 


WITHOUT  DAMAGE 


WITH  DAMAGE  -  RESOLUTION  REDUCED 


The  Gaussians  Slili  mix,  bm  aciily  i.,  I'unhcr  reduced  by  ihe  proximity  of  the  response  to 
the  dead  zone.  Tins  lo.ss  in  resolution  is  manifested  by  a  lessening  of  one  of  the  gaussian  peak 
rc.sponses,  m-  ing  it  difficult  (for  example,  in  an  algorithm,  or  a  neural  network),  to  discriminate 
one  peak  from  another.  It  could  easily,  in  fact,  be  interpreted  as  a  single  peak. 

As  might  be  expected,  the  effect  is  worsened  as  the  light  ambient  intensity  increa,scs  -  the 
trough  IS  deepened,  thus  reducing  the  resolution  further; 
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pic/ihrugh  .‘'fr"iTr,"prrc“r  wm  r  r  '"li''’"  ^ 

■:).  The  “ni,ers"  arc  model,  for  neural  receplive  0^'  wh'”r  "" 

localized,  yet  WILL  nron-M’-iii'  (V  ik  «  i.  .  ’  localized,  so  the  effect  will  be 

depends  directly  on  the  effective  s'Zn  ("'  T\^''  significant  way.  The  extent  of  the  effect 

- 

arreeled.  The  cffec.  win  he  more  pronounced  when  ,l.c  amhieni  li,h.  LZ,y  1  h"l'.I“' 

2.  Multi-Channel  Model 


e  have  performed  a  preliminary  analysis  of  the  effects  of  hcpr  i 

or  filters.  The  multi-channels  mav  b(' •)  w  i  ix  r-  on  multiple  channels 

of  intcres.  consider  ,„e  e'Zs  1  d"  "  “  "  '> 

retina.  In  parlieular,  il  would  be  ikcIuI  h' ''  resolutions  over  the 

how  well  visual  tasks  can  he  t.cc„mplishel,sI'Ja,e"l!LTar' 
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x{t) 
x'(t  ) 


Experimental  Results: 


•  V  t  ■ ).  i.u,  ihc  Lena  image  reproduced  in  Figure  la. 

bandwid^roretri!!re  90°.  The 

r,.n.  f  •  ac^cordance  wilh  typical  Gabor-shaped  receptive  fields  THp 

n..r  frequences  vary  non,  ,„eyelesA.mase,„  UX)  cycle.s/in.asc  in  octave  incremenrs.  ' 

ihe  image  i  e  with  the  im-u,  ^  ^  with  a  square  “hole”  in 

.■or.hislpur;sat  :  r 

nuer  i.ea,  an.  fritere.  .a.a.e:' il^lnirr:  IZtft'rrT  ^T" 

r. ::  ““  -  --- 
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:  :p 


Original  Lenna 


Figure  1:a.  Unfiltered  256  x  256  Lenna  images  with  and  without  20 


X  20  mask. 
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Figure  1  ;c. 
difference. 


Unmasked  Lenna:  Channel  2 


Difference  image 


Channel  2  outputs  of  Gabor  filtered  images  and  their  absoiute 
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Figure  1;d. 
difference. 


Unmasked  Lenna;  Channel  3 


Channel  3  oulpuls  of  Gabor  fiiiered  images  and  their  absolute 
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Unmasked  Lenna:  Channel  5 


Masked  Lenna:  Channel  5 


Difference  image 


Figure  1  ;f.  Channel  5  outputs  of  Gabor  filtered 
difference. 


images  and  their  absolute 
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HUMAN  VISUAL  PROCESSING  OF  IMPULSIVE  NOISE 


James  R.  Dykes,  Jr. 
Susmitha  Baddipudi 


The  University  of  Texas  at  San  Antonio 
San  Antonio,  Texas  78285 


Abstract  of  Psychophysical  Experiments 


Any  iraiima  (such  as  iniravitrcal  blood  resulting  from  laser  damage)  which  adds  noise  to  the 
retinal  image  will  degrade  performance,  but  the  manner  in  which  the  human  visual  system  codes  the 
information  available  to  the  photoreceptors  will  have  a  profound  effect  in  how  resistant  performance  is  to 
that  noise.  Two  general  classes  of  models  of  this  coding  have  been  proposed.  The  first  (including 
threshold  theory  and  median  filtering)  predicts  minimal  effects  of  low  Ics  cls  of  noise,  but  p  discontinuity 
in  the  function  at  the  “break  point”.  The  second  (including  SDT  and  the  logistic)  predict  no  such 
disconitinuit)';  performance  will  gradually  decline  as  more  noise  is  added.  While  previous  researchers 
have  voiced  opinions  based  on  the  steepness  of  the  function,  those  data  are  insufficient  to  distinguish 
between  the  two  models.  The  current  research  had  observers  discriminate  between  a  pair  of  binary  target 
stimuli  as  increasing  levels  of  binary  noise  were  XORcd  w  ith  the  stimulus.  Extending  previous  research, 
two  additional  independent  variables  were  added;  the  number  of  distracters  (0,  1,  3,  or  15)  and  the  nature 
of  the  distracters  (blanks  or  hjbrids).  A  pilot  study  was  used  to  choose  the  appropriate  four  levels  of  noise 
and  to  refine  the  three  pairs  of  targets.  While  the  results  of  the  main  cx-periment  arc  limited  to  the  type  of 
stimuli  and  noise  chosen,  the  data  were  belter  fit  by  the  SDT/logislic  model  than  by  the  Ihrcshold/median 
filtering  model.  To  the  extent  dial  diffuse  iniravitrcal  blood  is  appropriately  modeled  by  binary  impulsive 
noise,  one  can  predict  that  increasing  levels  of  noise  will  lead  to  a  gradual  and  graceful  decline  in 
performance  as  opposed  to  level  performance  with  a  sudden  drop. 


Introduction 

Many  traumas  to  the  eye  can  produce  a  hemorrhage  causing  blood  diffusion  in  the  vitreous 
humor  of  the  eye.  This  in  turn  produces  noisy  images  on  the  retina  of  the  eye.  One  of  the  advantages  of 
using  a  median  filter  in  image  processing  is  that  it  is  resistant  to  noise.  As  the  anwont  of  noise  added  to  an 
image  is  mcnased,  the  output  of  a  median  filter  convolved  with  the  image  remains  high  in  quality  up  to  a  fairly 
high  level  of  noise.  Beyond  that  break  point,  the  output  quality  degrades  precipitously.  If  the  human  visual 
system  included  a  median  filler,  our  perception  should  remain  stable  over  low  levels  of  noise  and  then  should 
drop  rapidly  past  tlic  break  point.  This  prediction  is  consistent  with  classic  llueshold  theory.  When  image 
intensity,  duration,  and  quality  arc  above  threshold,  performance  should  be  almost  perfect.  Below  threshold, 
performance  should  be  at  chance.  Both  using  a  median  filter  and  threshold  theory  predict  tliat  the  performance 
function  should  have  a  discontinuity  at  tiic  break  point. 

In  contrast,  the  human  visual  sy  stem  might  cmplo)-  operators  that  produce  performance  functions  with 
no  such  discontinuity.  One  such  operator  (die  logistic)  has  been  widely  cmplojcd  in  PDF  models  of  visual 
function.  Assuming  the  logistic,  performance  would  degrade  gradually  and  gracefully  as  a  function  of  noise 
level.  Such  predictions  are  also  consistent  widi  Signal  Detection  Theory  (SDT).  According  to  SDT,  the  neural 
code  IS  a  continuous  function  (albeit  digitj7ed  b>-  neurons)  of  the  input.  There  is  no  threshold  or  di^ntinuity. 
Although  signal  detection  theory  was  developed  for  detection  problems,  it  can  be  extended  to 
disenmination  situation.  For  this  c.xpcrimcnl,  the  observer  has  to  say  which  of  two  stimuli  was  presented 
on  a  given  trial.  This  is  called  a  two-allcrnalivc  forccd-chicc  task  in  the  psychophysical  lilcralurc. 

The  cfTTccts  of  adding  noise  to  an  image  has  been  an  important  research  question.  As  one  e.xamplc. 
Barber  &  dc  la  Mahoticrc  (1982)  conducted  an  experiment  which  measured  the  case  with  which  words, 
degraded  by  noise,  could  be  recognized.  More  specifically,  the  motivation  for  this  research  came  from ' 
data  reported  by  Thompson  &  Slocsscl  (1987).  They  measured  their  drop  in  acuity  as  a  function  of 
viewing  the  Snellen  charts  through  different  dilutions  of  human  blood.  Blood  from  a  normal  subject  was 
diluted  with  saline  from  1/1(X)  to  1/12,800.  A  Iota!  of  seventeen  concentrations  were  prepared.  The  two 
opihamologisls  tried  to  read  the  Snellen  chart  through  a  4mm  aperture  and  a  c)linder  containing  the 
blood  solution.  Visual  acuity  was  20/20  when  tested  through  normal  saline  and  all  tests  with  blood 
dilutions  of  1/4,800  or  higher.  Visual  acuity  decreased  with  increasing  concentrations  of  blood  between 
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the  dilutions  of  1/3.200  and  I/IOO.  At  dilutions  of  1/400  of  less,  visual  acuity  reduced  to  hand  motions. 
Thus  between  dilutions  of  1/4,800  to  1/400,  acuity  dropped  from  20/20  to  hand  motions.  This  “rapid" 
change  in  performance  through  a  greater  than  10  fold  change  in  concentration  suggested  that  the  visual 
system  might  be  modeled  by  a  median  filter.  While  vaiyng  blood  corxxtration  is  the  most  direct  way  to 
mimic  the  density  of  an  intravitreal  hemorraghe,  it  might  be  possible  to  model  these  hemorraghes  varying  the 
level  of  binary  impulsive  noise  replacing  image  pi.xcls.  This  prosidcs  nrorc  trials  per  subject  and  generates  a 
broader  base  of  observen. 

Whether  using  concentrations  of  blood  or  implusivc  noise,  it  is  not  possible  to  distinguish  between 
these  two  nxxlcls  (threshold  and  SDTj  and  operators  (median  filter  arxl  logistic)  by  simply  looking  at 
pcrformaiKe  as  a  function  of  blodd  concentration  /  noise  level.  Since  there  is  always  some  error  variance  in  the 
data,  the  performance  function  would  not  contain  a  clear  discontinuity  at  the  break  point  even  if  the  human 
visual  system  djd  use  a  OKdian  filter  and/or  thresliold.  Simply  looking  at  the  performance  function  leaves  one 
asking  tlK  unanswerable  question;  "How  steep  is  steep?". 

(Tn  tire  other  hand,  it  might  be  possible  to  distinguish  between  the  two  models  if  the  data  were  reliable 
and  if  a  second  independent  variable  (in  addition  to  noise  level)  were  manipulated  that  made  diflcientiable 
predictions  between  the  two  models.  The  independent  variable  chosen  for  this  study  was  number  of  possible 
locations  (1, 2,  4,  or  16).  This  variable  has  an  established  research  history  within  SDT  and  the  two  models  make 
quite  different  predictions.  Burgess  (1985(a).  1985(b),  Burgess  &  Colbume  (1988),  Burgess  & 
Ghandcharian  (1984))  has  demonstrated  the  theoretical  importance  of  this  variable  in  signal  detection.  In 
order  to  bracket  performance,  tire  distractors  used  were  either  blank  cells  with  noise  added  or  hybrids  with  noise 
added 


Pilot  Experiment 


The  aim  of  the  pilot  e.xpcnmcnt  was  to  find  a  set  of  stimulus  pairs,  hybrids,  and  noise  levels  that 
generally  produced  performance  from  close  to  perfect  (for  1  alternative  and  low  noise  level)  to  chance  (for 

15  hybrid  distractors  and  high  noise  level).  Independent  noise  sample  replaced  pi.xcl  element  in  the  target 
and  distractor  locations.  ° 

Method 


Subjects 

Ninety  si.x  students  from  University  of  Te.xas  at  San  Antonio  parlicipaicd  in  the  experiment.  They 
were  given  credit  for  partial  fulfillment  of  the  requirements  for  their  freshman  year.  Visual  acuity  in  all 
subjects  was  normal  or  corrected  to  normal. 


Design  and  Procedures 

Tw  1  K  "  conducted  in  the  Psychology  lab  at  the  University  of  Te.xas  at  San  Antonio 

This  lab  had  six  rooms  and  each  room  had  an  IBM  PC  done  with  an  Intel  286  processor  and  mother 

Ti  T""'  (preventing  the  subjects  from  being 

distract^).  The  observers  were  led  into  the  room  which  contained  the  computer  and  were  asked  to  sit  on 

a  chair  facing  the  screen.  The  stimuli  were  presented  at  the  center  of  the  monitor.  Each  stimulus  display 

-/ "itvwT  ?'!!!  Observers  were  instnictcd  to  place  their  index  fingers  on  the  "z"  and 

/  keys  of  the  keyboard.  The  two  possible  targets  were  always  present  on  the  lower  Icfi  hand  and  right 
hand  comers  of  the  screen.  The  "z"  key  corresponded  to  the  target  on  the  left  hand  corner  of  the  screen 
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and  the  r  key  cxirrcsponded  to  the  right  hand  corner  of  the  screen.  For  a  gj\’en  stimulus  pair,  the 
position  of  the  two  possible  Urgets  was  counterbalanced  across  observers.  After  deciding  which  target 
was  present,  the  observer  pressed  the  appropriate  kc>'.  Stimuli  were  response  terminated.  Feedback  rvas 
provided  with  two  different  kinds  of  sounds:  a  "chirp"  for  the  correct  response  and  a  "buzzer"  for  the 
incorrect  response. 

Five  independent  variables  were  manipulated  in  the  e.xpcrimcnt:  (i)  stimulus  pair,  (ii)  Ijpc  of 
distraclor,  (iii)  percentage  of  noise  pixels,  (iv)  number  of  alternative  positions/distractors,  and  (v)  location 
of  the  target. 

Three  stimulus  pairs  were  chosen.  A  particular  observer  was  presented  with  either  two 
alphabetic  stimuli  (T  and  F),  a  pair  of  airplanes  inclined  at  45"  and  135"  angles,  or  a  pair  of  horizon 
angles  (  45  or  135  ).  The  three  sets  of  stimuli  arc  shown  below,  T  and  F  show  high  confusability 
(Gilmore  et  al.,  1979).  The  windows  containing  the  stimuli  were  16  pixels  wide  and  18  pixels  high 
including  a  4  pixel  border.  The  windows  consisted  of  black  pixels  (value  0)  against  a  white  background 
(value  255). 


Two  tjpies  of  disiraclors  were  cho.scn;  blank  windows  and  windows  with  hj-brids.  Blank  windows 
were  white  with  noise  superimposed  on  them  Hybrids  resembled  both  the  alternatives  and  were 
composed  of  parts  of  each.  The  same  noise  level  (but  independent  samples)  replaced  pixels  in  each 
distractor  as  in  the  target.  The  distractors  al.so  were  the  same  size  as  the  target. 


Given  binary  images,  impulsive  noise  involved  XORing  individual  noise  elements  with  image 
pixels.  This  IS  called  "salt  and  pepper",  because  XORing  flips  the  image  pixels:  a  white  pixel  becomes 
black  and  a  black  pi.xel  becomes  white.  The  proportion  of  image  elements  XORed  was  independenUy  but 
uniformly  distnbul^  across  each  line  of  the  stimulus  w  indow.  The  proportion  of  pi.xels  XORed  can 
.  wen  0  and  .5.  Between  .5  and  1,  the  "noisy"  output  becomes  more  recognizable  but  reverse  in  binary 
'^h  bridTr  "";‘’™'‘°ns  of  the  first  three  independent  values  ( each  of  the  three  stimulus 

A  .'l'”  7  "  °  f,g^.  Four  diirc^nl  M  of' 


impulsive  noise  {0.125.  0.188,  0.25.  and  0.313}  were  used  in  the  pilot . 

f  of  windows  or  possible  locations  be  denoted  by  the  letter  "m".  Four  different 

distractors  wa'sT  ^  '3^M5  I'  was  always  present,  thus  the  number  of 

0  f  l^or  Is  T  "1"  •"  =  1  one  of  the  windows  numbered 

0.  3.  12.  or  15  was  presented  on  the  screen  and  this  window  contained  the  target.  For  m  =  2.  two  windo^ 
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total  of  576  trials  broken  intoTblc^Ts  c^siqin'pTf '8^^"'^  T  ^istracicr.  Each  was  presented  with  a 

.hrcc  of  ,„c  a.™,c  hind  ap^,,:d“  a  ““  “  ^ 

\  A-  1 , 1 , 1  ,r  .r .  1 . ),  no  more  than  three  noise 
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lr^'clsofsamcv'alueappcaredinarow(cx:0.125,  0.125,  0.250.  0.188,  0  188  0  188  0  125  )  and  no 

more  than  three  trials  in  a  row  had  the  same  "m"  value  (ex:  tn=2.  m=2.  m=2'.  ni=l.’m=l,  m’=4'"..’.  ) 

The  allcmatives.  noise  Ics'els  and  m  size,  occurred  equally  oflen  within  each  block.  ’  All  lotions  (crossil 
With  Lhc  other  independent  variabics)  occurred  equally  often  across  the  8  blocks. 


Results  ancJ  Discussion 

For  each  stimulus  type,  the  percentage  of  correct  responses  are  displayed  as  a  function  of  type  of 
distractcrs^lank  versus  hybrid),  noise  level  (.125.  .188,  .250.  and  .313),  and  m  size  in  the  following 
paphs.  While  performance  generally  declined  as  a  function  of  noise  level  and  m  size,  the  data  arc  noisy 
In  spite  of  these  imperfections,  the  pilot  study  served  its  function.  The  two  types  of  distracters,  the  four 
noise  levels,  and  four  m  sizes  clearly  bracketed  performance.  On  the  other  hand  the  data  for  the  horizon 
angles  were  lowest  and  least  regular.  Consequently,  the  general  design  of  the  blocks  and  the  angle  and 
airplane  stimuli  were  mexiifed  for  (he  rnnin  experiment 


Main  Experiment 


Method 


Design  and  Procedures 

The  number  of  trials  in  each  block  of  the  pilot  was  too  long.  The  attention  span  of  the  subject 
diminished  towards  Uic  end  of  the  block.  The  number  of  limes  a  target  occurred  at  each  location  was  not 
sulTicicm  to  generate  stable  data  points.  Also  the  rapid  changes  of  m  sizes  within  a  block  caused  the 
resulu  to  ^  noi$>'  for  higher  m  sizes.  Tnese  problems  were  corrected  in  the  main  experiment  The 
num^r  of  Inals  in  each  block  was  changed  to  48  as  opposed  to  64  in  the  pilot  in  order  to  better  match  the 
at  cntion  span  of  the  subject.  The  number  of  blocks  was  increased  to  16  (compared  to  eight  blocks  in  the 
pilot  c.xpenmcnt).  Increasing  the  total  number  of  trials  to  768  (compared  to  576  in  the  pilot)  provided 
more  trials  for  each  target  at  each  location  with  each  level  of  noise  and  m  size.  This  should  provide  more 
^ble  wtimates  of  accuracy.  Each  block  consisted  of  only  one  m-sizc  to  also  increase  stability  in  (he  data. 

order  of  m  sizes  across  the  16  blocks  was  counter-balanced  using  a  digram-balanced  Latin  Square 
Each  block  was  preceded  by  the  observer  completing  two  perfect  runs  of  12  practice  trials.  These  trials 
were  the  same  m  size  as  in  the  test  trials  for  that  block,  but  contained  no  noise.  These  practice  runs  were 
designed  to  increase  the  onservers’  familiarity  with  the  task  and  m  size. 


nen  the  small  size  of  the  displays,  the  angles  and  airplanes  banking  at  45  degrees  and  135  degrees 
looked  N  cr)'  similar  and  this  caused  the  lowest  acairacy.  The  planes  were  shifted  to  vertical  and 


honzontjl  in  the  main  study  to  maximize  the  use  of  the  pixels  available.  The  horizon  angles  were 
changed  to  geometric  figures.  Tlicsc  are  shown  in  the  figure. 

Subjects 


One  hundred  and  forty  students  from  the  University  of  Texas  at  San  Antonio  fwrticipatcd  in  the 
experiment.  They  were  given  credit  for  partial  fulftilmcnt  of  the  requirements  for  their  freshman  year. 
Visual  acuity  in  all  subjects  was  normal  or  corrected  to  normal.  The  dau  for  44  of  these  observers  w-erc 
discarded  due  to  poor  baseline  performance  or  extreme  response  bias.  To  be  included  tire  observer  had  to 
have  at  least  50%  correct  responses  in  the  lowest  noise  condition  (.125).  No  effect  of  noise  could  be 
obtained  if  the  observer  was  at  chance  at  the  lowest  level.  Secondly  the  subject  had  to  have  at  least  20% 
responses  for  each  target  in  each  condition.  Since  the  targets  were  cquiprobabic,  a  more  extreme  than 
80%  /  20%  split  indicated  a  very  extreme  response  bias  and  suggested  that  the  observer  may  have  just 
repeatedly  pressed  the  same  button  without  attending  to  the  stimuli.  AH  but  5  of  the  observers  dropped 
were  given  the  airplane  targets.  The  remaining  96  observers  consisted  of  16  people  given  each  of  three 
stimulus  sets  by  two  distractor  combinations. 

Results 


For  each  subject,  the  proportion  of  correct  responses  was  computed  for  each  m  size  at  each  noise 
level.  As  can  be  seen  in  the  following  figures,  the  dramatic  drops  in  performance  as  a  function  of  noise 
lesel  and  m  size  were  more  stable  than  in  the  pilot.  On  the  other  hand,  the  new  airplane  stimuli  were  still 
very  dimcult  to  discriminate  (especially  from  the  hjbrid  distractors).  Even  for  the  subjects  meeting  the 
criteria  listed  above,  overall  performance  for  the  airplanes  was  low  and  unstable. 


Ihc  correct  location  of  the  stimulus  as  a  function  of  m  size  (d’m).  Parallel  to  the  assumptions  for  the  Luce 
model.  It  IS  assumed  that  the  observer  knows  the  location  (based  on  d’m)  and  is  accurate  based  on 
stimulus  discriminability  (d’)  or  thinks  another  location  contains  the  target  (also  based  on  d’m)  and 
posses  at  chance  (.5  in  the  2AFC  task).  Given  the  assumption  that  tlie  SDT  distributions  are  guassian,  it 
IS  ncces^ry  to  scale  proportion  correct.  A  similar  transformation  occurs  with  the  logistic.  The  proportion 
correct  for  a  given  d’  and  m  size  can  be  defined  as 

Pic)d,  =  [(p(c)/7/)  *  /?(c)^_  ]  +  [.5  ( 1  -  p{c)m] . 

where 


UKc)n\)  =  the  propbability  of  knowing  the  target  location  based  on  d'm 


probability  of  knowing  a  single  target  based  on  d 

Accuracy  for  m=  1  can  be  used  to  measure  d’.  d’  and  the  proportion  correct  for  m=  2  can  be  used  to 
determine  d  m.  As  with  the  Luce  model,  those  can  then  be  used  to  predict  accuracy  for  m=  4  and  16. 

chm  P'^'^'crions  for  each  model  (Luce  and  SDT)  were  separately  computed  for  each  observer 

s  own  each  stmulus  tyitc  (airplane.  T/F.  and  geometric)  at  each  noise  level  (.125.  188  250  and  3 13) 

...  .  locations.  It  should  be  noted  that  as  performance  in  the  m=  1  and  2  conditions 

pSSiiTfet,"  of  “"verge  lowerd  e  eommon 

Lr. '”spt.r„r  Lde  ri-rSre:^-  "  o“">' '« 

airplanes  at  “S 

r.te't  r  r  r  t'”"  '  ™ 

rOservers  are  shot  on  “  f*"  ““  "ff"" 

drop  in  ptotfntr L^n^o^r 

trr  > 

nerf.rm.*"  stimuli/hybrids  Condition 

pcifcct(l.O)  to  chance  (0.50)  almost 

Luce  vers^us'sDT)  X  4  (ToS^^IT"  iT'^^or  t D  X  ^  ^  (seometric  versus  TF)  X  2  (modcN 

to  sec  if  the  predictions  of  the  Luce  or  SDT  mrvtrv.  L  V*"  i  **  ANOVA  was  performed 

^rand  m'T  rct-eTm"  ?7‘"p< 
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The  main  goal  of  the  experiment  was  to  find  out  whether  visual  performance  with  impulsive 
noise  is  better  modeled  by  a  discontinuous  function  (as  predicted  by  threshold  theory  and  median  filters) 
or  by  a  smooth  cur\'c  with  no  break  point  (as  predicted  by  SDT  theory  and  the  logistic).  For  a  given 
subject  (corresponding  to  a  single  stimulus  and  distractor  type  combination),  it  is  possible  to  use  the  data 
for  m=  1  and  m=  2  for  a  given  noise  level  to  predict  what  performance  should  be  found  for  m=  4  arrd  m= 
16  at  that  noise  level.  The  models  make  different  predictions  for  those  m=  4  and  m=  16  conditions. 

Luce  proposed  a  model  of  threshold  theory.  On  a  given  trial,  the  stimulus  quality  is  either  above 
or  below  threshold.  In  terms  of  a  median  filter,  the  noise  level  of  the  image  is  either  above  or  below  the 
break  point.  When  the  stimulus  on  a  given  trial  is  above  threshold,  the  observer  is  certain  and 
performance  will  be  perfect  on  that  trial.  When  the  stimulus  on  a  given  trial  is  below  threshold,  the 
obsciver  will  perform  randomly  (at  chance)  on  that  trial.  At  the  threshold  (break  point),  performance 
should  rapidly  shifi  from  perfect  to  chance.  For  this  experiment.  Luce’s  model  needs  to  be  extended  in 
order  to  deal  with  the  distraclors.  To  be  sure  of  the  target  on  a  given  trial,  the  target  must  be  above 
threshold  and  all  of  the  distractors  must  be  below  threshold  for  the  incorrect  response.  Otherwise  the 
observer  must  guess.  Thus  two  probablities  must  be  known:  p  is  the  probability  that  the  target  is  above 
threshold  and  q  is  the  probability  that  a  distractor  is  below  threshold.  For  a  given  noise  level,  the 
proportion  correct  for  a  given  m  size  can  be  defined  as 

p{c)m  =  l[/7  -  (p  * 

When  sure,  performance  will  be  perfect  (1).  Otherwise  performance  will  be  at  chance  (.5  in  this  2AFC 
task).  p(c)l  provdics  a  measure  of  p  (tlK  probability  that  the  target  exceeds  threshold  for  a  given  noise 
level.  Given  this  estimate  of  p,  p(c)2  can  be  used  to  measure  q.  Thus  p(c)l  and  p(c)2  can  be  used  to 
measure  p  and  q  and  thus  predict  performance  for  m=  4  and  m=  16. 

In  contrast  to  the  two-state  threshold  theory.  Green  &  Swels  (1966)  proposed  SDT.  The  decision 
axis  is  continuous  for  an  interval  of  uncertainty  about  whcler  or  not  a  signal  is  present.  The 
representation  of  both  the  signal  and  noise  on  this  decision  axis  is  assumed  to  be  normally  distributed  with 
equal  standard  deviations.  Performance  is  determined  by  the  distance  between  the  means  of  the  signal 
and  noise  distributions  on  the  decision  axis  scaled  in  standard  deviation  units.  This  quantity  is  called  d'. 
The  observer  chooses  a  value  on  the  axis  called  the  criterion.  On  a  trial  the  stimulus  (either  signal  or 
noise)  will  sometimes  be  represented  as  a  value  exceeding  the  criterion  and  the  observer  will  respond 
“signal”.  Those  trials  will  produce  hits  and  false  alarms  for  signal  and  noise  trials  respectively.  If  the 
trial  value  is  below  criterion,  the  observer  responds  “noise"  producing  misses  or  correct  rejections.  For  a 
given  noise  lc%cl,  hit  rate  and  false  alarm  rate  can  be  used  to  measure  d’  and  the  criterion.  As  with  the 
Luce  model  discussed  above,  the  SDT  model  can  be  c.\tcnded  to  deal  with  (he  idcntirication  task  with 
various  m  sizes.  The  criterion  is  assumed  to  lead  to  the  target  1  versus  target  2  decision  as  opposed  to 
signal  versus  noise.  Table  1  in  Swels  (1964)  provides  a  conversion  between  d’  and  probability  of  knowing 
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overall  drop  in  performance  as  a  function  of  noise  level  is  not  prescriptive,  the  relevant  data  to  contrast  tlic 
models  involves  the  drop  in  px:rformance  as  a  function  of  m  size.  As  can  be  seen  in  the  next  two  graphs, 
the  data  for  both  noise  lc\'el=  .125  and  .188  arc  better  fit  by  the  SDT  model  than  by  the  Luce  model  for 
both  the  Geometric  and  T/F  stimuli.  Statistically,  the  superior  predcitions  of  the  SDT/logistic  models 
(relative  to  the  Lucc/median  filler  models)  arc  revealed  in  the  significant  interactions.  The  advantage  of 
the  SDT  model  over  the  Luce  model  interacted  with  level  of  noise  (F(3,51)=4.64;  p<  .006)  and  m  size 
(F(1,I7)=6.83;  p<  .019).  Finally,  the  Noise  Level  X  m  Size  intcraaion  (F(3.51)=3.02;  p<  .04)  and  the 
model  X  Noise  Level  X  m  Size  interaction  (F(3,51)=4.01;  b<  .012)  were  significant.  For  both  stimulus 
types,  the  advantage  of  the  SDT  model  was  greater  at  low  noise  levels  than  at  .3 13  and  was  larger  for  m= 
16  than  for  m=  4. 


Discussion 

In  summary,  the  use  of  a  median  filter  in  the  human  visual  S)’Slem  predicts  a  discontinuity  at  the 
break  point.  No  such  discontinuity  was  found.  Tlic  lack  of  such  a  discontinuity  in  such  an  overly 
simplistic  analysis  can  be  c.vplaincd  away  as  a  result  of  error  variance  in  the  data.  Thai  is  why  m  size  was 
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included  as  an  indepcndcnl  variable.  Where  there  is  a  discriminable  difTercncc.  these  data  arc  belter  f,t  bv 
a  continuous  modeI  (such  as  the  logistic  and  SDT)  than  by  a  discontinuous  model  (median  filterand  ^ 
threshold  theory).  For  these  stimuli,  the  human  \isual  system  docs  not  use  a  median  filter. 

This  inclusion  must  be  qualified  by  Mvo  limitations  in  these  e.x-periments.  First,  the  sUmuli  were 
inar>,  since  this  was  the  simplest  way  to  implement  impulsive  noise.  Retinal  images  are  coded  in  neural 
inng  rate  uhich  roughly  vanes  between  0  and  sescral  hundred  spikes  per  second.  It  is  possible  that  the 
^ta  might  have  ^n  different  if  the  images  had  been  grey-  scale  ratherVhan  binar>..  Secondly  I  is 

^  Produce  different  results.  Intravitreal  blood  might  be  better 

odcled  by  a  smoothing  operator  or  by  a  normal  distribution  of  impulsive  noise.  Endcnce  regarding  the 
importance  of  these  two  possible  limitations  awaits  future  research.  ^ 

conditinn^cy'"^  potential  limitations  of  this  research,  it  should  be  noted  that  the  data  from  these 

rndp  rpf  1  u  ^  ™SSOst  that  the  human  \isual  sj-stem  uses  a  continuous  function  (albeit  digitized)  to 
code  reunal  illuminance  into  neural  firing  rate.  Based  on  PDF  modeling,  the  best  asLmptionIs  that  this 
ffanslation  is  adeqi^tcly  fit  by  the  logistic.  These  data  provide  no  cv'idcncc  that  the  human  visual  system 

Z  Z"  ^  processing  advarl  h 

as  rfuJyr  '"^  f ^  biological  system  and  it  predicts  a  faster  drop  off  in  performance 
as  a  function  of  m  size  than  is  actually  found  vith  binary  images.  penormance 
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ABSTRACT: 


Purpose.  Studies  were  conducted  to  assess  the  effects  of  a  simulated  retinal  hemorrhage  on  the 
level  of  photopigment  bleaching  in  light  adapted  bovine  eyes.  Because  visual  sensitivity  relies 
on  photopigments  in  the  retina,  and  the  level  of  luminance,  the  effect  of  retinal  hemorrhage  on 
visual  sensitivity  is  inOuenccd  by  the  level  of  bleached  photopigment.  Methods.  Bovine  eyes 
were  obLimed  from  a  local  purveyor,  and  dark  adapted  by  storage  on  ice  in  a  light  tight  box, 
for  at  least  one  hour.  Whole  bovine  blood  was  injected  into  the  eyes  to  occlude  an  area  of  the 
retina  just  antenor  to  the  central  retina,  to  simulate  retinal  hemorrhage.  Parameters  such  as 
exposure  time,  light  intensity,  and  blood  volume,  were  varied  to  determine  optimal  conditions, 
nxjxvsurc  for  10  min  with  a  150W  .standard  light  bulb,  placed  100  cm  away  from  eyes  injected 
with  1ml  of  blood,  yielded  the  maximal  reduction  of  bleached  photopigment.  FoUowing 
bleaching,  retinals  were  extracted  by  the  formaldehyde  method,  and  the  U-cis-  and  ^lUrans- 
retinal  were  analyzed  by  high  performance  liquid  chromatography  (HPLC).  Results  Dark 
adapted  control  eyes  yielded  an  average  of  6%  all-.ans-rebn^  (i.e.  94%  I  l-c/s-retm^.  dte 
unbleached  pigment  chromophore).  Light  adaptation  increased  the  bleached  pigment  proportion 
10  about  80%.  The  addition  of  blood  (to  simulate  retinal  hemorrhage)  reduced  the  bleached 
photop.gment  proportion  by  5-12%.  Co.ichM.  Simulated  retinal  hemorrhage  resulted  in  a 
n  in  the  proportion  of  bleached  pholopigment  in  the  eye,  suggesting  that  the  level  of  light 
arriving  at  the  retina  is  reduced  by  blood  located  anterior  to  (he  retina.  The  ability  of  patients 
ith  retinal  hemorrhage  to  detect  light  due  to  a  combination  of  reduced  bleached  pigment  (or 
increased  unbleached  pigment)  and  reduced  luminance  will  be  explored  further. 
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INTRODUCTION: 


Clinical  manifestations  from  a  retinal  hemorrhage  resulting  from  injury  to  the  eye,  such 
as  during  laser  surgery,  would  be  of  concern  to  many  in  the  medical  field.  Simulation  of  retinal 
hemorrhage,  by  occlusion  of  light  to  an  area  of  the  central  retina  by  blood  and  tissue  fluid 
exudate,  would  attenuate  luminance  of  photoreceptors.  Tliis  may  result  in  a  decrease  in  visual 
sensitivity,  or  conversely,  an  increase  in  the  threshold  of  visual  sensitivity.  The  threshold  of 
visual  sensitivity  is  a  log  unit  measurement,  where  normal  vision  has  been  arbitrarily  designated 
as  1  log  unit,  (log  threshold  =  0),  and  impaired  vision  is  measured  in  log  increment  increases 
up  to  about  5  log  units  in  rats  and  18-20  log  units  in  man  (Dowling,  1960,  1963,  Dowling  .S' 
Wald,  1960,  and  Rushton,  19611)).  Visual  threshold  is  defined  as  tlic  lowest  luminance  of  a 
flash  of  light  (1/50-sccond)  necessary  to  evoke  a  barely  measurable  electroretinogram  (ERG) 
(Dowhng  &  Wald,  1960).  As  rhodopsin  content  decreases,  the  visual  sensitivity  threshold 
increases  (Dowling,  I960.  1963,  Dowling  .fcWald,  1960),  so  that  at  a  log  unit  of  3.3-3.5  (in 
rats)  only  1-5%  of  the  normal  rhodopsin  concentration  can  be  extracted  (Dowling  &  Wald). 
Rushton  (1961b)  found  that  in  humans,  for  each  4%  of  rhodopsin  that  is  bleached,  there  is  a  I 
log  unit  rise  in  visual  threshold.  Dowling  &  Wald  (1960)  concluded  that  the  increase  in  visual 
threshold,  could  be  completely  accounted  for  by  the  decrease  in  rhodopsin  concentration  in  rats; 
however,  it  has  since  been  found  that  although  there  is  a  correlaUon,  other  factors  must  be 
considered,(  Rushton,  1961a,  and  Dowling  1960)  such  as  reorganizaUon  of  nerve  synapses 
(Rushton,  1961a),  and  instances  of  photochemical  adaptation,  where  only  small  amounts  of 

bleaching  occur,  yet  human  visual  tlireshold  is  raised  2-3  log  units  instead  of  the  projected  0.5 
log  units  (Dowling,  1963). 
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Visual  sensitivity  is  lowered  after  light  adaptation,  and  is  recovered  during  dark 
adaptation  (Dowling,  1960,  1963).  The  dark  adapted  eye,  contains  1  l-cw-retinaJdehyde,  and 
when  the  eye  is  exposed  to  light,  the  retinal  chromophore  undergoes  conformational  changes  to 
become  all-rra/u-retinaldehydc  (retinal)  (Bongiomo,  1991).  This  principle  was  applied  to  the 
present  study,  where  changes  in  visual  sensitivity  due  to  simulated  retinal'  hemorrhage,  were 
assessed  by  extracting  retinals  from  light  adapted,  hemorrhaged  eyes  for  HPLC  analysis  of  11- 
cis-  aiid  reLinals. 

To  dale  most  of  ihis  type  of  work  is  performed  with  rals,  and,  when  possible,  on 

humans.  Difr,cull,cs  arise  in  me.isiiremenl  of  visual  pigment  concentrations  in  man  (Dowling, 

I960,  1963,  Dowling  &  Wald,  1960,  and  Rushton,  I961a&b).  Ncvertlieless.  some  studies  have 

been  conducted,  such  as  the  study  by  Rushton  (196ib),  with  a  human  monochromal,  or  cone 

vis, on  dencient.  subject.  The  study  cnploycd  retinal  densitometty.  in  which  visual  pign.ents 

■n  a  l.vmg  eye  can  be  measured  (Rushton,  I96la&b).  Brieny,  thodopsin  density  is  measured 

during  bleaching  and  regeneration  events,  by  calculating  relative  differences  in  red  and  green 

light  paths  which  arc  displaced  by  the  density  of  rhodopsin  through  which  the  lights  are  shone 

(Rushton,  1961a),  However,  absolute  rhodopsin  concentrations  cannot  easUy  be  obtained  in 
humans. 

The  literature  also  reflects  little  dau  on  vitamin  A  content  in  bleaching  adaptation  of 
eyes  (Bongiomo,  1991),  however,  since  tiiis  laboratory  studies  retinoid  metabolism  in 
bovine  eyes,  Ute  present  study  was  cond.tcted  to  measure  changes  in  visual  sensitivity  with 
emulated  reUn^  hemorrhage  in  bovtne  eyes.  Visu^  sensirivity  can  be  stated  as  a  funcUon  of 
bleached  photopigmenl  in  the  retina  as  determined  by  the  following  relationship;  log  I/I,  = 
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3.6[f}'-_yj/yj.  where  I  —  threshold  of  visual  sensitivity  at  different  stages  of  bleaching  (Dowling, 
1960,  Rushton,  1961a).  In  the  present  study,  we  took  y,  and  y  to  represent  rhodopsin  content 
in  light  adapted  bovine  retina  in  the  presence  and  absence  of  exogenous  blood,  respectively. 
Retinal  hemorrhage  was  simulated  by  injection  of  a  small  volume  of  blood  into  the  area  Just 

anterior  to  the  retina  of  the  enucleated  bovine  eyes.  Parameters  were  then  varied  to  optimize 
conditions. 

METHODS: 

The  present  study  was  conducted  in  order  that  the  relationship  between  photopigment 
conccntrauon  and  simulated  reunal  hemorrhage  could  be  assessed.  Fre^h  bovine  eyes  were 
obtained  from  a  local  purveyor  and  transported  to  the  laboratory  on  ice  in  a  light  tight  box, 
thereby  inducing  ,he  dark  adapted  stale.  A  total  of  15  bovine  eyes  were  used  in  the  present 
study.  ReUnaJ  hemorrhage  was  simulated  by  the  injection  of  whole  bovine  blood  (100  m1.  500 
Ml,  and  1  ml)  directly  in  front  of  the  central  retina.  The  eyes  were  then  exposed  to  standard 
white  light,  either  150  or  75  watts,  50  or  100  cm  away.  Exposure  time  was  varied  between 
5,  10,  and  20  min  .  The  maximal  effect  was  seen  at  10  min  of  exposure  with  150  watts,  100 
cm  away  from  eyes  injected  with  1  ml  of  blood.  Degree  of  dark  adaption  was  checked  by 
injecting  some  eyes  with  blood  but  not  exposing  them  to  light,  rather  they  were  processed 
omplelely  under  dim  red  light.  Each  retina  was  dissected  under  appropriate  lighting  conditions 
and  retinaldehydes  were  isolated  via  formaldehyde  exttacuon  (Suzuki,  1986).  Retinas  were 
homogenized,  glass/teOcn,  with  4ml  of  O.IM  phosphate  buffer,  pH  6.8,  and  incubated  for  5 
min.  Two  milliliters  of  37%  formaldehyde  was  added,  the  mixture  vonexed,  and  incubated  at 
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roc.  ror  5  .in,  ^,e„  2., 

m  temperature  for  5  min.  Eight  mtllilitcrs  of  dislilled  water  and  12ml  of  n-hexane  were 
then  added.  Tubes  were  vigorously  shahen  and  centrifuged  a,  2«»  rpm,  for  5  min,  at  4-C.  The 
hexane  extractron  was  repeated  twice  and  extracts  were  combined,  evaporated  with  ni^ogen  and 
rexuspended  in  .ml  of  n-hexane.  Rednaldehydes  from  each  sample  were  e.ut'ed  from  .  alumina 

m  the  study  were  resuspended  lOOp.  of  3%  D/H,  and  40p.  were  mtec.ed  into  a  Beckman 
HPLC  (23cm,  5pm,  microsorb  silica  column,  flow  mte  2  m./mi„,  mobile  phase  3^.  D/H). 
S.™ples  collected  later  were  resuspended  in  500pl  of  3%  D,H  and  2«)pl  from  ma,  were  injected 
Ibroiigh  a  Beck-man  502  Auto.sampier  into  a  Beckman  HPLC  with  a  photo  diode  army  detector 
Because  samples  were  analysed  by  two  different  me.Cs  using  two  separate  calibmtion  curves 

r  ----ehyde  desenbed  in  each  mbie,  were  reporied  a; 

erences  in  peak  area  and  not  molar  concentration. 


RESULTS: 


The  first  set  ofexpenments 


were  conducted  to  assess  the  effects  of 


,K„  ,■  exogenous  blood  on 

.na  and  to  es.binsh  the  hgh,  adaptation  procedure.  Bovine  eyes  were  irneced  wi.  .oopi 

r  :  "  ^o,  Show.' 

--  aimuialed  re.tna,  hemorrhages  were  conducted  wtth  .ml  of  blc^.  Percen.ges  of 
"""  -  ..ht  adaption  were  then 


Dificrences  in  ^\\-trQ/is-Tc\\n':\\  /'at  l 

(AT-al)  between  dark  adapted  eyes  and  those  exposed 


to 
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a  130W  standard  whits  light  for  10  min  were  determined  quart Utatjvely.  Retinas  extracted  from 
dark  adapted  eyes  contained  very  lilUe  AT-aJ  (see  Fig.  I),  whereas  the  light  adapted  retina 
contained  predominanUy  AT-al  (see  Fig.  2).  Once  it  was  determined  that  the  dark  and  light 
adaptation  procedures  yielded  expected  results,  the  effects  of  exogenous  blood  on  bleached 
rhodopstn  content  were  tested.  Light  adaptaUon  in  the  presence  of  blood  had  a  bleached 
rhodopsin  content  which  was  slightly  less  than  that  of  the  light  adapted  control  eyes  (see  Fig. 
2  &  3).  Changes  in  rhodopsin  content  simply  due  to  the  presence  of  blood  in  a  dark  adapted 
eye  were  also  tested.  This  treatment  had  the  effect  of  increasing  very  slighUy  the  AT-al  content. 


Subsequent  experiments  were  conducted  in  order  to  optimize  the  bleaching  conditions. 

Extracted  retinaldehydes  were  injected  into  HPLC  and  analyzed  quantiutively.  Results  from 

these  expenments  were  reported  in  molar  concentration  of  AT-al  (bleached  rhodopsin)  relative 

to  11^-al  (unbleached  rhodopsin).  Table  1  shows  the  results  of  dark  adapted  eyes  with 

exogenous  blood  injected.  They  were  exposed  to  I50W  of  standard  white  light  at  a  distance  of 

50  cm.  for  5  and  10  minutes.  Five  minutes  of  exposure  time  did  not  seem  adequate  to 

stgniftcanUy  bleach  the  photopigment,  for  the  difference  in  bleached  rhodopsin  in  experimental 

eyes  (Table  I,  column  II)  from  control  eyes  (Table  1,  column  1)  was  only  a  9%  decrease. 

However,  when  the  eyes  were  exposed  to  150W  of  light  for  10  minutes  this  difference  increased 

>0  22%.  A  third  treatment  was  also  attempted  (data  not  shown)  where  the  eyes  were  exposed 

for  20  min.  However,  the  Ussue  was  presumably  damaged  from  the  hea, -producing  light, 

because  the  reUna  had  deteriorated  so  that  it  was  difficult  to  extract,  hence  retinal  content  could 
not  be  assessed. 


Table  II  shows  the  results  of 


expenments  where  bovine  eyes  were  exposed  to  either 


I50W  of  standard  white  light  100  cm  away,  or  75W  of  standard  white  light  50  cm  away,  for  10 
min.  The  light-adapted,  hemorrhaged  eye  (Table  n,  column  U)  that  was  exposed  to  150W  ® 
100  cm  away,  had  12%  less  bleached  rhodopsin  than  did  the  light  adapted  control  eye  (Table 
n,  column  I).  The  light-adapted,  hemorrhaged  eye  aable  11,  column  n)  exposed  to  75W  ®  50 
cm  away  only  had  a  5%  difference  in  AT-al  when  compared  to  its  control  (Table  n,  column  1). 
This  suggests  that  the  first  treatment  in  this  series  of  expenments  offered  a  pathlength  that  was 
able  to  preserve  more  of  the  1  l-cit.  or  unbleached,  form  of  the  photopigment. 

Results  of  these  studies  taken  together  indicate  that  a  bovine  eyes  injected  with  1ml  of 
blood,  and  exposed  to  150W  of  white  light,  50  cm  away,  for  10  min  had  the  effect  of  inducing 
the  greatest  difference  in  bleached  rhodopsin,  when  compared  to  eyes  without  rednal 
hemorrhage.  In  other  words,  these  conditions  preserve  the  greatest  amount  of  photoreceptor, 
thus  have  the  effect  of  increasing  visual  sensitivity. 


DISCUSSION: 

The  present  study  was  conduced  to  assess  the  effects  of  exogenous  whole  bovine  blood 
in  dark-  and  lighl-adap.ed  bovine  eyes,  as  a  simulation  of  reUnal  hemorrhage.  It  has  been  weU 
established  .ha,  visual  sensitivity  is  a  funedon  of  pho.opigmen.  concenhatton  in  die  redna 
(Dowling  1960,  Rush, on  ,961).  We  adapted  ^e  Rushton-DcwUng  epuadon  for  measuring 
threshold  of  visual  sensidvity  in  „ur  study:  log  1/1.  =  Where  1  =  threshold  of 

visual  sensidvity  at  different  stages  of  bleachmg.  and  y.  and  y  =  Itgh,  adapted  bovine  retina  in 
Ihe  presence  and  absence  of  exogenous  blood,  respectively.  We  found  dia.  our  protocols  for 
tight  and  dark  adaption  were  optimal  a,  exposure  for  10  min  with  a  150W  standard  while  light. 
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placed  50  cm  away  from  eyes  injected  with  1  ml  of  whole  bovine  blood  (22%  difference,  Table 
I).  Such  treatments  resulted  in  a  significant  reduction  in  bleaching  of  photopigment  in  the  retina. 

In  the  first  set  of  experiments,  we  were  able  to  establish  a  baaeline  dark  adapted  control, 
this  sample  contained  only  a  small  percentage  of  all-rro/u-reUnaldehyde. ,  The  light  adapted 
sample  of  this  study  had  a  significantly  higher  all-rrc/u-retinaldehyde  content,  which  was 
decreased  slightly  by  the  addition  of  I  ml  of  blood  (see  Fig.  I,  0,  «&  IIT).  By  repeating  the 
expenments  under  different  conditions,  we  were  able  to  increase  this  difference  in  bleached 
rhodopsin  to  22%  (Table  I  and  II).  The  studies  reported  in  Table  U  included  another 

group  that  was  exposed  for  20  min;  however,  there  was  an  actual  increase  in  AT-al  content  in 
eyes  injected  with  blood.  This  negative  result  might  be  attributable  to  the  effects  of  this  more 
harsh  treatment.  Tlte  light  produced  much  heat,  which  combined  with  the  presence  of  exogenous 
blood  appears  to  cause  a  more  rapid  degradation  of  tissues.  The  retina  became  unattached  to 

the  eye  cup,  and  part  of  it  discarded  with  the  vitreous  humor,  making  it  impossible  to  fuUy 
extract  retinais  from  the  retina. 

Using  the  Rushton-Dowiing  equation,  we  cajculaled  that  treatment  under  optimal 
experimental  conditions  induced  a  decrease  in  visual  sensitivity  threshold  by  0.8  log  units  (or 
a  7.1  times  increase  in  visual  sensitivity). 

A  sample  calculation  using  the  Rushton-Dowiing  equation,  is  shown  below. 

log  I/k,  =  3.6[{y-yJ/yJ 
log  I/Io  =  3.6[(0.93-0.71)/0.93] 
log  I/I„  =  0.8516129  =  0.8 

According  to  this  equation,  the  threshold  of  visual  sensitivity  has  decreased  by  a  factor  of  0.8 
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log  units. 


Results  of  the  present  study  indicate  that  the  addition  of  exogenous  blood  to  the  central 
retina  of  the  eye  also  resulted  in  a  decrease  in  light  luminance  at  the  level  of  the  retina  based 
on  reduced  bleaching  of  photopigment.  However,  the  Rushton-Dowling  equation  suggested  a 
compensatory  increase  in  visual  scnsiuvity  based  on  an  increase  in  the  amount  of  unbleached 
pigment.  Thus  the  net  effects  of  reduced  luminance  and  increased  sensitivity  determines  the 
visual  perception  (such  as  contrast  sensitivity  function)  of  an  individual. 

Some  expenments  have  been  perform  on  humans,  to  determine  how  simulated  retinal 
hemorrhage  affects  Contrast  Sensitivity  Function  (CSF).  Seventeen  dilutions  of  blood  were 
added  to  a  chamber,  placed  in  front  of  the  eyes.  Subjects  were  asked  to  identify  the  clarity  with 
which  they  could  discriminate  objects  at  the  various  blood  concentrations  (Thompson  &  Stoessel, 
1987).  In  this  study,  12.5  ^1  of  whole  blood  was  enough  to  reduce  visual  acuity  to  mere  hand 

motions.  Using  a  Snellen  acuity  chart,  20/20  vision  was  only  achieved  with  1.40  y.\  of  blood 
or  less. 

Theoredcaliy,  ihese  expenments  could  also  be  conducted  in  patients  who  have  retinal 
hemorrhages  due  to  some  natural  causes,  or  by  injury.  In  certain  instances,  where  retinal 
hemorrhage  only  occurs  in  one  eye,  the  other  eye  could  be  used  as  a  control.  Obviously,  their 
retinas  cannot  be  extracted  (unless  the  injury  is  sufnc.enUy  senous  to  warrant  enucleation) 
therefore,  maybe  the  rhodopsin  content  can  be  assessed  with  a  retinal  densitometer,  as  described 
m  the  1961(a)  Rushton  paper.  Future  studies  regarding  measurement  of  contrast  sensitivity 
function  and  light  scallering  due  to  retinal  hemorrhage,  will  provide  more  dau  in  explanation 
of  the  effects  of  retinal  hemorrhage  on  visual  perception. 
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FIGURE  H;  Chromatographic  analysis  of  retinaJdehydes  from  a  light  adapted  bovine  eye. 
Retention  times  for  11-m-  and  all-;ra/w-  relinaldehyde,  5.86  min  and  7.41  min,  respectively. 
See  methods  for  H?LC  conditions. 


Absorbance  @  340nm 


Percentages  of  Bleached  Rhodopsin 
(AE-^mni’-Retinaldeliyde) 
from  Eyes  Light  Adapted  for  5  and  10  min 


Time 

Light  Adapted 

Light  Adapted 

%  Difference 

Without  Blood 

With  Blood 

_ 0) _ 

rin 

5  min: 

77% 

68% 

Q 

10  mm: 

93% 

71% 

y 

22 

Tnd  Pi^molar  concentrations  between  11-ds- 

white  lieht  nia  bovine  eyes  were  exposed  to  150  watts  of 

jTr  f  ^  5  and  10  min.  Experimental  eyes  received  1 

ml  of  whole  bovine  blood.  The  10  min  exposure  time  with  I  ml  of  blood,  allowed 
r  the  least  amount  of  bleaching,  thus  blocking  the  most  light  from  the  retina 
Reunals  were  processed  as  described  in  Methods,  and  analyzed  by  HplcQScm 

bZZTIZT  "-"a 


TABLE  n 

Percentages  of  Bleached  Rhodopsin 
(AlL/'ra/z^-Retinaldehyde) 
at  Two  Different  Pathlengths 


Pathlength 

IJght  Adapted 
Without  Blood 

Light  Adapted 
With  Blood 

%  Difference 

(I) _ 

(ID 

150W  @100  cm: 

82% 

70% 

12 

75W  @  50cm: 

82% 

87% 

5 

Studies  were  conducted  with  variation  of  paihlength.  The  table  shows  the 
difference  in  relative  peak  areas  between  11-cis-  and  all-trans-re  final  when  dark 
adapted  bovine  eyes  were  exposed  to  either  150W  of  standard  white  light  100  cm 
away,  or  75  W  of  standard  white  light  50  cm  away,  for  10  min.  In  this  study,  the 
first  condition  allowed  for  the  greatest  reduction  in  bleached  pigment  upon  light 
adaptation,  thus  blocking  more  light  from  arriving  at  the  retina.  Retinols  were 
extracted  as  described  in  Methods,  and  analyzed  by  HPLC  (25cm,  5fim  microsorb 
column,  flow  rate  Iml/min,  mobile  phase  3%  D/IT)  using  a  Beckman 
Autosampler. 
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The  Effect  of  Disease  on  Signal  Processing 
Characteristics  of  the  Human  Visual  Systems 

FINAL  REPORT 

Dmitry  Gokhman 
April  29, 1995 


My  work  on  this  project  included  a  background  study  of  the  visual  system  with  all  co-pi's, 
theoretical  work  and  preparations  for  an  experimental  study  of  signal  decomposition  in 
the  eye  with  Dr.  Jim  Dykes,  and  a  preliminary  investigation  of  available  approaches  to 
the  problem  of  scattering  in  the  eye  due  to  occlusions. 


Background; 

During  the  initial  stages  (^f  the  project  (he  team  made  a  literature  search  and  conducted 
a  background  study  of  the  visual  system.  We  had  a  seminar  series  and  carefully  went 
through  'Spatial  Vision'  by  R.L.  De  Valois  and  K.K.  De  Valois. 

Signal  decomposifion: 

in  tl>e  visual 

ILPEM  Tl  ^  1  r""'  ”  fo'  eventual  incorporation  into 

^a  and  ^  S'-!  ^Srcetnent  with  experime  t  ” 

data  and  computational  fadlity.  '^^pcrimentai 

We  designed  experiments  to  test  hypotheses  related  to  what  kind  of  decomoosiUon 
nrnatt"  ^  '"=’sking  technique  with  patterns 

““Sir— 

We  started  by  preparing  a  pattern  display  software  for  the  running  of  subjects  Also  we 
ftanualiy  generated  several  patterns  to  test  the  display. 

We  wrote  software  for  faster  and  n.orc  systemaHc  generation  of  patterns  (PTRGEN) 
and,  as  the  hrst  step,  generated  patterns  to  ditferenbate  Fourier  d^compos” 
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n  1  laar  wavelet  decomposition.  We  deconijx>sed  patterns  into  Fourier  components  as 
well  as  into  wavelet  components  with  respect  to  two  orthogonal  bases:  D2  (Maar)  and 
D4  (Daubechies).  We  selected  pairs  of  [lalterns  witli  similar  Fourier  content,  yet  very 
different  wavelet  ccK'fficients. 

We  began  work  on  incorporating  Fourier  decomposition  into  pattern  generation  softvv'are 
and  adding  wavelet  decomposition  with  respect  to  various  bases.  The  idea  was  to  make 
it  easier  to  generate  patterns  for  testing  that  would  disUnguish  various  wavelet  bases 
—  a  more  subtle  task  than  Fourier  versus  wavelet.  In  addition,  we  started  looking  into 
decomposition  with  respect  to  nonorthogonal  bases.  Concurrently  we  were  making 
preparations  to  start  the  running  of  subjects. 

After  discussions  at  the  Septe.nber  2, 1994  meeting  of  the  whole  group  and  consultants 
at  CMl  it  was  decided  to  halt  work  on  these  experiments.  In  partiailar.  Dr.  Randy 
Cdickman  said  that  Fourier  decomposition  has  been  abandoned  by  most  researchers  in 
favor  of  wavelet  decomposition.  Further,  we  got  the  impression  from  him  that  it  is  of 
little  comsequence  as  to  what  kind  of  wavelet  basis  is  appropriate. 

ScaHering; 

After  it  was  decided  to  abandon  the  signal  decomposition  project  1  started  lcx)king  into 
the  problem  of  scattering  in  the  eye  due  to  occlusions.  Tliis  was  not  an  area  of  expertise 
for  me,  so  I  looked  into  the  theoretical  literature  on  the  subject  and  consulted  with  several 
people  more  knowledgeable  in  the  field. 

In  November  1994  1  attended  .several  talks  on  scattering  at  the  23rd  Midwest  Differential 
Equations  Conference  at  the  University  of  Oklahoma,  Norman.  One  of  the  talks  was  by 
Dr.  Klibanov  (UNC  Charlotte)  with  whom  I  had  a  lengthy  discussion.  On  December  2, 
1994  I  attended  a  lecture  by  Dr.  Tom  Roberts  (Brtxrks  AlTl)  'Direct  &  Inverse  Scattering 
for  a  Hyperbolic  System  of  Integro-PDEs'.  With  Dr.  Roberts'  kind  assistance  I  compiiled  a 
list  of  available  MIE  software  for  scattering  in  a  sphere.  On  December  19, 1994 1  attended 
a  meeting  with  Dr.  Ray  Aiqilegate  (UH  1SC5A)  were  we  discussed  occlusions  in  eye. 

Deliverable  items: 

□  Pattern  generation  .software  PIRGEN  in  Perl. 

D  Lists  of  available  software  for  Fourier  and  wavelet  decomposition  and  scattering. 
D  Bibliography  reference  lists  on  Fourier  and  wavelet  decomjxKition,  and  scattering. 
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Conclusions: 


D  I  am  not  completely  convinced  that  it  is  not  of  consequence  as  to  what  type  of 
decomposition  occurs  in  the  visual  system.  While  it  is  generally  accepted  that  some 
sort  of  wavelet  decomposition  is  taking  place,  the  guess  as  to  whtit  type  of  basis  is 
involved  is  largely  a  matter  of  empirical  matching  of  basis  functions  to  CRF  curves. 
Many  researches  choose  the  Gabor  wavelet  for  this  reason.  However,  for  computer 
mcxieling  it  may  be  advantageous  to  select  a  basis  with  better  computational  aspects 
which  acceptably  matches  experimental  data. 

D  The  ability  to  compute  scattering  in  the  eye  due  to  occlusions  is  obviously  of 
great  importance  to  mcxieling  the  effects  of  disease  or  injury  on  vision.  This  is  a 
nontrivial  undertaking  considering  that  scattering  is  a  vast  area  of  research.  My 
own  expertise  in  this  area  is  limited  to  numerical  algorithms  that  are  often  used  to 
compute  scattering,  such  as  the  finite  element  method.  I  would  say  that  work  on 
this  project  would  require  significant  partidpation  of  a  spedalist  in  scattering. 
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Appendix  IV 

VEP  Data  Collection  Effort 


Title: 

Final  Report: 
To: 


Date: 


Introduction 

The  visual  evoked  potential  (VEP)  has  been  used  as  an  objective  measure  of  visual 
functional  parameters  such  as  visual  acuity  and  contrast  sensitivity  (Tyler  et  al.,  1979;  Nelson  et 
al.,  1984;  DeBruyn  et  al.,  1986;  Glickman  et  al.,  1991).  Accordingly,  the  interaction  between 
the  VEP  and  the  actual  viewing  or  stimulus  conditions  has  been  examined  in  numerous  studies. 
For  example,  contrast  sensitivity  measured  by  the  lowest  contrast  producing  a  detectable  VEP 
agreed  closely  with  psychophysical  contrast  sensitivity  for  all  spatial  frequencies,  and  the  change 
in  sensitivity  with  change  in  orientation  of  the  gratings  was  the  same  for  contrast  sensitivity  (CS) 
determined  psychophysically  and  by  VEP  measures  (Campbell  and  Maffei,  1970).  Also,  the 
psychophysically  determined  visual  acuity  depended  on  the  target  area  and  luminance,  such  that 
for  a  given  spatial  frequency  target,  there  was  a  critical  target  area  below  which  the  measured 
acuity  depended  on  the  number  of  cycles  (for  a  periodic  grating  stimulus)  in  the  target  aperture 
(Pokorny,  1968).  Similarly,  for  the  overall  modulation  transfer  function,  i.e.  contrast  sensitivity, 
of  the  retina,  the  observed  contrast  sensitivity  for  sine-wave  grating  stimuli  depended  on  the 
number  of  cycles  contained  in  the  stimulus  when  that  number  fell  below  a  critical  value.  Early 
studies  found  that  the  contrast  sensitivity  function,  when  corrected  for  the  number  of  cycles  in 
the  stimulus,  showed  no  maximum  (i.e.  no  low-frequency  roll-off)  (Hoekstra  et  al.,  1974;  Savoy 
&  McCann,  1975).  Later  work,  however,  modified  this  view  by  showing  that  even  when  at  least 
seven  cycles  are  used  at  each  spatial  frequency  (in  order  to  prevent  the  artifact  due  to  too  few 
cycles)  there  was  still  a  low  spatial  frequency  roll-off  in  central  CS,  as  well  as  dependence  on 
retinal  location  (Rovamo  et  al.,  1978;  Kelly,  1984,  1985).  Both  visual  acuity  and  contrast 
sensitivity  decline  with  retinal  eccentricity  (Ludvigh,  1941;  Kelly,  1985).  These  findings  have 
implications  for  the  determination  of  an  overall  retinal  contrast  sensitivity  function;  for  example, 
due  to  technical  limitations,  an  insufficient  number  of  cycles  in  a  low  spatial  frequency  grating 
or  checkerboard  may  be  displayed  on  the  stimulus  device  to  determine  the  response  adequately. 

The  visual  electrophysiological  responses  exhibit  a  dependence  on  target  diameter  and  retinal 
location  in  a  manner  comparable  to  that  of  the  psychophysical  measures.  The  classic  VEP 
"tuning  curve",  i.e.  a  graph  of  VEP  amplitude  with  respect  to  stimulus  spatial  frequency, 
typically  has  been  described  as  having  an  inverted  "U"  shape  (Harter,  1970),  although  the 
decline  in  amplitude  at  low  spatial  frequencies  may  well  have  resulted  from  the  presentation  of  a 
less-than-optimal  number  of  contrast  changes  (cycles)  in  the  stimulus  field  (Katsumi  et  al., 

1988).  Nevertheless,  by  the  use  of  spatially-restricted  stimuli,  e.g.  an  annular  check-pattern  that 
left  the  fovea  unstimulated,  a  larger  contribution  by  the  macular  region  than  by  the  retinal 
periphery  to  the  VEP  amplitude  (Harter,  1970;  Armington,  1968,  1981;  Sakaue  et  al,  1990),  as 
well  as  threshold  (Cannon,  1983),  was  demonstrated  .  The  peak  of  the  spatial  frequency  tuning 
curve  shifted  to  lower  frequencies  in  the  peripheral  retina  (Harter,  1970;  Armington,  1981; 
Sakaue  et  al.,  1990),  consistent  with  lower  visual  acuity  in  the  peripheral  retina.  The  VEP 
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amplitude  dropped  sharply  as  the  size  of  a  central  scotoma  (centered  on  the  fovea)  increased  up 
to  a  diameter  of  about  6  degrees  (Sakaue  et  al.,  1990).  Beyond  this  size,  the  scotoma  exerted 
much  less  effect  on  the  VHP  amplitude.  This  finding  illustrated  the  macular  dominance  of  the 
VHP,  and  also  pointed  to  the  value  of  the  VHP  in  assessing  the  effect  of  macular  lesions  (Sokol, 
1972).  Recent  investigations  found  abnormalities  in  the  transient  and  steady-state  VEPs  in  up  to 
78%  of  retinal  disorders  (mostly  maculopathies)  (Celesia  &  Kaufman,  1985).  This  corroborated 
the  potential  utility  of  the  VEP  in  predicting  acuity  or  contrast  sensitivity  after  discrete  lesions  to 
the  fovea  or  macular  area.  It  may  be  noted  that  the  macular  dominance  of  the  VEP  may  be  due 
to  a  combination  of  cortical  magnification  on  the  central  retinal  projection  (a  physiological  and 
anatomical  phenomenon  affecting  to  visual  function),  and  the  location  of  the  foveal  projection 
on  the  occipital  pole  where  it  is  most  accessible  to  recording  by  scalp  electrodes. 

The  VEP  was  initially  used  in  this  research  program  to  assess  visual  acuity  in  monkeys  after 
focal  laser  lesions  to  the  fovea  and  parafovea  (Schmeisser,  1992;  Glickman  et  al..  In  Press).  Tlie 
VEP  proved  to  be  a  useful  method  for  objectively  tracking  and  estimating  acuity  after  threshold 
and  suprathreshold  laser  lesions.  The  data  from  this  initial  study  was  used  to  validate  a 
concentric  retinal  model  implemented  on  a  microcomputer  (referred  to  as  the  "CMI  model") 
which,  given  the  size  and  location  of  a  retinal  lesion,  predicted  the  lesion’s  effect  on  the  VEP 
amplitude,  and  by  extrapolation,  on  the  resulting  visual  acuity.  The  model's  basic  assumption  is 
that  retinal  acuity  varies  only  with  eccentricity  (as  hypothesized  by  Kelly,  1985),  so  that  by 
adjusting  a  scaling  factor  for  each  concentric  retinal  region,  the  contribution  of  that  region  to  the 
VEP  response  is  brought  into  accordance  with  the  physiological  response,  and  then  used  to 
predict  the  response  to  a  particular  lesion  input  by  the  user.  While  this  approach  has  proven  to 
be  suitable  for  predicting  the  effects  of  foveal  (central  lesions),  additional  VEP  data,  specifically 
obtained  from  more  peripheral  regions,  was  needed  to  validate  predictions  for  other  retinal 
regions.  Moreover,  some  meridional  differences  in  sensitivity  have  been  reported  (Regan  & 
Beverley,  1983),  which  could  be  responsible  for  variability  in  the  VEP  response  not  properly 
accounted  for  by  the  model.  Therefore,  we  collected  VEP  data  from  normal  subjects  using  a 
variety  of  stimulus  fields  covering,  in  different  experimental  runs,  the  central  macula,  annular 
regions  (i.e.  parafovea),  and  superior  and  inferior  hemifields  (as  an  attempt  to  measure 
meridional  sensitivity).  This  dataset  was  then  compared  to  similar  VEP  data  collected  from 
patients  with  selected  retinal  disorders.  In  addition  to  validating  the  CMI  model,  this  study  was 
also  performed  to  demonstrate  the  capability  of  the  model  in  predicting  visual  function  in  the 
presence  of  naturally  occurring  pathology  limiting  function  primarily  in  central  or  peripheral 
regions  of  the  visual  field. 

Methods 

Visual  Stimulation 


The  visual  stimuli  were  displayed  on  a  19"  diagonal,  high-resolution  color  CRT  display  with 
a  resolution  of  1024  by  1024.  The  stimulus  patterns  themselves  were  produced  by  the 
"VisionProbe"  stimulus  generator  (Conceptual  Mind  Works  WisionProbe  Enterprises,  San 
Antonio,  TX).  The  monitor  was  maintained  at  a  viewing  distance  of  1  m  by  an  adjustable  chin- 
head  rest,  so  that  the  full  extent  of  the  visible  raster  on  the  CRT  subtended  about  a  12  by  12 
degree  stimulus  field.  Subjects  viewed  the  stimulus  screen  monocularly  with  the  other  eye 
patched.  The  VEPs  were  elicited  by  sine  wave  luminance  gratings  counterphased  in  a  square 
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wave  temporal  profile  at  7  Hz  (14  reversals/sec).  Space-averaged  luminance  of  the  gratings  was 
about  6  fL  (20.8  cd/m^)  (green  CRT  gun  only).  Surrounding  areas  in  the  room  had  a  luminance 
of  <2  fL.  The  VisionProbe  was  set  to  sweep  the  contrast  of  the  sine  wave  gratings  from  0%  to 
50%  contrast  in  a  12  sec  period,  with  the  contrast  increasing  logarithmically.  Logarithmic 
sweeps  were  used  to  increase  the  "dwell"  lime  at  the  lower  contrasts  around  threshold.  The 
contrast  threshold  was  determined  at  spatial  frequencies  of  0.5, 1,  3, 5,  8,  12,  and  16  cycles  per 
degree  (cpd),  depending  on  the  field  condition.  Table  1  shows  the  stimulus  matrix  used  in  the 
experiments. 

The  grating  patterns  were  presented  in  fields  of  different  sizes  and  configurations  (see  Table 
1).  Different  field  configurations  were  produced  by  masking  the  full  12  degree  field  with  opaque 
templates  (cf.  pictorial  representation  of  different  stimuli  in  fig.  1).  During  the  recording  epoch, 
the  subject  was  instructed  to  look  at  a  small  red  fixation  spot  centered  in  the  stimulus  field.  The 
chin-head  rest  also  assisted  the  subject  in  maintaining  stable  fixation,  as  well  as  reducing 
artifacts  arising  from  head  or  neck  movements  and  muscle  contractions. 


Table  I.  Stimulus  Matrix 


Field  Type 

Spatial  Frequencies 

Contrast  Range 

Subjects  Tested 

Whole  Field  ( 12x12“) 

0.5,  3,8, 12, 16’ 

0-50% 

Normals  and  Patients 

3“  Central  Field 

0.5,  3,8,12,16’ 

0-50% 

Normals  and  Patients 

6“  Central  Field 

0.5,  3,8,  12, 16’ 

0-50% 

Normals  and  Patients 

6“- 12“  Annulus 

1,  5,  8, 12’ 

0-50% 

Normals  and  Patients 

6“-12“  Upper 

Hemifield 

1,  5,  8,  12’ 

0-50% 

Normals  only 

6“-12“  Lower 

Hemifield 

1,  5,  8,  12’ 

0-50% 

Normals  only 

’These  spatial  frequencies  were  only  tested  if  conditions  warranted,  i.e.  a  response  was  observed. 
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Figure  1.  Pictorial  representation  of  stimulus  masks  used  in  sweep-VEP  experiment 


"5>weep-VEP''  Recording 

Because  the  VEPs  were  recorded  with  a  dynamic  stimulus,  i.e.  the  contrast  of  the  stimulus 
was  increased  gradually  during  data  acquisition  (i.e.  "swept")  over  the  range  of  0  to  50%,  these 
responses  are  called  "sweep- VEPs".  The  VEPs  were  recorded  with  a  type  E5GH  gold  cup 
electrode  (Grass  Corp.,  Quincy,  MA  -  now  Grass  Division  of  Astro-Med,  Inc.,  West  Warwick, 
RI)  placed  on  the  scalp  2  cm  above  the  inion.  Good  contact  between  electrode  and  scalp  was 
ensured  by  first  cleaning  the  recording  site  with  Redux  Paste  (Hewlett-Packard  Medical  Products 
Group,  Waltham,  MA),  placing  the  recording  electrode  on  the  scalp  with  conductive  electrode 
paste  (Grass  Corp.,  type  EC-2),  and  covering  the  electrode  with  a  piece  of  paper  tape  to  ensure 
stability  of  the  lead  during  recording.  Gold-plated  clip  electrodes  (Grass  type  E34D-G)  on  the 
ears  served  as  reference  (on  the  contralateral  ear  lobe  -  the  subjects  viewed  the  stimulus  screen 
monocularly)  and  as  ground  (on  the  ipsilateral  earlobe).  Electrode  impedance  in  all  cases  was 
less  than  lOK  ohms  and  was  usually  below  5K  ohms  before  VEP  recording  began. 

The  stimulus  presentation  and  data  acquisition  system  is  shown  schematically  in  Fig  2.  The 
electrical  signals  were  amplified  with  a  Grass  Model  P5,  AC-coupled  amplifier  with  bandpass 
set  from  1  to  1000  Hz.  The  raw  VEP  signal  was  then  further  low-pass  filtered  with  a  Wavetek 
model  452  active  filter  with  the  band  limit  set  to  100  Hz.  The  output  of  this  filter  was  fed  into  a 
SR530  lock-in  amplifier  (Stanford  Research  Systems,  Palo  Alto,  CA).  The  lock-in  amplifier 
was  referenced  to  the  stimulus  counterphase  frequency  by  an  appropriate  signal  from  the 
VisionProbe  generator.  The  output  of  the  lock-in  consisted  of  DC  voltages  proportional  to  the 
amplitude  and  phase  of  the  input  signal  whose  frequency  corresponded  to  that  of  the  reference 
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frequency  (cf.  Nelson  et  al,  1984;  Glickman  et  al.,  1991).  The  raw  VHP  signal,  the  reference 
signal  (an  analogue  signal  -  corresponding  to  the  stimulus  contrast  -  also  generated  by  the 
VisionProbe),  and  a  voice  commentary  were  all  recorded  on  standard,  1/4"  audio  tape.  The  VHP 
and  control  signals  were  converted  to  audio  frequencies. by  voltage  to  frequency  conversion 
(models  2D  and  3  FM  converters,  Vetter  Corp,  Rebersburg,  PA).  The  lock-in  amplifier  outputs 
were  also  captured  with  the  "Experimenter's  Workbench"  (DataWave  Corp.,  Longmont,  CO) 
data  acquisition  system,  and  stored  on  computer.  These  computer  files  were  analyzed  to 
determine  the  point  in  the  contrast  sweep  at  which  a  recordable  VEP  occurred. 


Figure  2.  Schematic  of  experimental  layout  used  to  record  sweep-VEP  responses. 


Analysis  of  VEP  data 

The  original  plan  for  this  project  called  for  developing  a  computer  program  for  automatically 
determining  contrast  sensitivity  from  the  sweep-VEP  data  files.  Because  the  data  from  tlie 
normal  subjects  were  obtained  first,  it  was  used  to  test  the  analysis  software  as  it  was  developed. 
The  following  approach  was  adopted  for  the  analysis  algorithm.  Due  to  the  variability  observed 
in  the  subjects'  responses  to  a  given  stimulus,  each  subject's  responses  to  a  fixed,  repeated 
stimulus  were  averaged,  point  by  point,  to  get  a  "typical"  response  for  that  individual  and  a  given 
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stimulus  condition.  Thus  all  of  the  sweeps  in  a  file  (typically  7-10  sweeps)  were  averaged,  and 
the  response  detection  algorithm  was  run  on  the  average  sweep  response  dataset. 

The  response  detection  algorithm  used  linear  regression  to  fit  a  line  to  all  the  points  in  a 
window  of  given  size.  Lines  were  fit  to  all  the  points  of  the  average  sweep  by  using  overlapping 
windows.  Starting  with  the  first  window  at  the  beginning  of  the  sweep,  the  slope  of  the  fitted 
line  was  compared  to  an  arbitrary  threshold  value.  If  more  than  two  consecutive  windows  had 
slopes  above  the  threshold  value,  a  response  was  assumed  to  have  occurred.  The  response  point 
was  determined  by  calculating  the  intersection  of  the  line  in  the  first  response  window,  with  the 
line  in  the  previous  window  (whose  slope  was  below  the  threshold).  The  sample  number  where 
this  intersection  occurred  was  correlated  to  a  stimulus  parameter  reference  signal,  also  contained 
in  the  experimental  records,  the  amplitude  of  which  was  proportional  to  the  display  contrast. 

This  contrast  value,  ranging  from  0  to  50%,  was  taken  to  be  the  contrast  of  a  threshold  stimulus 
for  the  particular  spatial  frequency  and  field  size  of  the  luminance  grating  on  the  display. 

Due  to  noise  spikes  in  the  original  data,  a  robust  form  of  linear  regression  was  used;  the  least 
median  of  squares  (LMS)  regression.  This  was  predicted  to  give  a  better  line  fit  in  the  presence 
of  impulsive  noise  than  ordinary  least  squares  regression.  The  size  of  the  line  fit  window  and 
the  minimum  slope  threshold  were  experimentally  determined  using  sweep- VHP  data  from 
normal  subjects.  As  it  turned  out,  however,  the  algorithm  developed  for  the  normal  subjects' 
dataset,  proved  to  overestimate,  by  a  large  margin,  the  contrast  sensitivity  for  the  patient 
subjects.  Therefore,  the  entire  dataset  collected  in  this  study  was  analyzed  by  visual' inspection 
by  two  observers,  i.e.  the  observers  looked  at  the  individual  and  averaged  VHP  traces  and 
determined  if  and  at  what  contrast  a  VHP  occurred  (See  Results,  below).  This  manual  analysis 
was  performed  in  a  masked  fashion,  i.e.  the  observers  did  not  know  the  stimulus  conditions  for 
the  VEP  record  being  analyzed. 

Statistical  differences  in  the  mean  CS  determined  in  the  various  subject  groups  were  tested 
first  by  analysis  of  variance  followed  by  three  post-hoc  tests  (Least  Significant  Difference  test, 
Tukey  test,  and  Unequal  N  test)  of  the  difference  between  means.  Statistical  tests  were  from  the 
CSS  Statistica  Package.  The  raw  data  from  all  subject  groups  are  tabulated  in  Attachment  2. 

Subjects 

Normals.  Eleven  subjects  with  no  history  of  ophthalmic  disease,  other  than  corrected  refractive 
errors,  made  up  this  group.  They  were  mostly  recruited  from  the  University  community  and  had 
their  visual  acuity  measured,  to  ensure  they  had  acceptable  vision  for  the  study.  Visual  acuity 
was  measured  using  a  Bailey-Lovie  chart  #4.  Subjects  were  required  to  read  the  letters  on  the 
chart  until  they  made  a  total  of  7  errors  or  until  they  misidentified  all  5  letters  on  one  line.  A 
numerical  score  was  calculated  based  on  the  number  of  letters  read  correctly,  and  the  score  was 
converted  into  a  conventional  measure  of  visual  acuity.  Subjects  wore  their  usual  corrective 
eyeglasses  (if  any)  during  the  experimental  session.  They  received  $25  upon  completion  of  tlie 
session.  The  visual  acuity  of  the  normal  subjects  are  shown  in  Table  II. 

Patients.  There  were  fourteen  patients  referred  from  the  UTHSCSA  Eye  Consultants  Clinic,  the 
Brady-Green  Eye  Clinic^  and  the  Audie  Murphy  Veterans  Administration  Hospital  with  a 
diagnosis  of  a  retinal  disorder  such  as  retinitis  pigmentosa  (RP),  cystoid  macular  edema  (CME) 
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with  laser  photocoagulation,  ocular  histoplasmosis  syndrome,  glaucoma  with  documented  field 
loss,  and  one  case  of  accidental  laser  exposure.  In  addition  to  the  clinical  history,  all  patients 
had  their  acuity  tested  on  the  day  of  the  VEP  recording,  and,  if  necessary,  wore  their  usual 
correction.  Exclusion  criteria  for  patients  for  this  study  were  the  presence  of  moderate  to  dense 
cataracts  or  visual  acuity  worse  than  20/100.  Patients  received  a  $25  reward  upon  completion  of 
the  VEP  recording  session.  The  diagnoses  and  visual  acuity  of  the  patient  subjects  are  shown  in 
Table  II. 


Table  II.  Subjects 


Subject  ID 

Visual  Acuity 

Log  M.AR 

Visual  Condition 

BBl 

20/18 

-0.05 

Normal 

JHl 

20/20 

0 

Normal 

CGI 

20/19 

-0.01 

Normal 

DFl 

20/15 

-0.11 

Normal 

RGl 

20/15 

-0.11 

Normal 

BB2 

20/15 

-0.13 

Normal 

VHl 

20/22 

0.05 

Normal 

DRl 

20/14 

-0.15 

Normal 

LVl 

20/20 

0 

Normal 

MMl 

20/16 

-0.11 

Normal 

BPl 

20/16 

-0.09 

Normal 

omniiii 

20/18 

-0.05 

Patient:  RP 

VH2 

20/20 

0 

Patient:  RP  -  peripheral  defect 

NGl 

20/68 

0.53 

Patient:  RP  -  peripheral  defect 

FCl 

20/16 

-0.09 

Patient:  laser  injury  -  macular  defect 

WSl 

20/12 

-0.21 

Patient:  retinal  degeneration  (probable 
RP)  -  peripheral  defect 

20/24 

0.07 

Patient:  Toxoplasmosis,  retinal 
degeneration,  RP  -  peripheral  defect 

JWl 

20/54 

0.43 

Patient:  RP  -  peripheral  defect;  field 
constricted  to  10" 

LBl 

20/98 

0.69 

Patient:  RP  -  peripheral  defect 
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MWl 

20/18 

-0.05 

Patient:  ARMD  with  subretinal  scar 
adjacent  to  fovea  -  macular  defect 

KTl 

20/47 

0.37 

-Patient:  RP  &  Usher's  syndrome;  field 
constricted  to  5° 

HPl 

20/43 

0.33 

Patient;  epiretinal  membranes,  drusen, 
RPE  changes,  &  probable  SRN  - 
macular  defect 

YBl 

20/44 

0.34 

Patient;  parafoveal  telangiectasis  with 
SRN,  laser  ablation  -  macular  defect 

MGl 

20/50 

0.4 

Patient:  mild  diabetic  retinopathy, 

CME,  post  focal  laser  photocoagulation 
-  macular  defect 

CMl 

20/85 

0.65 

Patient:  solar  maculopathy  -  macular 
defect 

Abbreviations; 

ARMD;  Age-Related  Macular  Degeneration 
CME:  Cystoid  Macular  Edema 
RP:  Retinitis  Pigmentosa 
SRN:  Sub-Retinal  Neovascularization 


Psychophysical  Determination  of  Contrast  Sensitivity 

In  order  to  compare  the  electrophysiological  and  psychophysical  estimates  of  contrast 
sensitivity,  we  determined  the  contrast  sensitivity  of  each  of  the  normal  subjects 
psychophysically.  We  measured  monocular  contrast  sensitivity  using  stimulus  characteristics 
closely  matching  those  of  the  6°diameter  targets  producing  pattern  VEPs.  The  screen  size  was 
5.6'  at  the  1.6  m  viewing  distance  which  was  maintained  constant  by  an  adjustable  chin-head 
rest.  The  circular  target  was  produced  by  a  5.6°  circular  aperture  in  a  diffusing  Lucite  mask. 

The  spatially  sinusoidal  gratings  were  modulated  with  a  square  wave  temporal  profile  at  7.5  Hz 
(15  reversals/s).  The  average  luminance  of  the  screen  was  23  cd/m^  (~  195  td).  The  ambient 
luminance  was  kept  dim  (<2cd/m')  to  avoid  loss  of  screen  contrast  due  to  ambient  light  reflected 
off  the  screen.  The  black  and  white  monitor  (15  MHz  non-interlaced  RCA  TC1214  with  a  P4 
phosphor)  had  a  60  Hz  frame  rate  and  displayed  240  lines  vertically  in  22.5  cm.  The  spatial 
modulation  was  synchronous  with  the  raster  modulation. 

We  used  a  two  temporal  interval  forced  choice  adaptive  staircase  technique  to  measure 
contrast  sensitivity.  The  two  observing  intervals  were  denoted  by  tones  occurring  before  each. 

A  blank  screen  was  presented  during  one  of  each  pair  of  observing  intervals.  The  subject's  task 
was  to  indicate  by  a  button  press  which  of  the  two  intervals  contained  the  stimulus.  The  contrast 
decreased  by  0.3  log  unit  in  the  initial  testing  sequence  until  the  subject  made  two  errors.  The 
contrast  decrement  and  increment  size  was  then  changed  to  0.15  log  units.  Two  correct 
responses  caused  a  0.15  log  unit  decrease  in  contrast.  One  incorrect  response  caused  a  0.15  log 
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unit  increase  in  contrast.  These  decision  rules  result  in  an  estimate  of  detection  probability 
■(p(C))  of  71%.  Data  were  collected  until  there  were  7  reversals  of  contrast  change  with  the  0.15 
log  unit  steps.  The  final  sensitivity  was  the  mean  of  the  sensitivity  at  the  7  reversal  points.  The 
duration  of  each  observation  interval  was  0.75  s  with  0.4  s  between  intervals  and  0.5  s  between 
trials.  Tlie  stimuli  were  generated,  the  testing  sequence  was  controlled,  and  the  data  were 
collected  using  a  NeuroScientific  Corporation  8010  System.  The  NeuroScientific  Corporation 
8010  System  is  an  eight  bit  PC  compatible  board  controlled  by  NeuroScientific  software.  The 
board  was  driven  by  a  386  DX33  PC. 


Results 

Computer  vs.  Manual  Analysis  of  the  Sweep-VEP. 

As  described  in  Methods,  above,  the  computerized  algorithm  developed  for  this  project  was 
accepted  based  on  its  performance  with  the  sweep-VEP  data  recorded  from  the  normal  subjects. 
That  the  contrast  sensitivities  returned  by  the  automated  method  were  reasonable  -  for  the 
normal  dataset  only  -  was  demonstrated  by  subsequent  comparison  to  the  results  obtained  from 
manual  analysis.  These  results  are  shown  graphically  in  Figure  3.  For  the  3,  6,  and  12  degree 
field  conditions,  the  CS  determined  by  computer  and  that  measured  manually  either  was 
identical  or  largely  overlapped.  There  was  a  consistent  difference  between  the  manual-  and 
computer-returned  data  in  the  6-12  degree  annulus  condition  (manual  >  computer,  see  Figure  3), 
however,  there  was  at  least  some  overlap  in  the  data  for  all  spatial  frequencies  tested,  i.e.  the 
difference  was  not  large. 

In  contrast,  the  computer  analysis  returned  CS  determinations  from  the  patient  data  that  were 
much  higher  (more  sensitive)  than  were  those  of  the  normals,  in  some  cases  by  a  log  unit  or 
more.  This  was  patently  incorrect,  because  for  many  of  the  patients  with  macular  disease,  there 
were  little  or  no  VHP  responses  observed  with  many  stimulus  conditions.  Therefore,  the  CS 
measurements  should  have  been  lower  compared  to  the  normals.  By  examining  the  patient 
datasets  processed  by  the  computer  algorithm,  we  hypothesized  that  in  the  absence  of  a  robust 
response,  the  algorithm  reported  the  first  noise  variation  in  the  record  as  a  response  inflection, 
thus  returning  an  erroneous  CS  value.  It  was  clear  that  further  work  would  be  required  to  make 
the  algorithm  suitable  for  processing  the  weak  responses  often  obtained  from  the  patients,  e,g.  by 
including  additional  "rules"  to  provide  minimum  response  criteria,  or,  perhaps,  by  increasing  the 
length  of  the  analysis  windows  used  to  detect  excursions  from  baseline.  In  any  event,  without 
the  time  or  resources  to  develop  the  computer  algorithm  further,  we  decided  to  re-analyze  the 
entire  dataset  by  hand,  simply  determining  a  response  inflection  by  visual  inspection  of  the 
dataset. 

Psychophysical  Estimate  of  Contrast  Sensitivity  in  Normal  Subjects 

The  results  of  the  psychophysical  tests  of  contrast  sensitivity  are  shown  in  Table  HI. 

Although  the  intent  had  been  to  determine  CS  at  3  and  8  cpd,  two  of  the  spatial  frequencies  used 
in  the  evoked  potential  studies,  several  of  the  subjects  were  actually  tested  at  1.5  and  6  cpd 
instead  before  the  error  was  noted.  This  error,  however,  permitted  us  to  examine  a  wider  range 
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Figure  3.  Comparison  of  CS  determined  by  computer  and  visual  inspection  (Man.)  for  the 
normal  subjects,  for  the  3“-,  6°-,  12°-central,  and  6-12°- annular,  field  size  condition.  Error  bars 
in  this  and  all  following  figures  indicate  ±1  standard  error  of  the  mean. 
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of  spatial  frequencies,  bracketing  those  used  in  the  VEP  study.  All  the  measurements  are  shown 
in  Table  III 


Table  HI.  Contrast  Sensitivity  of  Normal  Subjects 


Subject  I.D. 

1.5  cpd 

6  cpd 

3  cpd 

8  cpd 

0 

277.24 

107.05 

n.d 

n.d 

JHl 

333.01 

51.05 

n.d. 

n.d. 

CGI 

333.01 

107.05 

n.d. 

n.d. 

DFl 

273.37 

124.13 

n.d. 

n.d. 

RGl 

247.68 

107.05 

n.d. 

n.d. 

BB2 

158.87 

175.34 

n.d. 

n.d. 

VGl 

213.60 

68.67 

n.d. 

n.d. 

DRl 

n.d. 

n.d. 

184.21 

145.98 

LVl 

n.d. 

n.d. 

158.87 

65.36 

ESOillllllllllll^ 

n.d. 

n.d. 

175.34 

6.43 

n.d. 

72.14 

118.16 

n.d. 

Note;  "n.d."  =  not  done 


Comparison  of  Electrophvsiological  and  Psychophysical  Measurements  of  Contrast  Sensitivity. 

In  the  normal  subjects,  over  all  spatial  frequencies  tested,  the  CS  determined  by  the  sweep- 
VEP  was  lower  than  that  determined  psychophysically.  The  difference  was  about  a  log  unit  at 
spatial  frequencies  of  3  cpd  or  less,  and  decreased  to  less  than  one  half  log  unit  at  8  cpd.  A 
comparison  of  the  two  types  of  CS  measurement  are  shown  in  Figure  4  (automated  sweep-VEP 
analysis)  and  in  Figure  5  (manual  sweep-VEP  analysis). 

Comparison  of  CS  in  Normal  and  Patient  Subjects. 

For  this  analysis,  the  results  from  patients  were  separated  into  two  groups:  patients  with 
disorders  of  the  central  retinal  field  (termed  "macular"  subjects,  N=6),  and  those  with  disorders 
of  the  peripheral  retinal  fields  (termed  "peripheral"  subjects,  N=8).  The  CS  measurements  made 
by  sweep-VEP  and  determined  manually  from  the  data  are  shown  in  Figure  6.  Results  for  each 
stimulus  field  size  are  shown  in  separate  insets  in  the  Figure. 
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Figure  4.  Comparison  of  psychophysical  and  automated  sweep-VEP  measures  of  CS  for 


all  field  sizes  tested  (normal  subjects  only). 


Figure  5.  Comparison  of  psychophysical  and  manual  sweep-VEP  measures  of  CS 
determined  for  all  field  sizes  (normal  subjects  only) 
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1  2  Degree  Field 


3  Degree  Field 


spatial  Ffoquoncy  (cpd) 

6  Degree  Field 


Normals 

Macular 

Periph. 


Spatial  Frequency  (cpd) 

6-1 2  Degree  Annulus 


Normals 

Macular 

Periph. 


Figure  6.  Comparison  of  CS  data,  determined  by  manual  analysis  of  the  sweep-VEP,  of  patient 
and  normal  subjects  groups.  Each  inset  shows  the  results  from  the  indicated  stimulus  field  size. 
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As  shown  in  the  four  insets  in  Figure  6,  corresponding  to  the  3"-,  6°-,  12°-central,  and  6-12° 
annular  fields,  the  mean  CS  of  the  "macular"  subjects  was  much  lower  than  the  CS  of  normal 
subjects  in  all  conditions.  By  analysis  of  variance  followed  by  post-hoc  tests  on  differences 
between  the  means,  this  difference  was  significant  in  all  conditions  (p<.05).  The  mean  CS 
measured  in  the  "peripheral"  subjects  did  not  differ  significandy  from  the  normals  or  from  the 
"macular"  subjects  in  any  of  the  central  field  sizes  tested,  in  spite  of  the  apparent  difference 
between  peripherals  and  normals  in  the  12°  field  condition,  although  it  should  be  noted  that  there 
is  low  statistical  power  to  detect  small  differences  because  of  the  small  sample  size  and  large 
variability.  With  the  6-12°-annular  field,  however,  the  differences  between  normal,  macular,  and 
peripheral  patients  were  all  significant  (p<.05). 

Comparison  of  CS  in  upper  and  lower  hemifields. 

In  the  normal  subjects,  the  CS  was  measured  in  upper  and  lower  hemifields,  consisting  of 
half  of  the  6-12°  annular  mask,  oriented  so  that  the  appropriate  field  was  stimulated  (see  Figure 
1).  The  results  are  shown  in  Figure  7.  No  significant  difference  was  found  in  the  CS  measured 
in  these  two  hemifields. 


6-12  Degree  Hemi-annulus 


upper 

Lower 


Spatial  Frequency  (cpd) 


Figure  7.  Comparison  of  upper  and  lower  hemifield  CS  determined  in  the 
nomial  subjects. 


Discussion 

Suitability  of  Sweep- VHP  for  clinical  testing. 

The  sweep-VEP  technique  is  no  more  difficult  to  use  with  human  subjects  in  clinical  or 
laboratory  settings  than  are  other  evoked  potential  recording  techniques,  because  the  actual 
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recording  procedures,  e.g.  electrode  preparation  and  application,  are  essentially  similar  to  other 
(primarily  flash  and  steady-state  pattern  VEP)  recording  procedures.  The  main  difference  is  in 
the  stimulus  presentation,  which  does  require  a  specialized  pattern  generator  (cf.  Glickman  et  al., 
1991).  However,  in  the  usual  clinical  setting,  we  recommend  that  fewer  conditions  be  tested, 
probably  using  only  one  field  size,  for  example  only  the  "full-field"  (12“)  size,  or  perhaps  the  6° 
field  if  only  the  macular  response  is  to  be  tested.  Additional  conditions  may  not  add  useful 
clinical  information,  and  may  cause  the  testing  session  to  be  unduly  protracted.  In  our  hands,  the 
patient  testing  period  was  usually  about  three  hours,  which  is  far  too  long  for  the  typical  clinical 
test.  Indeed,  if  the  goal  is  to  distinguish  macular  from  peripheral  dysfunction,  testing  2  or  3 
spatial  frequencies  with  the  12°  field  will  probably  be  sufficient. 


Contrast  Sensitivity  measured  with  the  sweep-VEP  and  bv  psychophysics. 

In  our  study,  the  CS  returned  by  the  sweep-VEP  was  lower  than  that  measured  by  a  standard 
psychophysical  test.  The  difference  was  about  1  log  unit  or  less.  This  difference  was  not 
unexpected  and  is  probably  due  to  the  difficulty  of  distinguishing  an  electrophysiological 
response  near  the  threshold,  as  well  as  the  possibility  that  the  recorded  electrophysiological 
response  is  not  identical  with  the  response  of  the  "perceptron",  i.e.  the  neural  unit  or  units 
mediating  the  psychophysical  response.  An  additional  factor  is  that  in  all  the  psychophysical 
data  cited  here,  except  for  Morris  et  al.  (1987),  tliresholds  were  estimated  by  an  up-down 
contrast  adjustment  series.  The  VEP  thresholds  and  the  psychophysical  thresholds  of  Morris  and 
Harrison  (1987)  were  estimated  from  only  increasing  contrast  series.  The  difference  between 
these  two  estimates  was  only  0.5  -  0.75  log  units.  In  any  event,  our  electrophysiological 
measures  of  CS  were  similar  to  those  reported  by  Seiple  et  al  (1984),  and  our  psychophysical 
measures  of  CS  were  also  similar  to  published  values  of  psychophysically-determined  CS 
(Kelly,  1977;  Seiple  et  al.,  1984;  Morris  et  al.,  1987).  Several  datasets  from  the  literature  of 
psychophysical  CS  are  graphed  in  Figure  8  along  with  the  psychophysical  and 
electrophysiological  data  obtained  in  the  present  study. 

Sweep-VEP  measures  of  CS  and  retinal  pathology. 

It  has  been  shown  in  several  studies  that  the  VEP  recorded  from  the  occipital  pole  principally 
reflects  activity  in  the  central  6  degrees  of  the  visual  field  (Meredith  &  Celesia,  1982;  Yiannikas 
&  Walsh,  1983).  This  may  be  due  to  (1)  the  projection  in  humans  of  the  fovea  to  the  cortex  on 
the  occipital  pole,  putting  the  electrical  dipole  of  the  foveal  projection  in  closest  proximity  to  a 
recording  electrode  on  the  occipital  scalp,  and  (2)  the  disproportionately  large  representation  of 
the  fovea  in  the  visual  cortex  (related  to  the  greater  cortical  magnification  factor  of  the  central 
visual  field).  Thus,  pathology  impairing  foveal  function  is  expected  to  reduce  the  amplitude  or 
other  parameters  of  the  VEP  recorded  from  the  scalp.  Indeed,  in  the  present  study,  the 
"macular"  subjects  had  the  lowest  CS  at  all  spatial  frequencies  and  in  all  field  sizes  tested.  In 
comparison,  the  "peripheral"  subjects  had  CS  measurements  that  were  similar  or  identical  to  that 
of  the  normal  subjects,  except  when  the  stimulus  field  was  restricted  to  a  6-12°  annulus  centered 
in  the  parafoveal  region.  With  this  restricted  field,  some  reduction  in  the  CS  was  noted  in  the 
"peripheral"  subjects,  and  marked  reduction  was  noted  in  the  "macular"  subjects.  Overall,  the 
measurement  of  CS  by  the  sweep-VEP  was  able  to  distinguish  retinal  pathology  among  the 
different  subject  groups,  both  by  the  relative  robustness  of  the  swept- VEP  response  in  the 
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normals  compared  to  the  patients,  as  well  as  by  the  relative  CS  measured  in  the  various  testing 
condition. 


6  Degree  Field 


Normals 

Mor&Har 

Kelly 

PsychoCMI 

PsychoSeip 


Spatial  Frequency  (cpd) 


Figure  8.  Comparison  of  CS  data  from  the  present  study  with  selected  data  from  the 
literature.  Identification  of  datasets:  "Normals":  CS  measured  with  sweep- VHP  from  the 
present  study;  "Mor&Har":  pyscho physical  CS  from  Morris  et  al.  (1987);  "Kelly"; 
psychophysical  CS  from  Kelly  (1977);  "PsychoCMI":  psychophysical  CS  from  the 
present  study;  and  "PsychoSeip";  psychophysical  CS  from  Seiple  et  al.,  (1984).  All  the 
psychophysical  CS  measures  were  obtained  under  similar  testing  conditions. 


The  marked  CS  loss  in  the  larger  fields  and  annulus  in  the  "macular"  subjects  was  rather 
surprising;  the  clinical  impression  in  these  patients  was  that  there  was  more  or  less  normal  retina 
beyond  5°  eccentricity.  One  possible  explanation  of  their  poor  performance  is  that  these  subjects 
may  have  used  parafoveal  retina  for  fixation,  thus  positioning  their  defective  macular  area  closer 
to  the  annulus.  Small  errors  in  this  maneuver  may  have  resulted  in  part  of  the  annular  field 
being  occluded  by  their  scotomatous  field,  thus  reducing  the  apparent  VHP. 

The  similarity  in  the  CS  of  the  "peripheral"  subjects  compared  to  the  normals  was  expected. 
Most  of  these  subjects  had  retinitis  pigmentosa  with  visual  field  sizes  >10°  in  diameter.  Thus, 
with  central  stimulus  fields  up  to  12°  in  diameter,  these  subjects  were  able  to  use  essentially 
normal  retina  to  observe  the  stimulus  screen.  With  the  annular  field,  however,  the  field 
constriction  of  some  of  the  subjects  may  have  started  to  encroach  on  the  visual  subtense  of  the 
stimulus  screen,  resulting  in  loss  of  VHP.  This  reduction  in  VEP  may  account  for  the  lower  CS 
of  the  "peripheral"  subjects  in  the  6-12°  annular  field  test  condition. 
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Using  the  sweep- VEP  data  to  validate  the  CMT  mode.!. 


The  data  from  the  normal  human  subjects  can  be  compared  directly  to  the  output  of  the  CMI 
retinal  contrast  sensitivity  model  in  order  to  assess  how  reasonable  was  the  selection  of  the 
scaling  factors  assigned  to  each  of  the  concentric  retinal  segments  of  the  model  which  adjust  the 
model  output  to  match  the  physiological  data.  For  the  patient  data,  the  site  and  extent  of  the 
lesions  may  be  inferred  from  their  clinical  data  (see  attachment  1),  and  this  lesion  information 
entered  into  the  model.  Then,  the  subject's  measured  CS  data,  as  determined  from  the  sweep- 
VEP  data,  may  be  compared  to  the  CS  (at  any  selected  retinal  location)  predicted  by  the  CMI 
model.  If  discrepancies  occur  between  the  model  and  the  experimental  data,  they  may  have  been 
due  to  a  mismatch  between  the  physical  extent  of  the  "observed"  lesion  and  the  actual  extent  of 
the  pathological  process,  or  to  larger  than  expected  effects  of  retinal  lesions  upon  the  VHP. 

Some  adjustment  of  the  scaling  factors  to  reduce  any  discrepancies  may  have  to  be  considered. 
Nevertheless,  these  comparisons  will  help  improve  the  application  of  the  CMI  model  to  the 
prediction  of  visual  function  in  the  presence  of  various  stressors  or  retinal  lesions,  as  well  as  its 
anticipated  use  in  visualizing  visual  disorders  as  an  aid  to  educating  patients  and  their  families. 
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Attachment  2 
(Appendix  B) 

Spatial  Frequency  (cpd) 


0.5 

3 

8 

12 

16 

Normals:  1 

2  Degree  Field,  Log  Contrast  Sensitivity 

BB1 

1 .3872 

1 .5376 

2.1549 

BB2 

1 .3872 

1 .4202 

1 .2757 

0.9245 

0.6536 

BP1 

1 .3979 

1.4815 

1.0410 

0.8386 

1 .6990 

CG1 

1 .3872 

1 .3279 

0.6615 

0.5654 

DF1 

1 .2291 

1 .5229 

0.9208 

0.5072 

DR1 

1 .9586 

1.4815 

1.1487 

0.6840 

0.7825 

JH1 

0.5638 

0.5513 

0.4949 

LV1 

0.8827 

1 .4559 

1 .0269 

0.4522 

0.6596 

MM1 

1 .9208 

2.1549 

3.0000 

1 .5850 

RG1 

1.1871 

0.9626 

0.4711 

1 .4685 

VH1 

1 .2366 

1.5086 

1 .0044 

1 .0555 

1 .5850 

Avg 

1.3583 

1.4844 

1.1741 

0.9050 

1.1160 

Std 

0.2789 

0.1770 

0.4997 

0.7251 

0.4774 

s.e.m. 

0.1054 

0.0590 

0.1580 

0.2417 

0.1804 

Spatial  Frequency  (cpd) 

0.5 

3 

8 

12 

16 

Normals: 

3  Degree  Field,  Log  Contrast  Sensitivity 

BB1 

1 .9208 

2.2218 

1 .3872 

BB2 

0.8601 

1 .2757 

0.7799 

0.4330 

BP1 

2.3010 

1 .4437 

1 .9208 

1 .0223 

CGI 

1 .4685 

0.9586 

0.6757 

0.6716 

DF1 

2.0969 

0.8386 

0.8894 

0.7258 

DR1 

0.6737 

1 .6021 

1 .3768 

0.6478 

0.7423 

JH1 

0.6038 

0.6478 

0.5467 

LV1 

0.9431 

0.5884 

0.4949 

0.5058 

0.7212 

MM1 

RG1 

0.5229 

1 .6778 

0.5935 

1 .2676 

VH1 

0.9547 

1 .0086 

1 .6021 

0.5735 

Avg 

1 .2656 

1 .2846 

0.9509 

0.9157 

0.6175 

Std 

0.6534 

0.4782 

0.4517 

0.3631 

0.1248 

s.e.m. 

0.2310 

0.1691 

0.1707 

0.1210 

0.0721 

B-1 


Spatial  Frequency  (cpd) 


0.5 

3 

8 

12 

16 

Normals:  6  Degree  Reid,  Log  Contrast  Sensitivity 

BB1 

1 .2007 

0.5686 

0.8356 

0.7852 

BB2 

0.6946 

1.4815 

0.7620 

1.7447 

0.6799 

BP1 

0.9318 

1 .3665 

1 .4949 

2.0000 

2.0000 

CGI 

0.7447 

0.6478 

1.2291 

DF1 

0.8697 

1 .3279 

2.0000 

2.3979 

DR1 

0.8416 

1 .5686 

1 .0757 

0.5560 

0.6308 

JH1 

0.5638 

0.6676 

0.4522 

LV1 

1 .6021 

0.4868 

0.5817 

1.0177 

MM1 

2.0000 

1.1367 

0.7932 

RG1 

0.9706 

1 .5229 

0.4389 

1 ,0458 

0.7077 

VH1 

1 .8239 

1 .0757 

0.9393 

2.5229 

Avg 

1.1499 

1.1461 

0.9024 

1 .4293 

1 .0072 

Std 

0.4678 

0.3526 

0.4551 

0.7182 

0.5146 

s.e.m. 

0.1559 

0.1175 

0.1439 

0.2539 

0.2573 

Spatial  Frequency  (cpd) 

1 

5 
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Normals:  6-12  Degree  Annulus,  Log  Contrast  Sensitivit 
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Normals;  6-12  Degree  Upper  Hemifield,  Log  Contrast  Sensitivity 
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Normals:  6-12  Degree  Lower  Hemifield,  Log  Contrast  Sensitivity 
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Spatial  Frequency  (cpd) 

0.5 

3 

8 

Peripherals:  12  Degree  Field,  Log  Contrast 
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Peripherals:  3  Degree  Field,  Log  Contrast  Sensitivity 
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Peripherals:  6  Degree  Field,  Log  Contrast  Sensitivity 
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Spatial  Frequency  (cpd) 
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Peripherals:  6-12  Degree  Annulus,  Log  Contrast  Sensitivity 
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Maculars: 
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Log  Contrast  Sensitivity 

0.0000 

0.5638 

0.6308 

0.0000 

0.0000 

0.0000 

0.6440 

1 .3768 

0.0000 

0.0000 

1.3188 

1.3468 

0.4323 

0.5479 

0.4883 

0.6091 

0.2184 

0.2724 

12 

1 .3979 
0.0000 
0.0000 
0.0000 
0.0000 
1 .0362 
0.4057 
0.5832 
0.2608 


Spatial  Frequency  (cpd) 


0.5 


Maculars:  3  Degree  Field 
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Maculars:  6  Degree  Reid,  Log  Contrast  Sensitivity 
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Maculars:  6-1 2 

Degree 

Annulus,  Log  Contrast  Sensitivity 
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Appendix  V 

The  Visualization  Tool  User’s  Guide 

(withdrawn  from  published  report) 


Appendix  VI 


The  Inhomogenous  Retinal  Cortical  Model  User 

Guide 

(withdrawn  from  published  report) 


Appendix  VII 

The  VisionProbe™  User’s  Guide 

(withdrawn  from  published  report) 


